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RESEARCH  SUMMARY 

The  National  Fire-Danger  Rating  System  (NFDRS), 
implemented  in  1972,  was  revised  during  a  3-year  pro- 
ject (1975  to  1978)  and  reissued  as  the  1978  NFDRS. 
This  report  describes  the  developmental  history  of  the 
NFDRS  and  its  technical  foundation. 

Detailed  information  is  provided  on  modeling  forest 
fuels  and  fuel  moisture,  and  on  development  of  the 
NFDRS  components  and  indexes.  The  report  presents 
equations  used  in  the  1978  NFDRS  and  an  extensive 
bibliography. 
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HISTORY  AND  BACKGROUND 

Fire  control  conferences  called  by  the  Forest  Service, 
U.S.  Department  of  Agriculture,  in  Ogden,  Utah,  in  1940 
and  1954  highlighted  the  need  for  a  uniform  fire-danger 
rating  system  that  could  be  applied  nationwide.  Con- 
ference committees  recommended  that  the  system  focus 
on  the  environmental  factors  controlling  the  moisture 
content  of  fuels. 

In  1954,  several  different  fire-danger  rating  systems 
were  in  use  across  the  Nation.  Improved  communica- 
tions and  transportation,  however,  were  making  possible 
mutual  assistance  agreements  among  fire  control 
organizations.  State  compacts  and  agreements  among 
Federal  agencies  and  regions  brought  together  fire  con- 
trol teams  from  widely  dispersed  sections  of  the  country. 
A  uniform  system  of  rating  fire  danger  and  fire  behavior 
was  therefore  essential  for  efficient  communication 
among  all  those  concerned  with  wildland  fires. 

In  1958,  a  committee  of  Forest  Service  fire  research 
and  wildfire  control  personnel  decided  that  development 
of  a  national  fire-danger  rating  system  was  feasible.  In 
June,  the  Washington  Office  Division  of  Fire  Research 
organized  a  team  headed  by  John  Keetch,  Washington 
Office,  Aviation  and  Fire  Management,  to  formulate  and 
develop  the  system.  FuU-time  work  began  a  year  later. 

By  1961,  the  basic  structure  of  a  four-phase  fire-danger 
rating  system  had  been  outlined,  but  only  the  first 
phase,  fire  spread,  was  ready  for  field  testing.  The 
spread  phase  provided  two  indexes  to  predict  the 
relative  forward  spread  of  a  fire— one  for  fires  burning  in 
a  comparatively  closed  environment  under  a  timber 
canopy  (timber  spread  index),  and  the  other  for  fires 
burning  in  open  areas  of  fine  fuels.  A  third  number,  the 
buildup  index  (BUI),  was  designed  to  indicate  the 
cumulative  drying  of  the  heavier  fuels.  The  BUI  was 
used  in  the  computation  of  the  timber  spread  index. 
Following  field  testing  in  1962  and  1963,  the  Forest 
Service  issued  a  handbook  in  1964  (FSH  5109.11)  cover- 
ing the  spread  phase  of  the  planned  development.  By  the 
next  year,  most  fire  control  organizations  in  the  United 
States  were  using  at  least  a  modified  version  of  the 
spread  index. 

The  regional  adaptations  that  quickly  followed  in- 
dicated that  the  spread  index  was  not  uniformly  ap- 


plicable across  the  country.  Furthermore,  because  the 
remaining  phases— ignition,  risk,  and  fuel  energy— were 
not  available,  many  fire  control  agencies  failed  to  adopt 
the  new  system,  preferring  instead  to  continue  using  in- 
place  fire-danger  methods.  The  Keetch  project  was 
closed  before  these  three  phases  were  developed. 

In  1965,  a  research  unit  at  Seattle,  Wash.,  led  by 
Donald  F.  Flora,  took  a  new  look  at  the  needs  and  re- 
quirements of  a  national  system.  A  research  forester 
working  for  Flora,  James  E.  Hefner,  surveyed  fire  con- 
trol agencies  throughout  the  country,  analyzed  their 
requirements,  and  recommended  that  research  leading  to 
the  completion  of  the  National  Fire-Danger  Rating 
System  be  resumed. 

In  1968,  the  Forest  Service  established  the  National 
Fire-Danger  Rating  System  research  work  unit  at 
Fort  Collins,  Colo,  (appendix  A).  Led  by  research 
meteorologist  Mark  J.  Schroeder,  the  unit  formulated 
the  following  goals: 

1.  A  1972  target  date  was  set  to  complete  develop- 
ment of  the  system  for  field  use.  Fire  researchers  agreed 
that  a  fire-danger  rating  system  superior  to  any  current- 
ly in  use  could  be  developed  from  knowledge  at  hand;  it 
was  not  necessary  to  wait  until  all  pertinent  research 
was  completed. 

2.  The  system  would  be  structured  to  enable  informa- 
tion such  as  better  prediction  equations  and  improved 
fuel  models  to  be  readily  incorporated.  Such  refinements 
would  take  the  form  of  updated  computer  programs  or 
new  tables  supplied  to  users;  the  basic  format  and  defini- 
tions were  to  remain  unchanged. 

3.  The  system  would  be  introduced  as  a  complete, 
comprehensive  package,  not  index  by  index. 

4.  The  complete  system  would  include  a  subjective 
evaluation  of  "risk."  The  development  of  an  objective 
method  would  be  deferred  until  the  physics  of  fuel 
moisture  relationships  and  fire  behavior  had  been 
developed  sufficiently  to  meet  the  needs  of  the  system, 

5.  Ultimately,  the  system  would  be  purely  analytical, 
based  on  the  physics  of  moisture  exchange,  heat 
transfer,  and  other  known  aspects  of  the  problem. 
Although  laboratory  and  field  checks  would  be  made  and 
experimentation  needed  to  establish  basic  relationships, 
the  fire  behavior  based  aspect  of  the  system  would  not 
be  based  on  empirical  studies  or  statistics. 
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6.  The  system  would  be  evaluated  and  updated  by 
1978. 

In  1970,  a  preliminary  version  of  the  system  was 
tested  in  Arizona,  New  Mexico,  and  Georgia.  Eight 
National  Forests,  one  Bureau  of  Land  Management 
(BLM)  district,  two  National  Park  Service  (NPS)  units, 
and  the  Georgia  Forestry  Commission  participated  (see 
appendix  A).  In  1971,  an  improved  version  was  tested  in 
the  Southwest.  The  Forest  Service,  BLM,  NPS,  Bureau 
of  Indian  Affairs  (BIA),  and  State  agencies  conducted 
field  trials  at  nearly  150  locations  in  the  continental 
United  States  and  Alaska. 

The  1972  National  Fire-Danger  Rating  System 
(NFDRS)  (fig.  1)  provided  three  fire  behavior  com- 
ponents and  three  indexes  for  rating  fire  danger.  It  had 
five  fuel  classes  (three  dead,  two  Hve),  nine  fuel  models, 
and  three  slope  classes.  Live  fuel  moisture  was 
estimated  by  measuring  the  ratio  of  green  material  to 


total  "fine"  plant  material  in  ten  1-foot-diameter  ground 
samples  along  a  300-foot  transect  (Fosberg  and 
Schroeder  1971).  Separate  but  crude  risk  factors  for  both 
man-caused  and  lightning-caused  fires  were  provided. 
These  were  combined  with  a  single  ignition  component 
to  provide  an  occurrence  index  (a  number  related  to  the 
potential  fire  incidence  within  a  rating  area).  All  of  the 
components  and  indexes  of  the  1972  NFDRS  were 
normalized  on  a  scale  of  0  to  100. 

Forest  Service  Research  Paper  RM-84  (Deeming  and 
others  1972,  revised  1974)  was  the  summary  publication 
of  the  development  effort.  RM-84  was  written  for  the 
field  user  and  contained  complete  instructions  for  com- 
putation and  application  of  the  components  and  indexes 
for  its  nine  fuel  models.  It  was  nontechnical  and 
qualitative,  but  a  technical  development  paper  was 
prepared  as  an  office  report  (Schroeder  and  others  1973). 
Much  of  the  material  here  is  from  that  report. 
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Figure  1. -Structure  of  the  1972  National  Fire-Danger  Rating  System  (from 
Deeming  and  otfiers  1972,  revised  1974,  p.  3). 
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In  1973,  the  computer  program  FIRDAT  (Furman  and 
Helfman  1973)  was  made  available  at  the  U.S.  Depart- 
ment of  Agriculture's  Fort  Collins  Computer  Center 
(FCCC)  to  process  historical  fire-weather  data  through 
the  NFDRS  algorithms.  In  1975,  the  National  Fire- 
Weather  Data  Library  (NFWDL)  for  archival  and 
retrieval  of  fire-weather  data  was  developed  and  installed 
at  FCCC  (Furman  and  Brink  1975).  FIRDAT  and  the 
NFWDL  provided  a  systematic  means  of  developing  fire- 
weather  and  fire-danger  cHmatologies,  which  are  critical- 
ly important  for  fire  management  planning. 

The  same  year,  after  4  years  of  development  and 
testing,  the  entire  fire-danger  rating  process  was  com- 
puterized in  the  AFFIRMS  processor  (Helfman  and 
others  1975,  revised  1980).  AFFIRMS,  an  interactive 
time-share  program,  computes  fire-danger  ratings  from 
daily  fire-weather  observations  and  forecasts  and  creates 
fire-weather  data  tapes  that  are  incorporated  into  the 
NFWDL  at  FCCC. 

By  the  summer  of  1976,  data  from  more  than  800  fire- 
weather  stations  across  the  United  States  were  being 
processed  through  AFFIRMS,  and  the  system  indexes 
and  components  were  being  calculated  manually  each 
day  for  some  400  stations.  By  the  spring  of  1977,  all 
Federal  agencies  and  35  State  agencies  charged  with 
forest  and  rangeland  fire  protection  responsibilities  were 
using  the  1972  version  of  the  National  Fire-Danger 
Rating  System. 

Responsibility  for  the  1978  NFDRS  update  was  given 
to  John  E.  Deeming  who  had  joined  the  NFDRS  project 
at  Fort  CoUins,  Colo.,  in  1970.  The  project  was  relocated 
to  the  Intermountain  Forest  and  Range  Experiment 
Station,  Northern  Forest  Fire  Laboratory,  Missoula, 
Mont.,  during  the  summer  of  1975.  Deeming  was  joined 
that  summer  by  Robert  E.  Burgan  and  the  following 
spring  by  Jack  D.  Cohen,  both  research  foresters.  The 
1978  update  reflects  changes  made  possible  by  feedback 
from  1972  system  users  and  advances  in  fire  science  and 
fuels  technology. 

The  1978  National  Fire-Danger  Rating  System  update 
addressed  several  major  problems  identified  in  the  1972 
system: 

Response  to  drought.— The  1978  update  has  improved 
response  to  short-term  drought,  with  the  addition  of  a 
1,000-hour  timelag  dead  fuel  class,  and  the  inclusion  of 
live  fuels  in  all  but  the  slash  fuel  models.  Live  fuel 
moisture  models  emulate  plant  moisture  response  to 
phenological  cycles,  replacing  the  fuel  moisture  transects 
used  in  the  1972  NFDRS. 

Seasonal  sensitivity.— Calculation  of  fuel  moisture  for 
the  100-hour  and  1,000-hour  timelag  fuels  was  altered  to 
account  for  a  day's  drying  power  as  affected  by  day 
length.  Day  length  is  calculated  from  the  date  and  a  fire- 
weather  station's  latitude. 

Component  sensitivity.— The  1972  system  had  its  com- 
ponents and  indexes  normalized  to  a  0  to  100  scale, 
which  often  caused  moderate  fire  climates  to  never  see  a 
rating  greater  than  10.  The  spread  and  energy  release 
components,  and  the  burning  index  now  have  open-ended 
scales  that  yield  a  threefold  to  fivefold  increase  in 
sensitivity. 


Occurrence  indexes,— The  1978  NFDRS  has  separate 
occurrence  indexes  for  man-caused  and  lightning-caused 
fires  (the  1972  system  had  only  one  occurrence  index  for 
both  types  of  fires).  The  occurrence  index  models  in  the 
1978  update  reflect  major  improvements  over  the 
simpUstic  model  used  in  the  1972  system. 

Ignition  component.— The  1978  NFDRS  ignition  com- 
ponent is  a  function  of  the  probabihty  of  ignition  and 
spread  component.  The  1972  system's  ignition  compo- 
nent was  simply  the  probability  of  ignition. 

More  fuel  models.-The  1978  NFDRS  offers  20  fuel 
models  to  describe  a  fuel  situation;  the  1972  system 
had  9. 

Better  slope  definition.— The  1978  NFDRS  has  five 
slope  classes  to  describe  an  area's  topography;  the  1972 
system  had  three.  Two  slope  classes  were  added  to  cover 
steep  terrain  (greater  than  50  percent  slope). 

Live  fuel  moisture  models.-The  1978  NFDRS  intro- 
duced models  for  computing  the  fuel  moisture  of  live 
herbaceous  and  shrub  fuels.  The  model  also  causes 
herbaceous  material,  when  cured,  to  be  transferred  to 
the  1-hour  dead  fuel  loading  class  to  better  emulate 
actual  fire  danger  prior  to  spring  green-up  and  after 
autumn  curing. 

A  summary  pubHcation  (Deeming  and  others  1977) 
describes  the  1978  NFDRS  and,  for  NFDRS  users  not 
on  the  AFFIRMS  network,  there  is  a  manual  version  of 
the  1978  system  (Burgan  and  others  1977).  The  1978 
National  Fire-Danger  Rating  System  is  now  being  used 
by  most  Federal,  State,  and  private  agencies  charged 
with  wildland  fire  protection.  The  1978  changes  have 
been  included  in  new  releases  of  FIRDAT  and 
AFFIRMS  (Main  and  others  1982;  Helfman  and  others 
1975,  revised  1980). 

PHILOSOPHY 

The  development  and  application  of  both  the  1972  and 
1978  NFDRS's  are  grounded  in  six  major  principles: 

1.  The  system  would  consider  only  the  "initiating 
fire."  This  is  defined  as  a  fire  that  is  not  behaving 
erratically;  it  is  spreading  without  spotting  through  fuels 
that  are  continuous  with  the  ground  (no  crowning).  The 
"state  of  the  art"  does  not  yet  extend  to  fires  that 
exhibit  erratic  behavior  other  than  to  show  that  extreme 
behavior  is  correlated  with  increasing  fire  danger. 

2.  The  system  would  provide  a  measure  of  that 
portion  of  the  potential  job  of  containment  that  is 
attributable  to  fire  behavior.  The  concept  of  containment 
as  opposed  to  extinguishment  is  essential  because  it 
limits  the  fire  behavior  prediction  to  the  head  of  the  fire. 
Those  portions  of  the  containment  job  dealing  with 
accessibility,  soil  condition,  and  resistance  to  Hne  con- 
struction must  still  be  evaluated  by  other  means. 

3.  The  length  of  the  flames  at  the  head  of  the  fire  was 
assumed  to  be  directly  related  to  the  contribution  of  fire 
behavior  to  the  containment  job. 

4.  The  system  would  attempt  to  evaluate  the  "worst" 
conditions  on  a  rating  area  by  using  meteorological 
measurements  taken  (a)  when  fire  danger  is  normally  the 
highest  (usually  in  tne  early  afternoon),  (b)  at  sites  in  the 
open,  and  (c)  where  possible,  at  sites  on  drier  (southerly 
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or  westerly)  exposures.  This  means  that  extrapolation  of 
fire-danger  values  to  areas  other  than  those  immediately 
in  the  vicinity  of  the  fire-danger  station  would  involve 
scaling  down,  not  up. 

5.  The  system  would  provide  ratings  that  would  be 
physically  interpretable  in  terms  of  fire  occurrence  and 
behavior.  These  evaluations  could  then  be  used  alone  or 
in  combinations,  giving  the  user  the  flexibility  needed  to 
deal  with  the  entire  spectrum  of  fire  control  planning 
problems. 

6.  Ratings  would  be  relative,  not  absolute.  The  ratings 
would  be  linearly  related  to  the  particular  aspect  of  fire 
danger  being  evaluated.  This  means  that  when  a  compo- 
nent or  index  doubles,  a  doubling  of  the  rated  activity 
relative  to  what  has  previously  been  observed  should  be 
anticipated.  Because  of  the  many  variables  in  the  com- 
putations, the  low  spatial  and  temporal  resolution  of 
fuels  and  weather  data,  and  incomplete  understanding  of 
some  relationships,  fire  danger  can  only  be  broadly 
defined  within  a  rating  area. 


STRUCTURE 

The  1978  NFDRS  (fig.  2)  provides  four  indexes  to 
facilitate  the  planning  of  fire  control  activities:  the  man- 
caused  fire  occurrence  index  (MCOI);  the  lightning- 
caused  fire  occurrence  index  (LOI);  the  burning  index 
(BI);  and  the  fire  load  index  (FLI). 

The  MCOI  is  derived  from  man-caused  risk  (RmC)- 
assessment  of  man-caused  fire  sources  in  the  rating  area, 
and  the  ignition  component  (IC),  the  likelihood  that  a 
firebrand  will  cause  a  reportable  fire. 

The  lightning-caused  fire  occurrence  index  (LOI), 
similar  in  concept  to  the  MCOI,  is  derived  from  the  igni- 
tion component  (IC)  and  lightning  risk  (Rl),  an  indicator 
of  thunderstorm  and  lightning  activity.  After  being 
scaled  to  an  area's  experience,  both  the  MCOI  and  LOI 
can  be  used  to  predict,  on  the  average,  the  total 
expected  number  of  reportable  fires  that  will  occur  on  a 
rating  area. 


1978     NATIONAL  FIRE -DANGER  RATING  SYSTEM 
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BASIC  OBSERVATION  TI/VIE 


LIGHTNING  MAX.  &  MIN.  RELATIVE  AIR  CLOUDI-  WIND- 
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Figure  2.— Structure  of  the  1978  National  Fire-Danger  Rating  System  (fronn 
Deeming  and  otiiers  1977,  p.  6). 
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The  burning  index  is  derived  from  the  spread  compo- 
nent (SC),  a  relative  index  of  rate  of  fire  spread,  and  the 
energy  release  component  (ERC),  a  relative  index  of  the 
amount  of  heat  released  per  unit  area  in  the  flaming 
zone  of  an  initiating  fire.  Considered  together,  they 
indicate  the  difficulty  of  containment.  The  combined 
index,  BI,  is  linearly  related  to  the  length  of  flames  at 
the  head  of  the  fire. 

The  difficulty  of  containing  a  single  fire  (the  BI)  and 
the  expected  number  of  fires  projected  by  the  MCOI  and 
LOI  combine  to  produce  the  fire  load  index  (FLI),  a 
measurement  of  the  total  fire  containment  job.  The  FLI 
is  the  NFDRS  cumulative  index,  integrating  risk,  igni- 
tion probability,  and  fire  behavior  potential. 

The  risk  ratings  (RmC        Rl),  man-caused  fire  occur- 
rence index  (MCOI),  lightning-caused  fire  occurrence 
index  (LOI),  ignition  component  (IC),  and  fire  load  index 
(FLI)  are  expressed  on  a  scale  of  0  to  100.  The  scales  of 
the  spread  component  (SC),  energy  release  component 
(ERC),  and  burning  index  (BI),  are  open  ended. 

Eleven  elements  of  a  fire-weather  observation  drive  the 
various  models  that  make  up  the  National  Fire-Danger 
Rating  System.  Observation  time  (early  to  midafternoon) 
elements  are: 

1.  Temperature,  °F, 

2.  Relative  humidity,  percent, 

3.  State  of  the  weather, 

4.  Ten-minute  average  20-ft  windspeed  (mi/h),  and 

5.  Fuel  stick  moisture,  percent. 

Elements  for  the  24-hour  period  ending  at  the  observa- 
tion time  are: 

6.  Duration  of  precipitation,  hours, 

7.  Amount  of  precipitation,  inches, 

8.  Maximum  24-hour  temperature,  °F, 

9.  Minimum  24-hour  temperature,  °F, 

10.  Maximum  24-hour  relative  humidity,  percent,  and 

11.  Minimum  24-hour  relative  humidity,  percent. 
NOTE:  AFFIRMS  will  also  accept  metric  environmental 
inputs. 

Dead  fuels  are  stratified  by  moisture  response  timelag 
classes  (1,  10,  100,  and  1,000  hours);  live  fuels  by  type  of 
vegetation  (grass-forbs  or  woody  shrubs).  Fuel  classifica- 
tions are  discussed  in  the  section  titled  "Classification  of 
Fuel  Components,"  under  "Forest  Fuels";  fuel  moisture 
calculations  are  covered  in  sections  titled  "Dead  Fuel 
Moisture  Models"  and  "Live  Fuel  Moisture  Models." 
Precipitation  duration  affects  fuel  moisture  more 
significantly  than  does  precipitation  amount.  The 
24-hour  precipitation  amount  is  recorded  because  it  is  a 
standard  climatological  element. 

The  traditional  fuel  moisture  sticks  are  used  as  an 
analog  of  the  10-hour  timelag  fuel  class.  A  set  of  sticks 
is  an  array  of  three  1/2-  by  18-inch  ponderosa  pine 
dowels.  The  ovendry  weight  of  fuel  sticks  decreases  with 
extended  weathering,  so  a  correction  for  age  is  applied. 

Fuel  particle  and  fuel  bed  property  values  are  quan- 
tified in  the  system's  20  stylized  fuel  models.  Fuel 
models  consist  of  fuels  information  required  for  input  to 
the  system's  SC  and  ERC  models.  A  fuel  model  contains 
fuel  characteristic  values  typical  of  fuel  descriptions  for 
a  general  cover  type. 


Wind,  slope,  fuel  moistures,  and  fuel  descriptors  (via 
fuel  models)  are  required  to  compute  the  SC  and  ERC 
using  a  modification  of  the  Rothermel  (1972)  fire  spread 
model  (see  section  titled  "Fire  Behavior  Model").  The 
burning  index  computations  are  based  on  Byram's  flame 
length  model  (Byram  1959).  The  specifics  of  the  1978 
NFDRS  indexes  and  components  SC,  ERC,  BI,  IC, 
MCOI,  LOI,  and  FLI  are  discussed  in  sections  titled 
"Fire  Behavior  Components,"  "Ignition  Component, 
Risk,  and  Occurrence  Indexes,"  and  "The  Fire  Load 
Index." 

FOREST  FUELS 

The  delineation  of  wildland  fuel  characteristics  for  fire 
behavior  modehng  and  resultant  fire-danger  rating  is 
based  on  quantitative  descriptions  of  fuel  particle  and 
fuel  bed  properties.  The  important  fuel  particle  proper- 
ties are  size,  density,  chemical  composition,  and  shape  (a 
cylindrical  shape  is  assumed).  Important  fuel  bed 
characteristics  are  fuel  amounts  (dry  weight  load  by  size 
class,  lb/ft2)  and  fuel  bed  depth.  (The  fuel  bed  is  assumed 
to  be  uniform  and  continuous,  and  the  various  fuel 
classes  uniformly  distributed  throughout  the  fuel  bed.) 
Stylized  fuel  models  are  discussed  in  the  section  titled 
"Fuel  Models." 

Classification  of  Fuel  Components 

DEAD  FUELS 

Dead  fuels  are  fuels  in  which  the  moisture  content  is 
exclusively  controlled  by  environmental  conditions- 
temperature,  radiation,  relative  humidity,  and 
precipitation. 

The  relationships  among  environmental  conditions  and 
dead  fuel  moisture  in  the  1978  NFDRS  draw  heavily 
upon  definitions  and  theory  proposed  by  Byram  (1963) 
and  expanded  on  by  Fosberg  (1970).  Byram 
demonstrated  that  the  moisture  content  of  dead  fuels 
drying  under  constant  conditions  follows  an  exponential 
decay  curve.  He  defined  the  timelag  interval,  r,  as  the 
time  required  for  fuels  to  lose  approximately  two-thirds 
of  their  initial  moisture  content  (the  actual  amount  is 
1  — 1/e,  where  e  is  the  base  of  natural  logarithms). 

Byram  defined  the  relative  moisture  content,  ft.,  as: 
ix  =  ix,-  exp(-rT)  (%)  (1) 

where 

T  is  the  period  of  environmental  stress  (h), 
T  is  the  fuel  particle  timelag  (h),  and 
IJ.Q  is  the  initial  moisture  content  of  the  fuel  particle  (%). 
The  relative  moisture  content  (%)  may  also  be  defined: 

fji  =  (EMC  -  ^)/(mo  -  m)  (2) 
where 

'ji  is  the  mean  moisture  content  (%)  over  the  time  period, 
EMC  is  the  equilibrium  moisture  content  (%),  and 
fiQ  is  the  initial  moisture  content  (%)  at  the  onset  of  the 
stress  period. 

Equilibrium  moisture  content  (EMC)  is  the  moisture  con- 
tent dead  fuels  would  obtain  if  left  in  a  steady-state  en- 
vironment long  enough  to  obtain  equilibrium  (no  net 
moisture  exchange).  It  is  computed  using  the  fuel- 
atmosphere  interface  dry  bulb  temperature  and  relative 
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humidity.  Computation  of  EMC  values  is  discussed  further 
in  the  section  titled  "The  1-Hour  Timelag  Fuel  Moisture 
Models,"  and  appendix  C. 

Fosberg  (1970)  then  introduced  similarity  theory  ex- 
pressed by  a  Fourier  number  (Fq).  Using  both  ex- 
perimental and  theoretical  methods,  Fosberg  and  others 
(1970)  showed  the  Fourier  number  to  be  a  universal  con- 
stant for  fuel  moisture.  For  cylindrical  fuels,  the  Fourier 
number  is 

Fo  =  (Tj^)/r2  (3) 
and  for  the  general  case 

Fo  =  (r.a2)/(4xs2)  (4) 
where 

T  is  the  particle's  moisture  timelag  (h), 
u  is  the  characteristic  moisture  diffusivity  of  the  fuel 
particle  (cm2/min), 
r  is  the  fuel  particle  radius  (cm), 

a  is  the  particle's  surface  area-to-volume  ratio  (see  sec- 
tion titled  "Fuel  Particle  Properties"),  and 

Xs  is  a  fuel  particle  shape  factor  (not  used  in  the 
NFDRS). 

In  both  cases,  Fq  =  0.18. 

Experimentally  determined,  diffusivity  inherently  con- 
tains information  on  the  species,  and  thus  automatically 
accounts  for  the  presence  of  waxes  and  resins  that 
hinder  fuel  moisture  transport.  Some  typical  values  of  v 
are  given  in  table  1  for  selected  species. 

Fosberg  (1971b)  then  used  the  Fourier  number  to  solve 
for  a  characteristic  timelag  similarity  coefficient  (f),  and 
defined  the  change  of  moisture  content  for  transient  sur- 
face conditions  as: 

5m/A/x  =  1  -  r-  e<-T/r)  (5) 
where 

8fi  is  the  actual  moisture  content  change  (%)  in  the 
stress  period  (T,  h), 

Afx  is  the  potential  moisture  content  change  (%)  in  the 
stress  period  (T), 

f  is  the  similarity  coefficient  (dimensionless), 

e  is  the  base  of  the  natural  logrithms,  and 

r  is  the  fuel  particle  timelag  (h). 
Afi  may  be  estimated  from  the  difference  in  the  EMC  at 
the  beginning  and  end  of  the  stress  period  as  calculated 
from  the  initial  and  ending  fuel-atmosphere  temperatures 
and  relative  humidities.  This  equation  provides  the  basis 
for  theoretical  calculation  of  dead  fuel  moisture. 

The  range  of  moisture  contents  over  different  time  in- 
tervals is  characteristic  of  fuel  particle  size.  Figure  3  il- 
lustrates a  normal  diurnal  cycle  in  which  the  moisture 
content  of  fine  fuels  vary  widely.  Intermediate-sized 
fuels  vary  within  a  narrower  range,  and  large  fuels 
within  a  very  limited  range  (Gisborne  1928;  Brackebusch 
1975). 

Using  equation  5,  Fosberg  (1971a)  simulated  the 
magnitude  of  the  range  of  moisture  responses  of  fuels 
with  different  timelags  to  weather  cycles  of  1-day  (diur- 
nal), 4-day  (synoptic),  30-day  (planetary),  and  1-year  (an- 
nual). He  used  climatological  data  from  10  locations  to 
represent  the  diurnal  variation  of  boundary  value 
moisture  content  and  data  from  2  locations  to  represent 
the  annual  variations  of  boundary  moisture  content. 


Table  1  .—Typical  diffusivity  values  for  selected  forest  fuels  (from  un- 
published 1972  NFDRS  documentation) 


Fuel  type 


Ponderosa  pine  dowel  (Pinus  ponderosa) 


0.8 
0.6 

0.4 
0.2 

0 

-0.2 
-0.4 
-0.6 
-0.8 
-1.0 


Diffusivity 


cm^/s 


0.159  cm  radius 

8.32  X  10-7 

0.318  cm  radius 

8.32  X  10-7 

0.635  cm  radius 

8.32  X  10-7 

1.0     cm  radius 

8.32  X  10-7 

1.27  cm  radius 

8.32  X  10-7 

2.54  cm  radius 

8.32  X  10-7 

5.08  cm  radius 

8.32  X  10-7 

Douglas-fir  dowel  {Pseudotsuga  menziesii) 

1.0  cm  radius  8.30X10-7 

Lodgepole  pine  dov\/el  (Pinus  contorta) 

1.0  cm  radius  8.42X10-7 

Alaska  cedar  dowel  (Chamaecyparis  nootkatensis) 

1.0  cm  radius  5.94  X  10-7 

Eastern  white  pine  dowel  (Pinus  strobus) 

1.0  cm  radius  10.3  X  10-7 

Ponderosa  pine  needles  (Pinus  ponderosa)  2.0  X  10-8 

Western  white  pine  needles  (Pinus  monticola)  1.18  X  10-8 

Eucalyptus  leaves  (Eucalyptus  obliqua)  4.77  X  10-^ 

Quaking  aspen  leaves  (Populus  tremuloides)  — 

Cheatgrass  (Bromus  tectorum) 

Leaves  — 
Stalks  — 

Plant  2.58  X  10-7 

Fescue  (Festuca  idatioensis) 

Leaves  — 
Stalks  — 
Plants  — 


1.0,- 


 0.6 

 0.4 

 0.2 

 0.1 

 /l-EMC 


10     12    14     16     18    20     22  24 


LOCAL  STANDARD  TIME (H) 


Figure  3.  — Diurnal  variations  of  relative 
moisture  content  (m)  for  fuels  of  increasing 
radii  (adapted  from  Fosberg  1971a,  p.  68). 
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Figure  4  shows  the  relative  ranges  of  moisture  content 
(A/x)  plotted  over  timelag  (t)  for  the  four  timelag  periods. 
The  range  of  moisture  content  in  each  case  period  is 

=  (Mmax  -  /^min)  (6) 

where 

Mmax  is  the  period  maximum  moisture  content  (%),  and 
Mmin  is  the  period  minimum  moisture  content  (%). 


CLASS  (H) 


TIMELAG  (H) 

Figure  4.— Moisture  recovery  curves  for  fuel 
particles  of  increasing  radii  (adapted  from 
unpublished  1972  NFDRS  documentation). 

The  A/i's  for  a  given  time  period  characterize  fuel  par- 
ticles of  a  given  timelag.  As  r  increases,  Afx  decreases  if 
the  period  is  constant.  As  time  periods  increase,  fuel  par- 
ticles with  a  constant  t  exhibit  a  larger  Ajx. 

This  relationship  provided  a  basis  for  dead  fuels 
classification  (Lancaster  1970).  The  recovery  curves,  all 
having  a  reverse  S-shape,  are  characterized  by  three 
straight-Hne  segments.  For  the  diurnal  cycle,  fuels  with 
a  timelag  less  than  2  hours  show  a  high  recovery;  those 
with  timelags  greater  than  20  hours  show  a  low 
recovery.  Fuels  with  timelags  less  than  200  hours  show 
an  intermediate  recovery  range  for  the  30-day  cycle  and 
a  high  recovery  range  for  the  annual  cycle. 

Based  upon  these  differences,  four  groupings  of 
timelags  were  delineated:  fuels  with  timelags  up  to  2 
hours;  from  2  to  20  hours;  from  20  to  200  hours;  and 
more  than  200  hours.  Each  group  responds  uniquely  to 
successively  longer  cycles.  The  approximate  midpoints  of 
these  timelag  ranges  are  used  as  the  basis  for  a 
classification  scheme:  1-hour,  10-hour,  100-hour,  and 
1,000-hour  timelag  fuels. 

Exposed  mosses,  lichens,  and  cured  grasses  and  herbs 
have  timelags  of  1  hour  or  less.  The  uppermost  layer  of 
weathered  conifer  needles  on  a  forest  floor  typically 
respond  to  environmental  stress  within  2  hours,  as  do 
dead  twigs  of  woody  plants  up  to  y4-inch  diameter. 


Fresh  conifer  needles  exhibit  timelags  ranging  from  1  to 
10  hours,  but  approach  1  hour  as  waxes  and  resins  are 
leached  away  (Van  Wagner  1969;  Anderson  and  others 
1978). 

The  10-hour  class  includes  dead  twigs  and  branches 
from  one-fourth  inch  to  1  inch  in  diameter.  Dead  branch- 
wood  in  the  1-  to  3-inch  diameter  class  falls  into  the 
100-hour  timelag  group,  and  dead  logs  and  branchwood 
from  3  to  8  inches  constitute  the  1,000-hour  timelag 
class.  The  dead  fuel  moisture  models  for  each  size  class 
of  dead  fuel  are  presented  in  the  section  titled  "Dead 
Fuel  Moisture  Models,"  with  subsections  on  the  1-hour, 
10-hour,  100-hour,  and  1,000-hour  timelag  fuel  moisture 
models. 

LIVE  FUELS 

In  living  fuels  moisture  content  is  controlled  by  the 
physiological  processes  of  the  plant.  Changes  in  moisture 
of  the  two  classes  of  live  fuels  considered  by  the  1978 
NFDRS  are  both  seasonal  and  short  term.  These  classes 
are  herbaceous  plants  (grasses  and  forbs)  and  woody 
shrubs.  Dynamic  live  fuel  moisture  models  have  been 
developed  to  simulate  the  greening  and  curing  process  of 
these  fuels  through  a  growing  season  and  also  short- 
term  moisture  content  fluctuations  due  to  extreme  en- 
vironmental conditions.  The  live  fuel  moisture  models 
are  discussed  in  the  section  titled  "Live  Fuel  Moisture 
Model." 

Herbaceous  Plants 

Plants  that  do  not  develop  ptersistent  woody  tissues 
such  as  grasses,  forbs,  and  ferns  make  up  the  NFDRS 
herbaceous  fuel  class.  The  herbaceous  class  is  further 
subdivided  into  annual  and  perennial  types.  When  the 
fuel  moisture  falls  below  30  percent,  these  plants  are 
considered  cured  and  the  moisture  content  defaults  to 
that  of  the  1-hour  timelag  fuels.  The  herbaceous  fuel 
moisture  model  is  presented  in  the  section  titled  "Her- 
baceous Fuel  Moisture." 

Woody  Shrubs 

The  second  category  of  live  fuel  in  the  1978  NFDRS  is 
the  perennial  woody  shrub.  These  fuels  are  considered 
dormant  when  the  moisture  content  falls  to  50  percent. 
Above  this  value  these  plants  are  allowed  a  maximum 
moisture  content  of  250  percent  during  the  growing 
season.  The  woody  shrub  fuel  moisture  model  does  not 
allow  the  estimated  shrub  moisture  content  to  fall  so  low 
that  the  shrubs  would  have  to  be  considered  dead.  But 
at  the  low  end  of  their  moisture  range  they  are  con- 
sidered dormant.  The  woody  fuel  moisture  model  is 
presented  in  the  section  titled  "Woody  Fuel  Moisture 
Model." 

Fuel  Particle  Properties 

The  physical  and  chemical  fuel  particle  properties 
described  in  this  section  were  originally  defined  by 
Rothermel  (1972).  These  properties  and  fuel  moisture 
response  to  environmental  stress  (via  fuel  size)  are  fun- 
damental to  rating  fire  danger.  Physical  characteristics 
such  as  fuel  particle  size  and  shape  affect  the  ease  of  ig- 
nition and  rate  of  moisture  exchange  with  the  environ- 
ment. The  energy  potentially  available  through  combus- 
tion is  a  chemical  property. 
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PHYSICAL  PROPERTIES 

Fuel  Particle  Surface  Area-to-Volume  Ratio  (a) 

The  gain  and  loss  of  heat  and  moisture  and  the  evolu- 
tion of  combustible  gases  occur  through  the  surface  of  a 
fuel  particle  and  occur  at  rates  directly  related  to  the 
amount  of  surface  area  per  unit  volume  of  the  particle. 
Thus  the  surface  area-to-volume  ratio  (o)  of  a  fuel  par- 
ticle is  a  very  important  fuel  particle  property.  Fuel  par- 
ticles with  large  o's  will  ignite  more  readily  than  those 
with  relatively  small  o's.  For  a  cylinder,  the  surface  area- 
to-volume  ratio  is  inversely  related  to  the  radius. 

The  surface  area  (neglecting  end  area)  of  a  cylinder,  A, 

is 

A  =  2  TrrL    (ft2)  (7) 


and  its  volume 

V  =  7rr2L  (ft3)  (8) 
where  r  is  the  radius  and  L  is  the  length.  The  surface 
area-to-volume  ratio  for  cylindrical  fuels  is 

a  =  A/V  =  (27rrL)/(7rLr2)  =  2/r  (1/ft).  (9) 

Fuel  Particle  Density  (gp) 

The  next  important  fuel  particle  property  is  the  densi- 
ty, Qp.  This  is  needed  to  determine  the  packing  ratio,  (3 
(see  section  titled  "Fuel  Compactness"),  of  a  fuel  bed.  A 
Qp  value  of  32  lb/ft3  is  used  for  both  live  and  dead  fuels 
in  all  of  the  fuel  models. 

Table  2  contains  typical  values  of  o  and  for  selected 
fuel  particles. 


Table  2.— Typical  fuel  particle  properties 
documentation) 


of  selected  forest  fuels  (fronn  unpublished  1972  NFDRS 


Fuel  Type 


Shape 
Factor 


Surface  area      Cross  sec-  Particle 

vol  ratio         tional  area  density  Timelag 

{a,  cm-1)  (A,  cm2)  (q,  g/cm3)  (r,  h) 


Ponderosa  pine  dowel 
Pinus  ponderosa 

0.159  cm  radius  1 

.318  cm  radius  1 

.635  cm  radius  1 

1.0     cm  radius  1 

1.27  cm  radius  1 

2.54  cm  radius  1 

5.08  cm  radius  1 

Douglas-fir  dowel 
Pseudotsuga  menziesii 

1.0  cm  radius  1 


12.6 
6.3 
3.15 
2.0 
1.57 
.787 
.394 


2.0 


0.079 
.317 
1.27 
3.14 
5.07 

20.3 

81.1 


3.14 


0.418 
.418 
.418 
.418 
.418 
.418 
.418 


.508 


1.5 
6.1 
24.0 
59.0 
97.0 
388.0 
1,551.0 


59.0 


Lodgepole  pine  dowel 
Pinus  contorta 
1.0  cm  radius 


1 


2.0 


3.14 


.432 


59.0 


Alaska  cedar  dowel 
Chamaecyparis  nootkatensis 
1.0  cm  radius  1 

Eastern  white  pine  dowel 
Pinus  strobus 

1.0  cm  radius  1 


Ponderosa  pine  needles 
Pinus  ponderosa 


1.3 


Western  white  pine  needles 

Pinus  monticola  1.3 

Eucalyptus  leaves 

Eucalyptus  obliqua  1 .43 

Quaking  aspen  leaves 

Populus  tremuloides  1.07 

Cheatgrass 

Bromus  tectorum 

Leaves  1.9 

Stalks  — 

Plant  — 

Fescue 

Festuca  idahoensis 

Leaves  2.2 

Stalks  — 

Plant  — 


2.0 


2.0 


57.6 


90.5 


61.8 


139.8 


189.0 
75.8 
145.0 


190.0 
68.6 
182.0 


3.14 


3.14 


6.5X10-4 


2.6X10-4 


6.8X10-4 


7.5X10-5 


1.3X10-4 


.458 


.362 


82.0 


47.8 


3.91 


2.7 


1.57 


.07 


1.7X10-4 
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CHEMICAL  PROPERTIES 

Fuel  particles  have  three  chemical  properties  that  af- 
fect combustion  and  hence  must  be  considered  by  fire- 
danger  rating:  heat  content,  total  mineral  content,  and 
effective  mineral  content. 
Heat  Content  (H) 

The  energy  available  per  unit  mass  of  fuel  through 
combustion  is  called  the  heat  content.  In  woody  fuels,  on 
a  dry  weight  basis,  it  is  relatively  consistent  from  one 
plant  species  to  another.  Heat  content  may  also  be  re- 
ferred to  as  the  heat  of  combustion.  In  the  NFDRS,  this 
value  ranges  from  8,000  Btu/lb  to  9,500  Btu/lb  and 
varies  by  fuel  model.  Table  3  summarizes  some  typical 
values  for  western  fuels. 
Total  Mineral  Content  (Sj) 

Total  mineral  content  (St)  is  the  fraction  of  a  fuel  mass 
composed  of  inorganic  minerals.  The  inorganic  mineral 


content  reduces  the  combustible  fuel  mass  because  only 
the  organic  portion  of  a  fuel  supports  combustion.  The 
tot£d  mineral  content  may  also  be  referred  to  as  total 
ash  content.  A  value  of  5.55  percent  is  used  for  both 
dead  and  live  fuels  in  all  NFDRS  fuel  models. 
Effective  Mineral  Content  (Sg) 

The  effective  (active)  mineral  content  (Se)  affects  fire 
behavior  by  interfering  with  the  chemical  processes  of 
combustion.  The  presence  of  certain  mineral  salts  alters 
the  pyrolysis  process  and  promotes  the  formation  of 
char  and  tar  at  the  expense  of  more  flammable  volatiles. 
Philpot  (1968,  1970)  showed  that  the  rate  of  thermal 
degradation  and  the  amount  of  volatiles  produced  are 
reduced  as  the  silica-free  ash  content  increases.  Table  4 
summarizes  the  total  ash  and  silica-free  ash  content  for 
selected  species.  A  value  of  1  percent  is  used  for  both 
dead  and  live  fuels  in  the  NFDRS  fuel  models. 


Table  3.— Average  heat  content  (heat  of  connbustion)  of  selected  forest  fuels  (from  Kelsey 
 and  others  1979) 


Species 


Wood 


Bark 


Twigs 


Foliage 


-Btu/lb- 


Western  redcedar 

9,700 

8,669 

8,708 

9,630 

Grand  fir 

8,300 

9,641 

8,894 

9,497 

Western  larch 

8,620 

9,162 

9,247 

8,703 

Western  white  pine 

8,620 

9,355 

9,464 

9,040 

Engelmann  spruce 

8,100 

9,616 

9,076 

8,806 

Lodgepole  pine 

8,600 

9,605 

9,371 

9,365 

Western  hemlock 

8,500 

9,943 

8,924 

9,491 

Douglas-fir 

9,200 

10,845 

9,113 

9,265 

Ponderosa  pine 

9,100 

9,452 

10,026 

9,527 

Average 

8,721 

9,461 

9,253 

9,315 

Sample  size 

35 

35 

43 

38 

Table  4.  — Mineral  content  of  selected  plant  materials  (adapted  from  Philpot  1970) 

 Mineral  content  

Species  Part  Total  ash  Silica-free-ash 


Percent  dry  weight- 


Cellulose'' 

0.01 

0.01 

Douglas-fir  (Pseudotsuga  menziesii) 

Wood 

.11 

.11 

Birch  (Botula  sp.) 

Wood 

.18 

.18 

Ponderosa  pine  (Pinus  ponderosa) 

Wood 

.22 

.22 

Poplar  (Populus  sp.) 

Wood 

.36 

.36 

Cheatgrass  (Bromus  tectorum) 

Leaves 

5.27 

1.04 

Medusahead  grass  {Elymus  caput-medusae) 

Leaves 

16.02 

1.17 

Ponderosa  pine  (Pinus  ponderosa) 

Needles 

3.87 

1.55 

Chamise  {Adenostoma  fasciculatum) 

Stems 

2.19 

1.75 

White  pine  (Pinus  monticola) 

Needles 

3.34 

2.54 

Chamise  (Adenostoma  fasciculatum) 

Leaves 

3.63 

3.33 

Aspen  (Populus  tremuloides) 

Leaves 

5.24 

5.24 

Guava  (Psidium  guajava) 

Green  leaves 

6.08 

5.58 

Guava  (Psidium  guajava) 

Dead  leaves 

6.24 

5.74 

Guava  (Psidium  guajava) 

Leaves2 

7.79 

7.29 

Saltbush  (Atriplex  canescens) 

Leaves 

12.89 

12.29 

Salt  tree  (Tamarix  apliylla) 

Leaves 

16.59 

14.53 

Saltbush  (Atriplex  polycarpa) 

Leaves 

15.39 

14.83 

Saltbush  (A.  lentiformis  v.  brewerii) 

Leaves 

19.26 

18.63 

Saltbush  (A.  nuttalli  v.  gardnerii) 

Leaves 

26.78 

23.57 

Saltbush  (A.  gardnerii) 

Leaves 

27.07 

25.27 

''Analytical  fiber  pulp.  Carl  Schleicher  and  Schuel  Co. 
^Treated  with  herbicide  (2-4-D  and  2-4-5-T). 
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Fuel  Bed  Properties 

Five  properties  of  fuel  beds  affect  combustion  and  fire 
behavior: 

1.  Amount  or  loading  of  fuels, 

2.  Compactness  of  fuels, 

3.  Continuity  of  fuels, 

4.  Uniformity  of  fuels,  and 

5.  Fuel  element  arrangement  and  distribution. 

Only  amount  and  compactness,  however,  are  allowed  to 
vary.  In  this  and  other  applications  of  the  Rothermel 
'1972)  fire  spread  model,  the  fuel  bed  is  assumed  to  be 
continuous  and  uniform,  and  the  elements  of  the  dif- 
ferent fuel  classes  are  assumed  to  be  thoroughly 
"mixed"  and  uniformly  distributed  throughout  the 
length,  width,  and  depth  of  the  bed. 

FUEL  LOADS  BY  SIZE  CLASS 

Fuel  loading  expresses  the  mass  of  fuel  per  unit  area 
in  a  fuel  bed.  Combined  with  the  heat  content,  H 
(Btu/lb),  the  net  fuel  load  is  a  factor  in  the  amount  of 
heat  available  per  unit  area  of  the  fuel  bed  (Btu/lb  X 
lb/ft2  =  Btu/ft2).  The  net  fuel  load  (W^)  is  the  total  fuel 
load  (Wo)  less  the  mass  represented  by  the  total  mineral 
content  (St)  of  the  fuel  particle.  In  combination  with  fuel 
bed  depth  (6),  it  determines  the  fuel  bed  bulk  density, 

=  WJ8    (lb/ft3)  (10) 
where 

Wo  is  the  total  dry  fuel  bed  load  (lb/ft2),  and 
8  is  the  fuel  bed  depth  (ft). 

FUEL  COMPACTNESS  (/3,0) 

The  compactness  of  a  fuel  bed  is  expressed  as  fuel  bed 
porosity(0),  or  packing  ratio  ((3).  Packing  ratio  is  the 
ratio  of  the  density  of  the  fuel  bed  (gb)  to  the  fuel  par- 
ticle density  (gp): 

0  =  QhlQp  (dimensionless)  (11) 
and  the  porosity  is  the  dimensionless  ratio  of  the  fuel 
bed  void  volume  to  the  total  fuel  bed  volume: 

0  =  Vv/Vb  =  (Vb  -  Vf)/Vb  (dimensionless)  (12) 
where 

Vy  is  the  bed  void  volume  (empty  space,  fts), 
Vf  is  the  bed  fuel  volume  (fts),  and 
Vb  is  the  total  bed  volume  (ft^). 
0,  (3,  and  Qb       interrelated  as  follows: 

^  =  (Qp  -  Qh)lQp'  (13) 
6  =  1-/3.  (14) 

MOISTURE  OF  EXTINCTION  (m^) 

The  fuel  moisture  content  at  which  a  fire  will  not  prop- 
agate in  a  fuel  bed  is  the  moisture  of  extinction.  In  the 
1972  NFDRS,  30  percent,  which  is  close  to  the  fiber 
saturation  value  for  woody  material,  was  used  as  the 
moisture  of  extinction.  But  Brown  (1972)  and  Blackmarr 
(1972)  found  that  moisture  of  extinction  values  vary 
with  fuel  bed  compactness,  fuel  particle  size,  windspeed, 
and  slope.  They  found  extinction  moisture  contents  rang- 
ing from  12  percent  in  grass  to  40  percent  in  conifer  Ut- 
ter. Current  research  at  the  Northern  Forest  Fire 


Laboratory  should  provide  better  guideHnes  for  future 
use. 

Lacking  better  information,  the  NFDRS  makes  several 
assumptions;  the  first  is  that  live  and  dead  fuels  have 
different  moistures  of  extinction,  and  that  a  specific  fuel 
complex  has  a  constant  moisture  of  extinction  for  its 
dead  fuel  components. 
Dead  Fuel  Moisture  of  Extinction  (nixd) 

The  dead  fuel  moisture  of  extinction  is  a  fuel-model- 
dependent  value,  with  values  ranging  from  15  percent  in 
the  grass  and  brush  models  to  30  percent  for  the 
southern  rough,  pine,  and  pocosin  models. 
Live  Fuel  Moisture  of  Extinction  (mxi) 

A  more  extensive  method  is  used  for  live  fuels  because 
their  moisture  contents  are  typically  much  greater  than 
fiber  saturation.  Fosberg  and  Schroeder  (1971)  assumed 
that  the  heat  produced  by  a  burning  mass  of  fuel  and 
transferred  to  an  equal  mass  of  unburned  fuel  having  a 
30  percent  moisture  content  is  just  sufficient  to  raise  the 
temperature  of  the  unburned  fuel  to  ignition.  As  the 
moisture  content  of  the  unburned  fuel  nears  zero  per- 
cent, only  a  very  small  amount  of  heat  is  needed  to  raise 
the  fuel  temperature  to  ignition.  When  a  linear  relation- 
ship is  assumed,  these  two  points  define  the  relationship 
between  moisture  content  and  effective  heat  energy  (line 
A  in  fig.  5). 


Figure  5.  — Relationship  between  live 
moisture  content  and  effective  heat  energy 
(adapted  from  Fosberg  and  Schroeder  1971, 
p.  3). 

If  the  burning  fuel  is  considered  dead,  and  the  un- 
burned fuel  live,  the  extinction  moisture  content  for 
various  proportions  of  live  and  dead  fuels  can  be 
estimated  over  the  range  of  dead  fuel  moisture  content 
from  zero  to  30  percent. 

Fosberg  and  Schroeder  (1971)  used  the  heat  of  pre- 
ignition  (Qjg)  per  unit  mass  (a  function  of  moisture  con- 
tent) from  the  Rothermel  (1972)  rate  of  spread  model  to 
estimate  the  effective  heat  energy. 
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In  metric  units: 

Qig  =  140  +  620mf  (15) 
where  mf  is  the  1-hour  fuel  moisture  (fraction). 

Effective  heat  transfer  (Qige)  is  a  function  of  fractional 
living  fuel  and  dead  fuel  moisture: 

10 


Qige  =  1,800((1  -  a),a)(l 


mf) 


16) 


where 

a  is  the  fraction  of  living  fuel,  and 
mf  is  the  1-hour  dead  fuel  moisture  (fraction). 
At  the  moisture  of  extinction  (assumed  m^  =  0.3), 

Qig  =  Qige: 

140  +  620mxi  =  1,800((1  -  a),a)(l  -  ^mf)  (17 

where  m^i  is  the  extinction  moisture  content  (fractional) 
of  the  unburned  fuel.  Solving  for  mxi: 

10, 


m^i  =  2.9((1  -  a)  a)a  -  ^mf) 


0.226  (18) 

gives  the  m^i  as  a  fraction.  It  is  multiplied  by  100  to 
convert  to  percent. 

m^l  =  290((1  -  a).  a)(l  -  ^mf)  -  22.6.  (19) 


This  expression  was  modified  by  Albini  (1976b)  using 
the  effective  heating  number  concept  introduced  by 
Frandsen  (1973).  The  heating  number,  a  fuel  load 
weighting  factor  of  the  form  e(— k/a)  is  used  to  compute 
a  dead-to-live  loading  ratio  (W),  and  a  weighted  "fine" 
dead  fuel  moisture,  mf^^  Albini's  final  expression  is 


mxi  =  290W(1  -  mf^v/mxd)  -  22.6 


(20) 


where 

mxi  is  limited  to  the  minimum  mxd  value, 
mxd  is  the  dead  fuel  moisture  of  extinction  for  the  fuel 
bed. 


W  = 


£  Wnij  exp(-138/a  ) 
J  =  i  

2 

E  Wn2j  exp(-500/y  ) 


(21) 


and 


2^mijWij  exp(-138;  a  ^.) 


2 

L  W 


(22) 


2^  W      exp(-138,  a,.) 

j  =  l     nl]       ^  li 


and  where 

Wnij  is  the  net  dead  fuel  load  (j  indicates  1-hour, 
10-hour,  or  100-hour  timelag  class), 

Wn2j  is  the  net  live  fuel  load  (j  indicates  woody  or 
herbaceous), 

mij  is  the  dead  fuel  moisture  content  (fraction), 

CTjj  is  the  surface  area-to-volume  ratio. 

exp(  — 138/CTij)  is  the  dead  fuel  effective  heating 
number,  and 

exp(  — 500/a2j)  is  the  live  fuel  effective  heating  number. 


This  expression,  in  conjunction  with  the  live  fuel 
moisture  model,  is  used  for  computing  moisture  damping 
coefficients  for  herbaceous  and  woody  fuels  in  the  fire 
model. 

FUEL  MODELS 

The  previous  section  described  the  classification  of  fuel 
particles  and  defined  eight  fuel  particle  and  bed 
characteristics  that  are  important  in  fire  behavior  model- 
ing and  resultant  fire-danger  rating.  These  eight 
parameters: 

H  heat  content,  Btu  lb 

g p.... fuel  particle  density,  lb  ft3 

St  total  mineral  content,  % 

Sg.... effective  mineral  content.  % 
Wq... .total  dry  fuel  bed  load,  Ib/ft^ 
5  fuel  bed  depth,  ft 

(7  surface  area-to-volume  ratio,  ft2/ft3,  and 

mx... moisture  of  extinction,  ^'c 
are  the  parameters  that  are  quantified  in  the  XFDRS 
fuel  models.  A  fuel  model  is  a  set  of  the  fuel  particle  and 
fuel  bed  descriptors  required  as  inputs  to  the  Rothermel 
spread  model,  plus  wind  reduction  and  rate-of-spread 
normalizing  factors  that  are  discussed  later. 

In  the  1972  NFDRS.  nine  fuel  models  were  provided. 
The  users  and  the  NFDRS  technical  committee  wanted 
the  number  of  fuel  models  increased  in  the  1978  XFDRS 
to  supposedly  improve  resolution  and  representativeness. 
The  researchers  assigned  to  the  project  did  not  agree 
and  recommended  that  the  number  of  models  be  reduced 
or,  at  the  very  least,  be  held  to  nine.  Their  position  was 
based  on  the  following: 

I.  Fire-danger  ratings  are  required  for  areas  on  the 
order  of  10^  to  10^  acres;  low  on  any  scale  of  resolution. 

2.  It  was  not  logical  to  incorporate  high  resolution 
fuel  descriptors  when  fuel  moistures  and  wind  are 
resolvable  at  such  low  levels  (one  observation  per  day 
per  105  acres). 

3.  The  information  required  to  construct  fuel  models 
was  (and  still  is)  very  limited. 

4.  The  fire  behavior  prediction  system  developed  by 
Rothermel  and  Albini  (FBO  package  reference)  and 
Rothermel  (1983)  was  designed  to  satisfy  the  "1-  to 
100-acre/lO  minutes  to  1-hour  resolution"'  fire  behavior 
prediction  needs  of  the  fire  manager. 

The  technical  committee  prevailed  and  the  result  was 
an  array  of  20  fuel  models. 

Eight  of  the  nine  1972  NFDRS  fuel  models  were 
retained— model  F  was  redefined  to  represent  in- 
termediate age  (7  to  15  years  old)  mixed  chaparral.  (The 
original  model  F  represented  very  young  stands  of  brush 
that  were  essentially  fireproof.)  Data  for  the  original 
nine  fuel  models  were  taken  principally  from  Byram 
(1959);  Philpot  (1968,  1970);  Countryman  and  Philpot 
(1970);  and  Brown  (1970a,b,c).  Twelve  totally  new  models 
were  developed. 

The  selection  of  the  twenty  fuel  situations  to  be 
modeled  was  done  by  the  researchers  subjectively.  The 
first  step  was  a  survey  of  fuel  types  entered  on  fire 
reports  by  nine  Forest  Service  Regions,  the  BLM,  and 
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several  State  fire  management  agencies.  A  consolidated 
list  was  produced  after  eliminating  overlaps  and  duplica- 
tions. From  that  list,  20  fuel  types  were  selected  for 
modeling  based  on: 

1.  Significance  as  a  fire  type, 

2.  Interest  expressed  by  fire  managers,  and 

3.  The  availability  of  fuels  and  fire  behavior  data  with 
which  to  develop  and  test  the  models. 

Some  of  the  fuel  situations  matched  the  assumptions 
of  uniformity,  continuity,  and  uniform  distribution  of 
size  classes  very  well,  and  enough  data  existed  for 
straightforward  modehng.  These  fuels  include  the  grass 
models  (A,  L,  and  N),  slash  models  (I,  J,  and  K),  and  the 
needle  litter  of  southern  plantations,  western  pines,  and 
short  needled  conifers  (P,  U,  G,  and  H).  There  was  heavy 
reliance  on  Hough  and  Albini's  work  (1978)  with  the 
palmetto-gallberry  association  to  put  the  southern  rough 
(D)  and  pocosin  (O)  models  together  (D  emulates  a  7-year 
rough).  The  B  and  F  models  (California  mixed  chaparral 
and  intermediate  brush)  approximate  the  models 
developed  by  Rothermel  and  Philpot  (1973)  and  adapted 
by  Van  Gelder  (1976)  for  FIRECAST. 

The  pine  grass  savanna  (C)  and  sagebrush-grass  (T)  are 
combinations  of  grass  and  litter,  and  grass  and  shrubs; 
they  violate  the  assumptions  of  uniformity  and  reflect 
art  more  than  science. 

The  models  constructed  for  Alaska— black  spruce  and 
tundra  (Q  and  S)— are  terribly  artificial,  but  the  outputs 
are  reasonable,  judging  from  the  limited  amount  of  fire 
behavior  data  available  from  Alaskan  and  Canadian 
sources. 

The  fuel  modeling  procedure  incorporated  a  great  deal 
of  repetitive  building  and  evaluating.  As  much  informa- 
tion as  possible  about  a  particular  fuel  situation  was 
gathered,  but  "best  guesses"  were  used  for  a  fuel  model 
when  nothing  concrete  was  available.  The  prototype  fuel 
model  was  then  tested  using  archived  fire-weather  data 
to  generate  daily  values  of  rate-of-spread  and  flame 
length.  The  predictions  were  compared  to  fire  behavior 
observations,  or  if  none  were  available,  the  results  were 
subjectively  evaluated  by  the  researchers. 

Once  a  preliminary  set  of  fuel  models  had  been 
developed,  the  fire-danger  estimates  produced  using 
these  models  were  evaluated  by  fire  managers.  The 
evaluation  included: 

1.  The  rate-of-spread  (SC)  and  flame  length  (0.1  X  BI) 
predictions,  and 

2.  The  seasonal  profile  of  the  NFDRS  indexes,  com- 
ponents, and  live  fuel  moistures. 

The  fire-weather  data  sets  for  the  field  tests  were 
specified  by  the  evaluators.  They  selected  a  wet  year,  a 
dry  year,  and  any  other  period  of  time  that  offered 
characteristics  a  fire-danger  rating  system  should  pick  up. 

Parameters  were  adjusted  and  the  performance 
reevaluated  several  times  in  some  cases  before  model 
parameters  were  settled  on.  Table  5  Hsts  the  20  fuel 
models  in  the  1978  NFDRS.  Fuel  model  parameters  are 
summarized  in  appendix  B.  Narrative  descriptions  and  a 
key  for  selecting  an  appropriate  fuel  model  are  included 
in  the  1978  NFDRS  summary  publication  (Deeming  and 
others  1977). 


Table  5.  — List  of  fuel  models  In  the  1978  NFDRS  (Models  A-E, 
G-l  were  included  in  1972  NFDRS;  model  F 
represents  a  different  fuel  than  in  1972  NFDRS) 


model 

General  description 

A 

Western  annual  grasses 

B 

California  mixed  chaparral 

C 

Pine  grass  savanna 

D 

Southern  rough 

E 

Hardwoods  (winter) 

F 

Intermediate  brush 

G 

Short  needle  pine  (heavy  dead) 

H 

Short  needle  pine  (normal  dead) 

1 

Heavy  logging  slash 

J 

Intermediate  logging  slash 

K 

Light  logging  slash 

L 

Western  perennial  grass 

N 

Sawgrass 

0 

High  pocosin 

P 

Southern  pine  plantation 

Q 

Alaskan  black  spruce 

R 

Hardwoods  (summer) 

8 

Tundra 

T 

Sagebrush-grass 

U 

Western  long-needled  conifer 

DEAD  FUEL  MOISTURE  MODELS 

The  dead  fuel  moisture  models  in  the  1978  NFDRS  are 
based  on  the  theory  described  in  section  titled 
"Classification  of  Fuel  Components,"  except  when  fuel 
moisture  sticks  were  weighed  as  part  of  the  fire-weather 
observation.  In  those  cases,  the  1-hour  and  10-hour  fuel 
moistures  are  computed  using  the  analog  stick  moisture 
content.  The  100-hour  and  1,000-hour  fuel  moistures  are 
computed  from  Fosberg's  theoretical  solutions,  using  an 
average  equilibrium  moisture  content  (called  boundary 
conditions)  for  24  and  168  hours,  respectively.  The  solu- 
tions are  based  on  two  equations,  one  from  Fosberg 
(1970): 

d/ilAli  =  (mc  -  mCo)/(EMC  -  mCo)  (23) 

and  equation  5  from  section  titled  "Classification  of  Fuel 
Components": 

6/x/Am  =  1  -f(exp(-T/r))  (5) 

where 

dfi  is  the  actual  moisture  content  change  (%)  in  the 
stress  period,  T, 

Afi  is  the  potential  moisture  content  change  (%)  for  the 
stress  period  T,  (mc^  —  mc), 

mCo  is  the  moisture  content  (%)  at  the  beginning  of 
the  stress  period  T, 

mc  is  the  moisture  content  (%)  at  the  end  of  the  stress 
period  (T+1), 

EMC  is  the  mean  equilibrium  moisture  content  (%)for 
the  stress  period  (T), 

f  is  the  dimensionless  similarity  coefficient, 

T  is  the  fuel  particle  timelag  (h),  and 

T  is  the  simulation  (stress)  period  timestep  (h). 
Repeating  again,  f  is  empirically  derived  and  valid  only 
for  particular  combinations  of  T  and  t  at  a  specific  time 
of  the  day. 
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The  1-Hour  Timelag  Fuel  Moisture  Model 

The  1978  NFDRS  gives  two  methods  to  calculate  the 
1-hour  fuel  moisture,  but  that  developed  for  the  Califor- 
nia wildland  fire-danger  rating  system  (USD A  1958, 
revised  1968)  is  preferable: 

mci  =  (4EMC  +  mcioi^)/5  (%)  (24) 
where 

mciojj  is  the  age-corrected  fuel  stick  moisture  (%,  see 
next  section  titled  "The  10-Hour  Timelag  Fuel  Moisture 
Model"),  and 

EMC  is  the  equilibrium  moisture  content  (%) 
calculated  using  the  temperature  and  relative  humidity 
at  the  fuel-atmosphere  interface. 

Observation  time  relative  humidity  and  temperature  at 
the  fuel-atmosphere  boundary  layer  are  estimated  from 
lapse  rates  determined  by  Haltiner  (1975)  that  are  depen- 
dent on  sky  cover. 

To  adjust  standard  exposed  instrument  readings  (4.5 
feet  from  ground  in  shelter)  of  relative  humidity  and 
temperature  to  fuel  level,  the  factors  in  table  6  are  used. 

Table  6.  — Fuel-atmosphere  interface  temperature  and  relative 
humidity  adjustment  factors  (from  Haltiner  1975) 


Fraction  sky  coversky  condition 

0.0-0.1:  0.1-0.5:  0.6-0.9:  0.9-1.0: 
Clear       Scattered    Broken  Overcast 


Dry  bulb  temperature, 

°F  (add)  25  19  12  5 

Relative  humidity 

(multiply)  0.75  0.83         0.91  1.0 


If  the  fuel  stick  moisture  (mcioj^  is  not  reported,  the 
method  described  by  Fosberg  and  Deeming  (1971)  is 
used.  Combining  equations  5  and  23  and  solving  for 
mci,  Fosberg's  basic  diffusion  equation  becomes 

mci  =  mco  +  (EMC  -  mCo)(l  -  r(exp(-T/r))  (25) 
where 

mci  is  the  1-hour  timelag  fuel  moisture  (%)  at  time  T 
and 

mCo  is  the  1-hour  timelag  fuel  moisture  at  time  (T  — 1). 
The  equation  was  solved  sequentially  for  forty-eight 
3D-minute  time  periods,  ending  at  a  time  in  the 
midafternoon. 

Using  diurnal  data  from  the  O'Neill,  Nebr.,  experiment 
(Lettau  and  Davidson  1957),  a  T  value  of  0.5  hours,  and 
T  and  f  values  of  1.0  yielded  a  simplified  expression  for 
1-hour  fuel  moisture  for  midafternoon  observation  time: 

mci  =  1.03  EMC  (%)  (26) 
where  EMC  is  the  calculated  equilibrium  moisture  con- 
tent at  the  fuel-atmosphere  interface. 

The  10-Hour  Timelag  Fuel  Moisture  Model 

The  preferred  method  for  obtaining  10-hour  fuel 
moisture  is  from  fuel  moisture  sticks,  corrected  for  stick 


aging.  When  analog  stick  weight  (Wf)  is  measured: 

mcio  =  AiC  +  BC(Wf  -  100)  (%)  (27) 
with 

Ai  =  a/60  (28) 
B  =  1  +  0.02(a/30)  (29) 
C  =  Cc/4  (30) 
where 

Wf  is  the  fuel  stick  weight  (g), 

a  is  the  number  of  days  since  the  sticks  were  set  out, 
Cc  is  the  1978  NFDRS  climate  class  (1  to  4),  and 
100.0  is  the  ovendry  weight  of  the  fuel  sticks  (g). 

The  age  correction  as  developed  by  Haines  and  Frost 
(1978)  was  changed  to  be  functional  on  climate  class 
because  the  original  function  over-estimated  weight  loss 
in  dry  climates.  NFDRS  climate  classes  are  discussed  in 
the  section  titled  "Live  Fuel  Moisture  Model." 

When  fuel  sticks  are  not  measured,  the  10-hour  fuel 
moisture  is  estimated  in  a  manner  similar  to  that  for  the 
1-hour  fuel  moisture  (also  reported  by  Fosberg  and 
Deeming  1971). 

The  diurnal  O'Neill,  Nebr.,  data  were  used  and  the 
basic  diffusion  equation  (eq.  25)  was  solved  in  4-hour 
time  steps  using  f  =  0.87,  T  =  4  hours,  and  t  =  10 
hours  to  produce  this  estimate  for  a  midafternoon 
observation: 

mcio  =  1.28  EMC  (%)  (31) 
where 

EMC  is  the  same  EMC  used  to  calculate  the  1-hour 
fuel  moisture. 

This  model  works  well  for  early  afternoons  in  strong 
continental  areas  at  the  approximate  latitude  of 
Nebraska  in  the  late  summer.  It  tends  to  underpredict 
stick  readings  under  other  conditions. 

FORECASTING  THE  10-HOUR  FUEL  MOISTURE  IN 
AFFIRMS 

Accurately  forecasting  the  next  day's  10-hour  timelag 
fuel  moisture  (TLFM)  requires  a  more  complex  model 
that  utilizes  24-hour  maximum  and  minimum 
temperatures  and  relative  humidities.  The  24  hours  be- 
tween observation  times  is  divided  into  two  periods:  the 
first  from  the  basic  observation  time  of  day  1  (1400 
hours)  to  0600  hours  of  day  2  (assumed  to  be  the  time  of 
minimum  temperature  and  maximum  relative  humidity 
for  the  24-hour  period);  and  the  second  from  0600  of  day 
2  to  the  basic  observation  time  of  day  2. 

The  average  temperatures  and  relative  humidities  for 
the  two  periods  are  estimated  from  predictions  or  obser- 
vations for  the  three  occasions  (times).  The  predicted 
duration  of  precipitation  during  the  two  periods  is  in- 
cluded in  the  calculation  of  average  boundary  conditions 
for  each  period: 

Di  =  [(16  -  Pdi)EMCi  +  (2.7pdi  +  76)pdi]/16  (32) 
D2  =  [(8  -  Pd2)EMC2  +  (2.7pd2  +  76)pd2]/8  (33) 
where 

Di  is  the  average  boundary  condition  for  the  16-hour 
period  from  1400  on  day  1  to  0600  on  day  2, 
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Pdi  and  pd2  are  the  precipitation  durations  for  the  two 
periods  (h), 

EMCi  and  EMC2  are  the  equilibrium  moisture  con- 
tents for  the  two  periods  derived  from  average 
temperature  andrelative  humidity  for  each  period,  cor- 
rected to  the  fuel-atmosphere  interface. 
The  precipitation  duration  effect  function  (2.7pd  +  76)  is 
from  Fosberg  (1972). 

The  next  step  is  to  predict  the  10-hour  fuel  moisture 
at  the  end  of  the  first  period  (0600  hours)  using  the 
observed  (or  computed)  10-hour  fuel  moisture  at  1400  on 
day  1  and  Di  as  initial  values.  The  predicted  value  at 
0600  (day  2)  then  becomes  the  initial  value  for  predicting 
the  10-hour  timelag  fuel  moisture  at  1400  on  day  2: 

mcio,  =  mcio  -  (Di  -  mcio  )(1  -  1.1  exp(-16/10)) 

(34) 

and 

mcio  =  mcio,-  (D2  -  mcioJd  -  0.87  exp(-8/10)) 

(35) 

where 

mcio^  is  the  10-hour  TLFM  at  1400  of  day  1, 

mcio°  is  the  10-hour  TLFM  at  0600  of  day  2, 

mcio  is  the  forecasted  10-hour  TLFM  at  1400  of  day  2, 

0.87  in  eq.  35  is  the  f  derived  for  eq.  31,  and 

1.1  in  eq.  34  is  the  f  derived  for  early  mornings  (0600). 

AFFIRMS  also  allows  direct  predictions  of  10-hour 
fuel  moistures  if  methods  such  as  that  by  Cramer  (1961) 
have  been  developed  for  local  use.  Thus,  if  the  predicted 
10-hour  moisture  content  is  entered,  the  computational 
model  is  skipped. 

The  100-Hour  Timelag  Fuel  Moisture 
Model 

The  100-hour  model  in  the  1978  NFDRS  first  com- 
putes a  24-hour  average  boundary  condition  as  an  input 
to  the  basic  diffusion  equation.  The  boundary  condition 
is  determined  from  precipitation  duration,  maximum  and 
minimum  temperature,  and  relative  humidity,  and  is  for 
the  24-hour  period  from  observation  time  on  day  1  to 
observation  time  on  day  2.  When  some  of  the  above  ob- 
servation elements  are  missing  (principally  in  FIRDAT), 
cruder  methods  are  used  for  100-hour  timelag  fuel 
moisture  computations  (Cohen  and  Deeming,  in 
preparation). 

The  similarity  coefficient  (f)  in  the  1978  model  was 
selected  to  produce  the  same  solution  as  would  the  1972 
model,  when  the  precipitation  duration  is  24  hours.  The 
boundary  condition  is 

D  =  [(24  -  Pd)EMC  +  (0.5pd  +  41)pd]/24  (36) 
and  the  basic  diffusion  equation  becomes 
mcioo  =  mciooo  +  <D  -  mcioo^) 

(1  -  0.87  exp(-24/100))  (37) 

where 

D  is  the  24-hour  average  boundary  condition  (%), 
EMC  is  a  weighted  24-hour  average  equilibrium 
moisture  content  (%),  calculated  from  the  day's  max- 
imum and  minimum  temperatures  and  relative 


humidities  (see  below), 

Pd  is  the  precipitation  duration  for  the  24-hour 
period,  and 

mciooo  and  mcioo  are  the  observed  and  predicted 
100-hour  timelag  fuel  moistures  for  observation  times  at 
days  1  and  2,  respectively. 

The  1978  computerized  100-hour  timelag  fuel  moisture 
model  calculates  the  24-hour  average  equilibrium 
moisture  content  (EMC,  %)  more  effectively  than  in  the 
1972  NFDRS.  The  manual  1978  NFDRS  (and  the  entire 
1972  system)  calculate  EMC  using  arithmetic  averages 
of  a  day's  maximum  and  minimum  temperature,  and 
maximum  and  minimum  relative  humidity^  

The  1978  computerized  NFDRS  derives  EMC  for  both 
the  100-hour  and  1,000-hour  timelag  fuels  from  a 
weighted  average  of  two  EMC's:  one  calculated  from  the 
day's  maximum  temperature  and  minimum  relative 
humidity  (EMCmin)'  and  a  second  calculated  from  the 
day's  maximum  relative  humidity  and  minimum 
temperature  (EMCmax)- 

Weighting  is  based  on  day  length,  which  improves 
seasonal  response  of  NFDRS  fuel  moisture  models.  As 
day  length  shortens,  nighttime  conditions  (EMC  ^ax)  are 
given  more  weight,  thus  promoting  moisture  recovery  in 
the  100-hour  and  1,000-hour  timelag  fuels.  Conversely,  as 
day  length  increases,  daytime  conditions  (EMCmin)  are 
given  more  weight,  thus  promoting  additional  moisture 
loss  in  the  heavier  fuels. 

Under  this  scheme  the  weighted  EMC  is  calculated  by: 

EMC  =  [(n)EMCmin  +  (24  -  n)EMCmax]/24  (38) 

where  n  is  the  hours  of  daylight  (refer  to  appendix  D  for 
computational  details  and  NFDRS  index  response). 

Manual  calculation  of  the  100-hour  timelag  fuel 
moisture  in  the  1978  NFDRS  is  simpler  than  in  the  1972 
method  because  no  arithmetic  is  necessary  (both  systems 
assume  a  constant  24-hour  mean  temperature  to  simplify 
the  calculations).  Using  average  relative  humidity  and 
precipitation  is  a  compromise,  but  it  does  retain  most  of 
the  calculated  moisture  response. 

There  are  three  types  of  missing  data  that  are  con- 
sidered in  the  100-hour  fuel  moisture  computations.  If 
only  precipitation  duration  is  missing,  it  is  estimated  by: 

Pd  =  (Pa  +  0.02)/pr  (39) 
rounded  to  the  nearest  whole  hour,  where 

Pd  is  the  precipitation  duration  (limited  to  a  maximum 
of  8  hours). 

Pa  is  the  24-hour  precipitation  amount  (inches),  and 

Pr  is  0.25  inch/h  for  cHmate  classes  1  and  2,  and  0.05 
inch/h  for  climate  classes  3  and  4. 

If  extreme  values  for  relative  humidity  are  missing, 
but  temperature  extremes  are  available,  maximum  and 
minimum  relative  humidities  are  estimated  by  assuming 
the  air  mass  does  not  change  over  the  24  hours  and  the 
specific  humidity  is  constant.  The  day's  observation 
time  relative  humidity  is  combined  with  the  day's 
temperature  extremes  to  estimate  relative  humidity  ex- 
tremes (Cohen  and  Deeming,  in  preparation). 

If  daily  temperature  extremes  are  missing,  EMCmin  is 
assumed  to  be  the  observation  time  EMC,  and  EMCmax 
is  a  climate  class  dependent  default  value.  For  climate 
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classes  1  and  2,  EMCmax  is  15  percent;  for  classes  3  and 
4,  it  is  23  percent. 

The  weighted  EMC  is  then  computed  from  equation 
38;  the  daily  boundary  condition  from  equation  36;  and 
the  100-hour  moisture  content  from  equation  37. 

The  1,000-Hour  Timelag  Fuel  Moisture 
Model 

The  1,000-hour  timelag  fuel  moisture  model  retains  the 
basic  format  of  the  100-hour  model,  but  requires  an 
assessment  of  the  average  boundary  conditions  for  7 
days  or  168  hours  instead  of  24  hours  (Fosberg  and 
others  1981).  A  7-day  interval  was  chosen  so  that  the 
calculation  would  be  done  the  same  day  each  week  in  the 
manual  version. 

The  mean  7-day  boundary  value  (D,  %): 

D  =  (Di  +  D2  +  ...  +  D7)/7  (40) 
is  computed  from  the  24-hour  average  boundary  values 
(day  length  weighted)  for  the  previous  7  days: 

Dn  =  [(24  -  PdJEMCn  +  (2.7pdn  +  76)pdn]/24  (41) 
where  n  denotes  the  nth  day  and  EMCn  is  weighted  as 
described  in  the  100-hour  model. 
The  diffusion  equation  then  becomes 

mciooo  =  mciooOo+  ^^^^lOOOo  "  D) 

(1  -  0.82  exp(- 168/1,000))  (42) 

where 

mciooo  is  the  predicted  1,000-hour  fuel  moisture, 
mciooOo  initial  value  from  7  days  prior, 

D  is  the  7-day  average  boundary  condition  (%),  and 
0.82  is  the  derived  f  value. 

In  manual  calculations,  the  7-day  change  in  1,000-hour 
timelag  fuel  moisture  is  calculated  with  a  nomogram  and 
manually  added  to  (or  subtracted  from)  the  1,000-hour 
timelag  fuel  moisture  from  7  days  previous.  This 
arithmetic  is  not  buried  in  the  nomogram  because  the 
7-day  change  is  a  variable  needed  for  input  into  the 
NFDRS  live  fuel  moisture  model.  As  a  result,  the 
1,000-hour  timelag  fuel  moisture  must  be  calculated  for 
all  fuel  models  in  spite  of  the  absence  of  1,000-hour  fuels 
in  many.  The  computer  version  does  this  bookkeeping 
automatically. 

RESPONSE  TO  DROUGHT 

By  itself  the  1,000-hour  timelag  fuel  moisture  is  a 
good  indicator  of  the  severity  of  a  medium  term  (4-  to 
6-month)  drought  event.  The  1,000-hour  timelag  fuel 
moistures  calculated  using  data  from  different  climate 
classes  showed  minimum  values  during  recent  drought 
years  to  be  characteristic  of  the  climate  class.  The 
following  tabulation  summarizes  some  cases  studied: 

Season  minimum 
1,000-hour  timelag 
Location  Year    fuel  moisture,  % 


Walport,  Greg. 

1973 

21 

Fort  Meyers,  Fla. 

1970 

14 

Ely,  Minn. 

1976 

11 

Bonners  Ferry,  Idaho 

1973 

9 

Lytle  Creek,  Calif. 

1972 

5 

As  another  example,  the  following  tabulation  portrays 
the  frequency  of  seasonal  minimum  1,000-hour  timelag 
fuel  moisture  values  at  Ninemile,  Mont.,  for  the  years 
1965  —  75.  In  1966  and  1967,  both  very  severe  fire  years, 
the  1,000-hour  TLFM  dropped  to  9  percent. 

Season  minimum 
1,000-hour 

TLFM  (%)  Frequency 

9  2 

10  0 

11  0 

12  3 

13  4 

14  1 

15  1 

LIVE  FUEL  MOISTURE  MODELS 

Improvements  in  the  Rothermel  (1972)  fire  model 
made  it  possible  to  treat  live  fuel  more  realistically  in 
the  NFDRS  update.  Live  fuels  can  now  act  as  either  a 
heat  sink  or  heat  source.  Live  fuels  become  a  heat  source 
when  their  moisture  content  drops  enough  to  allow 
dessication  and  ignition  during  dead  fuel  combustion.  If 
their  moisture  content  remains  above  a  critical  level  (live 
fuel  moisture  of  extinction),  live  fuels  do  not  burn  but 
act  as  a  heat  sink  (Albini  1976a). 

The  live  fuel  moisture  model  provides  more  analytical 
and  consistent  estimates  of  live  fuel  moisture,  replacing 
the  herbaceous  vegetation  transects  required  by  the 
1972  NFDRS.  The  following  section  draws  heavily  upon 
Burgan  (1979). 

Although  it  is  not  based  on  rigorous  principles  of 
plant  physiology,  the  live  fuel  model  broadly  approx- 
imates the  moisture  content  of  living  herbaceous  plants, 
and  the  leaves  and  twigs  of  small  woody  shrubs. 

In  the  original  live  fuel  model  proposed  by  Rothermel, 
plant  moisture  was  treated  as  a  function  of  the  Keetch- 
Byram  (1968)  drought  index.  Several  empirical  factors 
required  to  control  plant  response  to  drying  and  wetting, 
however,  could  not  be  derived  for  all  climates. 

Development  of  the  1978  NFDRS  led  to  the  discovery 
that  the  1,000-hour  timelag  fuel  moisture  function 
responded  to  wetting  and  drying  cycles  similar  to  those 
expected  for  live  fuel.  Thus,  the  1,000-hour  timelag  fuel 
moisture  serves  as  the  basic  meteorological  factor  for 
calculating  live  woody  fuel  moisture  and  is  used  in  a 
modified  form  to  emulate  the  wetting  and  drying  cycles 
of  live  herbaceous  fuels. 

Plants  adapted  to  various  moisture  regimes  respond 
differently  to  rainfall  anomalies:  those  adapted  to  moist 
environments  lose  moisture  faster  during  drought  than 
those  from  dry  environments.  The  1978  live  fuel 
moisture  model,  therefore,  provides  drying  rates  by 
climate  class. 

Although  the  United  States  can  be  divided  into  many 
climatic  zones,  four  classes  were  judged  adequate  to  pro- 
vide the  broadscale  plant  moisture  responses  needed  for 
rating  fire  danger.  These  four  climate  classes  are 
adapted  from  Thornthwaite's  (1931)  earliest  climate 
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classification  system.  Arid  and  semiarid  provinces  are 
grouped  into  a  single  class;  true  desert  is  of  little  prac- 
tical concern  to  fire  management.  And  when  considering 
fire  behavior,  the  subhumid  province  belongs  to  the 
humid  province  rather  than  the  dry,  subhumid  province. 
The  climate  class  descriptions,  general  geographic  areas 
to  which  they  apply,  and  the  vegetative  characteristic  of 
each  are  provided  in  table  7.  Figure  6  maps  the  chmate 
classes  for  the  contiguous  United  States  and  Alaska. 

Climate  class  defines  different  linear  drying  rates  for 
annuals,  perennials,  and  woody  plants,  but  within  a  par- 
ticular chmate  class,  a  single  drying  rate  is  assumed  for 
live  woody  plants  throughout  a  growing  season.  In  live 
herbaceous  plants,  drying  occurs  in  two  stages:  green 
stage  when  herbaceous  moisture  is  above  120  percent, 
and  transition  stage  when  moisture  is  between  120  and 
30  percent.  In  the  green  stage  both  annuals  and  peren- 
nials dry  at  the  same  rate,  but  in  the  transition  stage 
annuals  exhibit  faster  drying  rates  than  perennials. 


The  endpoints  of  the  drying  curves  are  required  to 
define  drying  rates.  First,  weather  data  from  several 
locations  within  each  climate  class  were  used  to  plot 
seasonal  1,000-hour  timelag  fuel  moisture  profiles. 
Typical  plant  moistures  for  each  climate  class  and 
vegetation  type  (woody,  annual,  and  perennial)  were  then 
matched  with  these  profiles.  Because  measurements  of 
seasonal  moisture  variation  in  woody  and  herbaceous 
plants  were  not  generally  available  from  different 
regions  of  the  United  States,  the  live  fuel  model  was 
calibrated  to  produce  reasonable  moisture  values.  Table 
8  shows  the  specific  performance  objectives. 

Maximum  live  fuel  moistures  were  chosen  to  be  250 
percent  for  herbaceous  fuels  and  200  percent  for  live 
woody  shrubs,  at  25  percent  1,000-hour  timelag  fuel 
moisture.  Towards  the  lower  end  of  the  moisture  scale 
herbaceous  plants  were  considered  cured  at  30  percent 
moisture  content  and  woody  plants  dormant  if  their 
moisture  content  dropped  to  50  percent.  The  minimum 


Table  7.— The  1978  NFDRS  climate  classes  (from  Deeming  and  others  1977) 


Thornthwalte 

NFDRS 

(1931) 

climate 

humidity 

Characteristic 

class 

province 

vegetation 

Regions 

1 

Arid 

Desert  (sparse 

Sonoran  deserts  of  western  Texas, 

grass  and  scat- 

New Mexico,  southwest  Arizona, 

tered  shrubs) 

southern  Nevada,  and  western  Utah; 

and  the  Mojave  Desert  of  California. 

Semiarid  Steppe  (short         The  short  grass  prairies  of  the  Great 

grass  and  shrubs)  Plains;  the  sagebrush  steppes  and 

pinyon-juniper  woodlands  of  Wyoming, 
Montana,  Idaho,  Colorado,  Utah, 
Arizona,  Washington,  and  Oregon;  and 
the  grass  steppes  of  the  central  valley 
of  California. 


Subhumid  (rain-     Savanna  The  Alasl<an  interior;  the  chaparral  of 

fall  deficient  in      (grasslands,  Colorado,  Arizona,  New  Mexico,  the 

summer)  dense  bush,  and     Sierra  Nevada  foothills,  and  southern 

open  conifer  California;  oak  woodland  of  California; 

forests)  ponderosa  pine  woodlands  of  the 

West;  and  mountain  valleys  (or  parks) 
of  the  Northern  and  Central  Rockies. 


3 

Subhumid  (rain- 
fall adequate  in 
all  seasons) 

Savanna 
(grasslands  and 
open  hardwood 
forests) 

Blue  stem  prairies  and  blue  stem-oak- 
hickory  savanna  of  Iowa,  Missouri,  and 
Illinois. 

Humid 

Forests 

Almost  the  entire  eastern  United 
States;  and  those  higher  elevations  in 
the  West  that  support  forests. 

4 

Wet 

Rain  forest  (red- 
woods, and 
spruce-cedar- 
hemlock) 

Coast  of  Northern  California,  Oregon, 
Washington,  and  southeast  Alaska. 
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NFDRS  CLIMATE  CLASS 


1  □ 

3  m 

Figure  6.  — 1978  NFDRS  climate  class  (from  Deeming  and  others  1977,  p.  5). 


Table  8.  — Minimum  allowable  moisture  content  of  live  fuels  in  live 
fuel  models  (from  Burgan  1979) 


Shrubs, 

Grasses  and  forbs  twigs,  and 

Type  of  season  Annuals  Perennials  foliage 


 Percent  

Wet  30  >80  >110 

(Late  cure) 

Normal  <30  50-80  80-100 

Dry  <30  <50  50-80 

(Early  cure) 
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woody  and  herbaceous  moistures  were  matched  with 
typical  minimum  1,000-hour  timelag  fuel  moistures  and 
minimum  XI, 000-hour  values,  respectively,  for  each 
climate  class.  Table  9  provides  the  slopes  and  intercepts 
of  the  drying  curves  for  each  climate  class. 

A  dynamic  live  fuel  moisture  algorithm  simulates 
greening  and  curing  of  herbaceous  fuels  by  transferring 
fuel  load  between  the  live  herbaceous  class  and  the 
1-hour  timelag  class  as  the  herbaceous  fuel  moisture 
varies  between  30  and  120  percent  during  the  growing 
season. 

Rothermel  provided  the  empirical  data  used  for 
development  and  initial  testing  of  the  live  fuel  moisture 
model.  He  used  the  xylene  distillation  technique  for 
determining  moisture  content  to  construct  profiles  of 
herbaceous  and  woody  plants  near  Missoula,  Mont.,  dur- 
ing the  1975  and  1976  fire  seasons.  Developed  to 
emulate  these  particular  moisture  profiles,  the  model 
was  then  adjusted  to  produce  reasonable  profiles  for  the 
rest  of  the  United  States. 

Herbaceous  Fuel  Moisture 

Although  relating  herbaceous  fuel  moisture  directly  to 
the  1,000-hour  timelag  fuel  moisture  proved  reasonable 
for  periods  of  drying,  excessive  fuel  moisture  recovery 
was  predicted  during  periods  of  precipitation.  A  function 
was  developed  to  decrease  at  the  same  rate  as  the 
1,000-hour  timelag  fuel  moisture,  but  to  increase  more 
slowly. 

During  the  growing  season  this  value  is 

Xiooo  =  Xyiooo  +  kikslAmciooo)  (43) 
where 

Xiooo  is  the  day's  live  fuel  moisture  recovery  value,  (%), 
Xyiooo  is  yesterday's  Xiooo  value,  (%), 

ki        is  a  drying/ wetting  factor  (dimensionless), 
ki  =  1  if  Amciooo  <  0  (drying  regime) 
ki  =  0.0333(mciooo)  +  0.1675  otherwise 

k2        is  a  temperature  factor  (dimensionless), 
k2  =  1.0, 

if  the  day's  average  temperature  >  50  F  (10°  C) 
k2  =  0.6, 

if  the  day's  average  temperature  <  50  F  (10°  C) 
Amciooo  is  the  24-hour  change  in  mciooo  i^^). 


ki  limits  the  increase  in  herbaceous  fuel  moisture  due 
to  precipitation.  It  is  scaled  to  allow  the  Xiqoo  value  to 
duphcate  the  1,000-hour  timelag  fuel  moisture  when  the 
1,000-hour  timelag  fuel  moisture  is  25  percent  or  more. 

k2  reduces  the  response  of  the  Xjooo  value  to  compen- 
sate for  plants'  slower  physiological  processes  during 
cool  weather.  Figure  7  compares  the  1,000-hour  timelag 
moistures  with  the  Xiqoo  value. 

Prior  to  spring  green-up,  herbaceous  fuel  is  assumed  to 
be  completely  cured,  letting  herbaceous  fuel  moisture 
(mch)  equal  the  1-hour  TLFM. 

During  spring  green-up,  Hve  herbaceous  fuel  moisture 
increases  gradually  from  the  1-hour  timelag  fuel 
moisture.  The  hve  fuel  model  accommodates  gradual 
green-up  at  the  request  of  eastern  United  States  users 
who  feared  that  the  instantaneous  green-up  originally 
proposed  would  not  properly  reflect  the  transition  from 
high  fire  danger  in  early  spring  to  low  fire  danger  in 
summer.  The  green-up  period  length  varies  Hnearly  from 
7  days  for  climate  class  1  to  28  days  for  climate  class  4. 
The  length  of  green-up  was  scaled  to  climate  class 
because  plants  growing  in  drier  cHmates  typically  re- 
spond more  quickly  to  favorable  growing  conditions  than 
plants  in  wetter  climates. 

When  a  spring  flush  of  growth  becomes  generally  ap- 
parent, the  user  specifies  the  beginning  of  green-up.  Her- 
baceous fuel  moisture  is  then  calculated  according  to  the 
equation: 

mch  =  mci  +  gu[(ah  +  bhXiooo)  -  mci]  (44) 
where 

ah  and  bh  are  climate-dependent  intercept  and  slope  for 
annuals  or  perennials  from  table  9, 
Xiooo  is  the  fuel  moisture  recovery  value  from  eq.  43, 
mci  is  the  1-hour  timelag  fuel  moisture,  and 
m  is  the  elapsed  fraction  of  the  green-up  period. 

If  a  second  green-up  occurs  during  the  growing  season, 
the  Xiooo  value  is  again  set  equal  to  the  1,000-hour 
timelag  fuel  moisture.  The  same  green-up  procedure  is 
followed,  except  herbaceous  moisture  increases  from  its 
current  value  instead  of  the  1-hour  timelag  fuel 
moisture. 

After  the  green-up  period  is  complete  (gu  =  1.0),  the 
herbaceous  fuel  moisture  is 

mCh  =  ah  +  bhXiooo  (y  =  a  +  bx).  (45) 


Table  9.  — Regression  slopes  and  intercepts  for  drying  rates  of  live  fuels  (from  Burgan  1979) 


Climate  class 


Inter  Inter-  Inter-  Inter- 

Fuel  Slope  cept       Slope  cept      Slope  cept       Slope  cept 


Woody  fuels  7.5  12.5  8.2  -5.0  8.9  -22.5  9.8  -45.0 

Herbaceous  fuels 
Annuals  and  perennials  - 

green  stage  12.8  -70.0  14.0  -100.0  15.5  -137.5  17.4  -185.0 

Annuals  -  transition  stage  18.4  -150.5  19.6  -187.7  22.0  -245.2  24.3  -305.2 

Perennials  -  transition  stage  7.4  11.2  8.3  -10.3  9.8  -42.7  12.2  -93.5 
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Figure  7.— Comparison  of  X1000  1,000-hour  values  in  live  fuel 
nnoisture  nnodel  (from  Burgan  1979,  p.  8). 

When  the  user  specifies  that  the  herbaceous  vegetation 
has  cured  phenologically,  or  frozen,  herbaceous  fuel 
moisture  again  equals  the  1-hour  timelag  fuel  moisture. 
A  freeze  occurs  at  or  below  a  minimum  temperature  of 
25°  F,  or  when  minimum  temperatures  between  32°  and 
26°  F  occur  five  times. 

Figure  8  illustrates  the  relationship  between  the  Xiqoo 
value  and  the  annual  and  perennial  herbaceous 
moistures. 


SEPT. 


AUG. 

Summer  Season  — 


Freeze 


Figure  8. 
moisture 


—An  example  season  trace  of  X-iooo  ''^e  herbaceous  fuel 
values  (from  Burgan  1979,  p.  9). 
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TRANSFER  OF  FUEL  BETWEEN  LIVE  AND  DEAD 
CATEGORIES 

Before  green-up  the  entire  herbaceous  fuel  load  falls 
under  the  1-hour  class.  During  green-up  the  live  her- 
baceous fuel  load  transfers  from  the  1-hour  to  the  live 
fuel  category  as  herbaceous  moisture  increases  from  30 
to  120  percent.  At  levels  greater  than  120  percent  the 
herbaceous  fuel  load  reaches  its  maximum,  and  the 
1-hour  timelag  fuel  load  its  minimum. 

As  herbaceous  plant  moisture  decreases  later  in  the 
growing  season  and  varies  between  30  and  120  percent, 
the  load  of  perennial  herbaceous  fuels  is  shifted  between 
the  live  and  dead  fuel  categories.  This  indicates  a  transi- 
tion stage.  One  hundred  and  twenty  percent  approx- 
imates the  upper  limit  for  transition,  roughly  defining 
the  moisture  content  at  which  new  growth  is  complete 
and  foliage  mature.  Thirty  percent  was  chosen  as  the 
lower  limit  because  it  approximates  the  fiber  saturation 
point,  below  which  herbaceous  plants  are  assumed  dead. 

The  process  differs  slightly  for  annual  herbaceous 
plants.  Following  green-up  the  model  does  not  allow  the 
moisture  content  of  annuals  to  increase.  The  fuel  load 
for  annuals  transfers  from  the  live  category  to  the  dead 
category,  never  in  the  reverse  direction,  as  allowed  with 


perennials.  At  30  percent  moisture  content  (after 
phenological  curing  or  following  a  freeze  in  the  fall),  all 
the  herbaceous  fuels  have  been  added  back  to  the  1-hour 
timelag  class. 
The  fuel  load  transfer  equations  are: 

Wihd  =  Wih  +  Whf  (46) 
Whg  =  Whd  -  f)  (47) 
where 

f  =  -O.Olllmch  +  1.33  and  0  <  f  <  1.0 
and 

W^hd  is  the  total  load  of  1-hour  timelag  fuel,  including 
dead  herbaceous  fuel  transferred  to  the  1-hour  timelag 
category, 

is  the  load  of  the  1-hour  timelag  fuel  before  inclu- 
sion of  any  cured  herbaceous  material, 
Wh  is  the  total  load  of  herbaceous  fuel, 
f  is  the  fraction  of  the  herbaceous  fuel  that  is  to  be 
transferred  to  the  1-hour  timelag  class, 

Whg  is  the  load  of  herbaceous  fuel  that  is  still  green, 
and 

mch  is  the  herbaceous  fuel  moisture. 
Figure  9  shows  herbaceous  fuel  load  changes  in  rela- 
tion to  moisture  content  changes. 


LIVE  FUEL  MODEL 

HERBACEOUS  FUEL  LOAD  CHANGES  IN  RaATION 
TO  PLANT  MOISTURE  CONTENT  CHANGES 


1-hTLFM  250  120  30        1-h  TL  FM 

HERBACEOUS  FUEL  MOISTURE  (PERCENT) 


Figure  9.— Stylized  live  fuel  moisture  model  load  transfer  function 
(from  Burgan  1979,  p.  11). 
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Figures  10  and  11  give  examples  of  fuel  transfer  be- 
tween live  and  dead  categories  for  perennial  and  annual 
grasses.  Herbaceous  moisture  will  typically  approach 


250  percent  at  completion  of  green-up,  but  given  a  dry 
spring  or  late  green-up,  herbaceous  moistures  peak  at 
lower  values. 
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Herbaceous  moisture  greater  than  30%,  live  fuel  load  increasing 
Figure  10.— Typical  seasonal  trace  of  live  fuel  moistures  and  resulting 
1-hour  load  transfer  between  dead  and  live  categories  for  perennial  her- 
baceous fuels  (from  Burgan  1979,  p.  12). 

2.4 


2.0 


Freeze,  live 
load  zero  again 


1.6 


o 

t  1.2 


0.4 


ANNUAL  HERBACEOUS 
MOISTURE  CONTENT 


MAY 
Pre-areen ,  live 


fuel  load  zero 


Green  ^' 


4-h  TL  LOAD 


,  HERBACEOUS  LOAD 
CURED 


 A.  I  I  I  I  


150 


125 


100 


75 


50 


25 


o 

UJ 

< 


AUG.  SEPT. 
Herbaceous  moisture  less  than  30%, 


live  fuel  load  zero 
■  Herbaceous  moisture  less  than  120%,  live  fuel  load  variable 
Herbaceous  moisture  greater  than  120%,  live  fuel  load  maximum 
Herbaceous  moisture  greater  than  30%,  live  fuel  load  increasing 

Figure  11.— Typical  seasonal  trace  of  live  fuel  moistures  and  resulting 
1-hour  load  transfer  between  dead  and  live  categories  for  annual  her- 
baceous fuels  (from  Burgan  1979,  p.  13). 
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Woody  Fuel  Moisture  Model 

Prior  to  spring  green-up,  woody  shrubs  are  assumed 
dormant,  and  the  woody  fuel  moisture  (mCw)  held  con- 
stant. Measurements  of  chamise  leaf  moisture  in 
southern  California  indicate  minimum  values  for  woody 
plants  in  climate  class  2  is  about  60  percent  (Dell  and 
Philpot  1965).  Similarly,  measurements  in  the  Southwest 
suggest  a  minimum  woody  moisture  content  of  70  per- 
cent for  climate  class  3  (Blackmarr  and  Planner  1968). 
Pregreen  mc^  values  are  defined  as  50,  60,  70,  or  80  for 
climate  classes  1,  2,  3,  or  4,  respectively.  During  spring 
green-up,  woody  moisture  gradually  increases  from  the 
pregreen  minimum: 

mcw  =  mcwo  +  gu[(aw  +  b^mciooo)  "  mc^o]  (48) 
where 

mCwo  is  the  pregreen  minimum  moisture, 
a^  and  b^  are  climate-dependent  intercept  and  slope 
from  table  9, 

mciooo  is  the  1,000-hour  timelag  fuel  moisture,  and 
gu  is  the  fraction  of  the  green-up  period  that  has 
elapsed. 

If  a  second  green-up  occurs  during  a  growing  season, 
woody  fuel  moisture  increases  from  its  current  value  in- 
stead of  from  the  pregreen  value. 

After  green-up  completion  (gu  =  1)  woody  fuel 
moisture  is 

mcw  =       +  b^mciooo.  (49) 

When  the  shrubs  become  dormant,  woody  fuel  moisture 
is  set  back  to  the  minimum  value  specified  by  the 
climate  class.  Figure  12  illustrates  a  typical  woody  fuel 
moisture  profile  and  its  relationship  to  1,000-hour 
timelag  fuel  moisture. 


FIRE  BEHAVIOR  MODEL 

The  Rothermel  (1972)  mathematical  fire  model  pro- 
vided an  excellent  method  for  systematically  integrating 
fuels,  weather,  and  topographic  information  to  determine 
potential  fire  behavior,  hence  fire  danger,  in  a 
homogeneous  fuel  bed. 

The  fire  behavior  portion  of  the  NFDRS  consists  of 
the  spread  component  (SC)  and  the  energy  release  com- 
ponent (EEC).  Using  relationships  developed  by  Byram 
(1959)  and  Albini  (1976b),  these  components  are  com- 
bined in  the  burning  index  (BI).  The  BI  reflects  the 
potential  containment  problem  presented  by  a  single  fire. 
It  is  linearly  related  to  the  predicted  flame  length  at  the 
head  of  a  fire. 

The  fire  model  predicts  fire  spread  and  intensity  in 
homogeneous  fuel,  based  solely  on  static  fuel  and  en- 
vironmental properties  assessed  prior  to  ignition,  not 
dynamic  fire  variables  such  as  temperature,  heat  fluxes, 
and  induced  winds. 

By  using  weighting  factors  assigned  to  the  various 
fuel  sizes,  the  fire  model  can  accept  heterogeneous  fuel 
arrays  composed  of  various-sized  fuel  particles.  By 
weighting  various  parameters,  attributes  of  hetero- 
geneous fuel  arrays  are  characterized  for  the  fuel  bed. 
Thus,  while  a  particular  fuel  bed  may  have  four  dead 
fuel  components,  each  with  its  own  surface  area-to- 
volume  ratio,  it  has  a  singular  "characteristic"  surface 
area-to- volume  ratio  and  "characteristic"  fuel  moisture, 
dependent  on  the  weighting  factors.  The  assumption  of 
horizontal  and  vertical  fuel  continuity  (and  thus  a 
uniform  packing  ratio  and  bulk  density)  throughout  the 
bed  is  maintained. 
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Fire  Model  Parameter  Weighting 

Fuel  arrays  consist  of  particles  of  various  sizes,  some 
live  and  some  dead.  A  method  of  weighting  surface  area- 
to-volume  ratios  (a)  and  fuel  moisture  content  (mf)  was 
developed  to  synthesize  these  varied  parameters  into  a 
single,  characteristic  value  for  use  in  the  fire  model. 

Rothermel  (1972)  introduced  surface-area  weighting  for 
calculation  of  rate  of  spread.  Weighting  fuels  by  surface 
area  eliminated  making  arbitrary  decisions  about  which 
fuel  sizes  to  include  and  which  not  to  include  as  fire 
model  inputs.  A  singular  characteristic  parameter  is 
determined  by  weighting  parameters  in  a  fuel  array  com- 
posed of  a  mixture  of  particle  sizes. 

The  concept  of  a  unit  fuel  cell  facilitates  understand- 
ing of  fuel  size-class  distribution.  A  unit  fuel  cell  is  the 
smallest  volume  of  fuel  within  a  stratum  of  mean  depth 
that  has  sufficient  fuel  to  statistically  represent  fuel  in 
the  entire  complex.  Now  define: 

At  as  the  mean  total  surface  area  of  fuel  per  unit  fuel 
cell, 

Aj  as  the  mean  total  surface  area  of  fuel  of  i^h 
category  perunit  fuel  cell,  and 

Ajj  is  the  mean  total  surface  area  of  fuel  of  jth  size 
class  and  i^h  category  per  unit  fuel  cell.where 

i  is  the  fuel  category  (dead  =  1,  hving  =  2),  and 

j  is  the  fuel  size  class  (1  =  1-hour,  2  =  10-hour,  etc.). 

The  mean  total  surface  area  per  unit  fuel  cell  of  each 
size  class  with  each  category  is  determined  from  the 
mean  loading  of  that  size  class  and  its  surface-area-to- 
volume  ratio  and  particle  density: 

Aij  =  io  ijWij)/  Qp.  (50) 
The  mean  total  surface  area  of  the  i^h  category  per 
unit  fuel  cell  and  the  mean  total  surface  area  per  unit 
fuel  cell  are  then  obtained  by  summing  the  areas  within 
each  category  and  within  the  fuel  cell  with  equations  51 
and  52: 


Ai  =  ^  A 


At  =  E  Ai. 

1=1 


(51) 


(52) 


Two  weighting  parameters  are  then  calculated  and 
used  to  weight  fuel  particle,  fuel  bed,  and  moisture  con- 
tent parameters  throughout  the  fire  spread  model: 

fij  =  Aij/Ai  (53) 
fi  =  Aj/At  (54) 


FIRE  BEHAVIOR  COMPONENTS 
Spread  Component  (SC) 

In  the  1978  NFDRS,  SC  is  numerically  equivalent  to 
the  predicted  rate  of  spread,  rounded  to  the  nearest 
whole  number.  The  SC  algorithm  is  the  unmodified 
Rothermel  fire  spread  model,  with  the  fuel  elements 
weighted  by  surface  area: 

SC  =  kR  (55) 
where 

k  is  a  SC  scaling  factor  (1.0  min/ft),  and 
R  is  the  predicted  rate  of  spread  (ft/min). 

WIND  AND  SLOPE  IN  THE  1978  NFDRS 

Windspeed  ajid  slope  enter  into  the  Rothermel  rate-of- 
spread  equation  through  coefficients  functional  upon 
wind,  slope,  and  fuel  bed  properties.  Windspeed 
measured  at  the  20-foot  level  (10-minute  average)  is 
reduced  to  midflame  height  (as  required  for  rate-of- 
spread  computations)  through  a  wind  reduction  factor. 

Wind  Reduction  Factors 

Ratios  of  the  windspeed  at  roughly  the  midflame 
height  to  the  20-foot  standard  height  were  calculated  for 
grass,  shrub,  and  timbered  areas  assuming  a  standard 
logarithmic  wind  profile.  Using  variable  roughness 
lengths,  the  ratios  fell  within  the  range  of  0.4  to  0.6. 

As  such,  wind  reduction  factors  were  made  a  fuel 
model  parameter.  For  the  grass-type  fuel  models,  0.6  is 
used.  For  the  shrub  and  brush  fuel  models  0.5  is  used; 
and  for  the  timber  fuel  models  the  midflame  windspeed 
is  0.4  of  the  20-foot  windspeed.  The  wind  factor  {(j)^)  is 
then  calculated  as  described  by  Rothermel  (1972,  eqs. 
47-50). 

Slope  Classes 

The  slope  classes  in  the  1978  NFDRS  were  selected  so 
the  slope  coefficient  (0s)  would  double  from  one  class  to 
the  next  higher  class  (slope  class  5  has  16  times  the  ef- 
fect [24]  on  rate  of  spread  as  slope  class  1).  Using  a  90 
percent  slope  for  the  slope  class  5  midpoint,  the  slope 
factor  (4>s)  for  class  5  was  computed  from  Rothermel's 
(1972)  equation  51.  The  class  4      was  then  computed  by 
halving;  class  3  by  quartering;  class  2  by  taking  one- 
eighth;  and  class  1  by  taking  one-sixteenth  of  the  class 
5  0s  value. 

The  1978  NFDRS  slope  classification  scheme  is  sum- 
marized in  table  10. 

Table  10.— 1978  NFDRS  slope  classes 


NFDRS 
Slope 
class 


Slope 
range 


Effective 

class 
midpoint 


Slope 
coefficient 


-Percent- 


0-25 
26-40 
41-55 
56-75 

>75 


22.5 
31.8 
44.5 
63.6 
90.0 


0.267 
.533 
1.068 
2.134 
4.273 
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Energy  Release  Component  (ERC) 

When  the  1972  NFDRS  was  being  developed,  it  was 
recognized  that  more  than  rate  of  spread  of  a  potential 
fire  had  to  be  considered.  Unfortunately,  there  are  no 
models  as  definitive  as  the  Rothermel  model  to  quantify 
the  effects  of  the  condition  of  large  fuels  and  the  lowest 
layers  of  duff  and  litter  on  the  fire  behavior  and  fire 
management  problem.  There  is  no  way  that  a  fire-danger 
rating  system  would  be  acceptable  if  the  ratings  were 
solely  based  on  the  condition  of  fine  fuels,  which  is  the 
effect  of  surface-area  weighting  of  fuel  model 
parameters. 

The  authors  decided  that  the  spread  component 
calculations  should  be  totally  consistent  with  rate-of- 
spread  calculations  in  the  Rothermel  model;  rate  of 
spread  is  dominated  by  the  condition  of  fine  fuels.  The 
approach  selected  to  bring  larger  fuels  into  play  was 
straightforward:  use  the  fuel  energy  computations  of  the 
spread  model— the  reaction  intensity— but  base  the  in- 
fluences of  the  different  classes  of  fuel  on  their  contribu- 
tion to  the  total  fuel  load.  Specifically,  for  the  energy 
release  component,  the  characteristic  surface  area-to- 
volume  ratio  and  weighted  fuel  moisture  of  the  fuel  bed 
would  be  calculated  using  fuel  class  weighting  factors 
based  on  loading,  not  surface  area  as  in  the  spread  com- 
ponent. This  solution  had  no  experimental  basis. 

For  the  1978  revision  it  was  hoped  that  work  being 
done  by  Frank  Albini  with  large  fuel  burnout  would  pro- 
vide a  concrete  approach  to  the  fuel  energy  problem.  Un- 
fortunately, Albini  was  not  satisfied  with  the  results  of 
his  effort  and  recommended  that  the  1972  loading- 
weighted  approach  be  retained.  Rothermel  agreed  with 
Albini,  so  with  that  counsel  the  1972  procedure  was 
retained. 

Another  decision  made  during  the  development  of 
1972  NFDRS  was  to  change  the  form  of  the  moisture- 
damping  coefficient  used  in  calculating  the  ERC.  The 
curve  form  of  this  coefficient  as  developed  by  Rothermel 
is  S-shaped.  It  declines  rapidly  with  increasing  fuel 
moisture  at  low  fuel  moistures  and  at  moisture  contents 
near  the  moisture  of  extinction;  it  is  relatively  flat  in  the 
midrange  of  fuel  moistures.  Again,  a  judgment  was 
made:  to  develop  and  use  a  moisture-damping  function 
that  did  not  have  the  flattened  "shoulder."  The  function, 

T,e  =  1  -  2(mc)  +  1.5(mc)2  -  0.5(mc)3 

(dimensionless)  (56) 

where 

mc  is  the  ratio  of  the  load-weighted  characteristic  fuel 
bed  moisture  to  the  moisture  of  extinction  (a  damping 
coefficient  is  computed  for  both  the  dead  and  live  fuel 
components), 

was  developed  by  Fosberg  to  assure  a  continuing  in- 
crease in  the  ERC  as  fuel  moistures  in  the  midrange 
decreased.  This  function  is  used  only  for  ERC  computa- 
tions, not  spread  component  which  uses 

t;  =  1  -  2.59(mc)  +  5.11(mc)2  -  3.52(mc)3  (57) 

which  was  developed  by  Rothermel.  This  function  was 
retained  in  the  1978  NFDRS. 


There  was  one  significant  change  made  to  the  ERC  in 
the  1978  version  that  should  be  noted.  The  1972  version 
computed  ERC  from  a  scaled  loading- weighted  reaction 
intensity  (I^e'  Btu/ft2-min): 

ERC72  =  j(lRe)- 

In  the  1978  version  the  residence  time,  t^,  as  defined  by 
Anderson  (1969)  was  included  to  make  the  ERC  relatable 
to  the  total  energy  released  per  square  foot  during  the 
flaming  combustion  stage  (residence  time): 

ERC78  =  k(lRe)(rr)  (58) 

where 

k  is  a  scaling  factor  (0.04  ft2/Btu), 
I  Re  is  the  loading- weighted  reaction  intensity 
(Btu/ft2-min),  and 

Tr  =  384/a  (min)  (59) 

where  a  is  the  surface  area-weighted  characteristic  sur- 
face area-to-volume  ratio  of  the  fuel  bed.  This  change 
was  made  for  several  reasons: 

1.  In  the  1972  NFDRS     was  introduced  at  the  point 
where  the  burning  index  (BI)  was  calculated.  Since  is 
fuel\nodel  dependent,  a  separate  table  to  calculate  the 
BI  was  required  for  each  fuel  model.  Combining  Tj.  with 
I  Re  made  only  one  BI  table  necessary  in  the  manual 
version. 

2.  Available  heat  per  unit  area,  E: 

E  =  TrdRe)  (Btu/ft2)  (60) 

is  more  understandable  and  easier  to  relate  to  than  reac- 
tion intensity. 

3.  Heat  per  unit  area  rather  than  another  intensity- 
based  index  is  more  consistent  with  the  original  idea  of 
a  fuel  energy  phase  as  proposed  by  Keetch. 

Burning  Index  (BI) 

The  1978  burning  index  is  the  scaled  predicted  flame 
length,  as  calculated  from  the  loading-weighted  reaction 
intensity  (iRe),  the  surface  area-weighted  rate  of  spread 
(R),  and  the  surface  area-weighted  residence  time.  A 
flame  length  relationship  developed  by  Byram  (1959): 
Fl  =  0.45(1)0.46  (61) 

where  I  is  the  fireline  intensity  (Btu/ft2-min),  was 
modified  by  Albini  (1976a),  and  used  to  estimate  flame 
length  from  which  the  BI  is  computed.  The  equation  for 
flame  length  now  incorporates  the  spread  component  and 
available  energy  (E): 

Fl  =  j[(R/60)(lRe)(7r)p-46    (ft)  (62) 
Fl  =  j[(SC/60)(25(ERC))]0-46    (ft)  (63) 
where 

R  is  the  rate  of  spread  (ft/min), 

I  Re  is  the  load- weigh  ted  reaction  intensity 

(Btu/ft2-min), 
Tj.  is  the  flaming  residence  time  (min), 
1/60  is  a  unit  conversion  (1  min/60  s),  and 
j  is  a  coefficient  (0.45  ft2-sec/Btu). 

Consequently, 

BI  =  JiFl  (64) 

where  jj  is  the  BI  scaling  factor  (10/ft). 
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IGNITION  COMPONENT,  RISK,  AND 
OCCURRENCE  INDEXES 

The  occurrence  indexes,  in  conjunction  with  the  burn- 
ing index,  compute  the  cumulative  fire  load  index.  The 
two  fire  occurrence  indexes,  man-caused  (MCOI)  and 
lightning-caused  (LOT),  give  daily  projections  of  the 
number  of  reportable  man-caused  and  lightning-caused 
fires  per  million  acres  of  protected  land.  Occurrence  in- 
dex values  range  from  0  to  100  and  are  scaled  such  that 
a  value  of  100  indicates  an  expected  fire  density  of  10 
fires/million  acres. 

The  MCOI  is  a  function  of  the  ignition  component  (IC) 
and  the  man-caused  risk:  LOI  is  a  function  of  a  weighted 
IC  and  lightning-caused  risk.  Risk  factors  are  the 
average  expected  number  of  firebrands  on  a  given  day 
by  source  (lightning-  or  man-caused).  Determined  from 
empirical  models,  risk  factors  allow  users  to  compute 
scaling  factors  to  adjust  the  risk  models  to  specific 
regions  (see  Deeming  and  others  1977). 

The  IC  is  the  probability  that  a  reportable  fire  will 
result  from  a  firebrand. 

MCOI,  tenfold  the  expected  number  of  reportable  man- 
caused  fires  per  million  acres  on  a  given  day,  is 

MCOI  =  IC(Rmc)  (65) 

where  Rmc  is  related  to  the  total  number  expected  man- 
caused  firebrands  for  the  day.  Similarly,  the  LOI,  or  ten- 
fold the  number  of  reportable  lightning  fires  per  million 
acres  for  a  given  day  is  defined  by: 

LOI  =  IC(Rl)  (66) 

where  Rl  is  related  to  the  number  of  expected  lightning 
strokes. 

Common  to  both  MCOI  and  LOI,  development  of  the 
IC  will  be  reported  first,  followed  by  a  discussion  of  the 
1978  NFDRS  risk  models. 

Ignition  Component 

The  ignition  component  is  calculated  from  the  prob- 
ability of  ignition,  P(I),  the  day's  spread  component,  SC, 
and  the  maximum  probable  spread  component  SC^ax- 
The  SCmax  is  a  fuel  model  parameter,  calculated  under  a 
set  of  severe  burning  conditions  (appendix  B).  This  dif- 
fers from  the  1972  NFDRS  ignition  component,  which 
was  simply  equal  to  the  probability  of  ignition.  The  1972 
version  of  ignition  probability  used  a  moisture  of  extinc- 
tion of  30  percent. 

PROBABILITY  OF  IGNITION,  P(I) 

P(I)  is  the  probability  that  a  firebrand  will  start  a  fire 
(reportable  or  not)  after  landing  on  receptive  fuels.  This 
differs  from  the  IC,  which  incorporates  burning  condi- 
tions (via  the  SC)  to  estimate  the  probability  of  a 
firebrand  becoming  a  reportable  fire.  The  P(I)  predicts 
only  whether  the  firebrand  has  sufficient  energy  to  pro- 
duce a  successful  ignition  and  is  taken  from  an  office 
report  prepared  by  Mark  Schroeder  in  1969.  Its  develop- 
ment is  included  here  because  of  its  importance  in  the 
system's  occurrence  indexes  and  general  unavailability. 
It  is  presented  in  original  form  of  metric  units. 


The  first  step  in  determining  P(I)  is  calculating  the 
heat  required  to  bring  a  fine  fuel  particle  with  a  given 
mci  from  its  initial  temperature  to  ignition  temperature, 
the  heat  of  preignition  (Qig). 

Heat  of  preignition  (Qjg,  cal/g)  is  calculated  by  sum- 
ming the  following  quantities: 

A.  The  heat  required  to  raise  the  temperature  of  the 
dry  fuel  from  its  initial  temperature,  T,,,  to  its  ignition 
temperature,  Tj  (assumed  to  equal  320  °C), 

B.  The  heat  required  to  raise  the  moisture  contained 
in  the  fuel  from  its  initial  temperature  to  the  boiling 
point, 

C.  The  heat  of  desorption, 

D.  The  heat  required  to  vaporize  the  moisture,  and 

E.  The  heat  required  to  raise  the  temperature  of  water 
vapor  contained  in  the  fuel  voids  from  the  boiling  point 
to  ignition  temperature. 

The  quantities  are  referred  to  as  Qa,Qb'  Qc  Qd-  Qe- 
respectively. 

To  compute  Qa  for  a  gram  of  fuel,  the  specific  heat  of 
dry  fuel,  Cf,  is  multiplied  by  the  temperature  range 
Tj  —  Tq.  According  to  Stamm  (1964),  Cf  varies  with  the 
temperature: 

Cf  =  0.266  +  0.00116T  (67) 

where  T  is  the  average  temperature  between  the  ignition 
and  initial  temperatures.  Thus 

Cf  =  0.266  +  0.00116(320  +  To)/2 


=  0.4516  +0.00058T, 


(68) 


and 


Qa  =  (Ti  -  To)(0.4516  +  0.00058To).  (69) 

Calculating  Qb  requires  the  temperature  change  from 
To  to  100°  C,  multiplied  by  the  mass  of  water  and  the 
specific  heat  (1.0  for  water): 

Qb  =  mfdOO  -  To)  (cal/g)  (70) 

where  mf  is  the  moisture  content  of  the  fuel  (fraction). 

Heat  of  desorption,  Qc,  is  the  heat  required  to  separate 
the  bound  water  from  the  fibers,  and  equals  the  heat 
given  off  (heat  of  adsorption)  when  water  vapor  is  ad- 
sorbed. From  Stamm's  (1964)  figure  12  —  1,  the  following 
equation  may  be  approximated: 

Qc  =  280  exp(-15.1mf)    (cal/g).  (71) 

The  total  heat  of  desorption^  is  obtained  by  integrating 
from  mf  to  mf=0: 


Qc  =  280  I  exp(-15.1mf)dmf 


Qc  =  -18.54  exp(-15.1mf)| 


(72) 
(73) 


Qc  =  18.54(1  -  exp(-15.1mf))    (cal/g).  (74) 

The  heat  required  to  vaporize  the  moisture  is  the  heat 
of  vaporization  times  the  mass  of  water: 

Qd  =  640mf.  (75) 

Qe,  the  heat  required  to  raise  water  vapor  in  the  voids 
from  the  boiling  point  to  ignition  temperature,  was 
calculated  to  be  negligible  compared  to  Qa,  Qb-  Qc 
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Qd,  and  is  omitted.  Thus 

Qig  =  Qa  +  Qb  +  Qc  +  Qd.  or 

Qig  =  144.51  -  O.266T0  -  0.00058Tq2  -  Toiiif 

+  18.54(1  -  exp(-15.1mf))  +  640mf  (cal/g). 

(76) 

Schroeder  then  proposed  that  if  a  firebrand  lands  on 
receptive  fuels  the  probability  of  ignition  as  a  function 
of  Qig  should  be  the  product  of  the  probability  that  a 
firebrand  of  a  specific  size  will  cause  an  ignition  and  the 
probability  that  the  firebrand  will  be  that  size.  The  lat- 
ter probability  is  dependent  on  the  size  distribution  of 
firebrands,  but  this  information  was  lacking.  Schroeder 
used  the  findings  of  Blackmarr  (1972)  for  the  probability 
that  a  specific  size  firebrand  will  cause  an  ignition.  This 
turned  out  to  be  reverse  S-shaped  curves  (fig.  13)  for 
slash  pine  litter  at  different  moisture  contents.  To 
generalize  this  solution,  Schroeder  first  defined  a  critical 
moisture  content  at  P(I)  =  0  for  a  specific  firebrand, 
above  which  no  ignition  will  take  place.  Using  equation 


100 


MINIATURE  MATCH 
FIREBRAND  TYPE  1 

ProbIt  of  100  -  (Ignition  %): 
1.925  +  0.359  (moisture  %) 


76,  he  converted  moisture  content  to  Qjg  at  a  constant 
temperature.  Defining  Qf  as  the  heat  from  a  firebrand  at 
the  critical  moisture  content  and  assuming  that,  at  the 
critical  moisture  content,  Qf  was  equal  to  Qjg,  he 
transposed  the  curve  to  show  P(I)  as  a  function  of 

Qf  -  Qig- 

Schroeder  next  reasoned  that  because  firebrand  size 
distributions  were  unknown  one  might  deduce  the 
distribution  by  knowing  the  shape  of  the  P(I/Qig)  curve. 
A  clue  as  to  the  shape  of  this  curve  could  be  obtained 
from  previous  studies.  John  Keetch,  in  a  circular  letter 
dated  April  19,  1960,  reported  on  his  survey  of  ignition 
studies  used  in  fire-danger  rating  systems  in  various 
parts  of  the  country.  These  were  based  on  man-caused 
fire  occurrence  frequency  and  fine  fuel  moisture.  He 
found  considerable  agreement  among  the  studies.  Figure 
14  is  a  mean  curve  on  log  paper  for  these  studies,  with 
the  highest  value  of  ignition  probability  set  at  1.5  per- 
cent moisture  content.  Some  information  on  the  size 
distribution  of  firebrands  was  desired  to  be  contained  in 
the  index,  but  not  the  entire  spectrum  of  firebrand  sizes. 
Rather  only  those  that  caused  ignition  were  included. 
Lacking  any  information  as  to  where  the  truncation 
point  might  be,  Schroeder  arbitrarily  set  the  lowest 
value  of  fuel  moisture  used  at  the  50  percent  cumulative 
probability,  replotted  the  curve  from  figure  14  on  log- 
probability  paper,  and  found  that  a  straight  line  fit  the 
data  fairly  well.  At  the  same  time  he  converted  moisture 
content  to  Qjg  at  constant  temperature  using  equation  76. 

The  mathematical  problem  of  backtracking  from  pro- 
bability of  ignition  for  a  specific  firebrand  and  the  Qig 
curve  to  find  the  size  distribution  of  firebrands  turned 
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SINGLE  KITCHEN  MATCH 
FIREBRAND  TYPE  2 

Probit  of  100 -(Ignition  %)  = 
0.662-1-0.202  (moisture  %) 


10      15     20     25    30    35     40    45  50 


THREE  KITCHEN  MATCHES 
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0.781+0.141  (moisture  %) 
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Figure  13.— Probability  of  ignition  curves  as 
a  function  of  moisture  content  for  tliree 
firebrand  sizes  (from  Blackmarr  1972). 
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Figure  14.— Ignition  probability  as  a  function 
of  dead  fuel  moisture  content  (from 
Schroeder  1969). 
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out  to  be  intractable,  and  the  effort  was  abandoned.  In- 
stead, the  useful  portion  of  the  straight  line  function  on 
the  log-normal  plot  was  solved  numerically,  yielding: 


P(I/Qig)  =  0.000048x4-3/50 
where 

X  =  (400  -  Qig)/10. 
The  general  form  of  this  equation  is 


(77) 
(78) 


P(I/Qig)  =  (kgx  4)/50.0  (79) 
where 

X  =  (Qigmx  -  Qig)/10, 
k3  and  £4  are  empirical  constants, 
Qigmx  is  the  heat  or  preignition  for  a  fuel  particle  at 
the  extinction  moisture  content,  and 


ks,  k4,  and  Q 
extinction. 


igmx 


are  functions  of  the  dead  moisture  of 


Table  11  shows  some  computed  values  of  k^,  k4,  and 
Qigmx  foJ*  various  values  of  mxd-  The  1978  NFDRS  uses 
the  ka  and  k4  values  associated  with  mx  =  25  and  Qigmx 
=  344,  although  Schroeder  developed  the  P(I)  using  mx 
=  30  percent  and  Qigmx  —  380.  The  reason  for  this 
change  was  because  the  30  percent  used  in  the  1972 
NFDRS  was  felt  to  be  too  high.  An  mx  of  20  percent 
was  first  proposed,  but  a  value  of  25  percent  was  finally 
used  for  the  1978  version. 

The  P(I)  function  is  calibrated  with  other  constants,  Pi 
and  P2,  such  that  P(I)  =  0  when  mc^  =  mx  =  25,  and 
P(I)  =  100  when  mci  =  1.5  percent.  Table  11  shows 
other  possible  values  of  Pi  and  P2. 

Ignition  component  is  then  computed  from  a  function 
developed  by  Bill  Main  at  the  North  Central  Experiment 
Station.  Using  pooled  data  from  Jackson  (1977): 

IC  =  P(I)(0.1)(SCn)0-5  (80) 
where 

SCn  =  SC/SCmax'lOO  (maximum  SCn  =  100),  (81) 

P(I)  =  (P(I)'  -  0.00232)(100)/0.99767,  (82) 

P(I)'  =  (0.000923x3-66)/50,  (83) 

X  =  (344  -  Qig)/10,  (84) 


and 

SC  is  the  day's  spread  component, 

Qig  is  computed  from  equation  76, 

Pi  =  0.00232, 

P2  =  0.99767, 

kg  =  0.000923,  and 

k4  =  3.66. 

This  IC  is  used  directly  in  calculating  daily  values  of 
the  MCOI,  and  weighted  with  a  rain-influenced  IC  for 
calculating  the  daily  LOl  (section  titled  "Lightning  Risk 
Model  and  Lightning  Occurrence  Index"). 

Man-Caused  Risk  Model 

Jackson  (1977)  presented  a  model  for  predicting  the 
average  number  of  reportable  man-caused  fires  for  a 
given  ignition  component  value: 

E(number  of  fires)  =  IC(d)(di)(d2)(A)  (85) 

where  d,  di,  and  d2  are  empirically  derived  coefficients 
that  tune  the  model  to  a  particular  protection  unit,  and 
A  is  the  protection  unit  area  in  millions  of  acres. 
Jackson  developed  the  model  using  fire  occurrence  data 
from  Forest  Service  Region  9  (Eastern).  Data  from 
Region  3  (Southwest)  were  used  for  testing  the  model. 
The  reported  R2  value  was  0.89  for  the  tested  area. 

Equation  85  is  normalized  to  millions  of  acres  by 
dividing  by  A: 

E(number  of  fires/106  acres)  =  IC(d)(di)(d2).  (86) 

The  MCOI  is  scaled  such  that  10  fires  per  million  acres 
results  in  an  MCOI  value  of  100,  or: 

MCOI  =  10  E(number  of  fires/lOe  acres)  (87) 

Combining  equations  86  and  87  yields 

MCOI  =  10(IC)(d)(di)(d2)  (88) 
and  by  letting  do  =  lOd,  the  expression  becomes: 

MCOI  =  IC(do)(di)(d2).  (89) 
But  the  desired  final  form  of  the  equation  is 

MCOI  =  IC(Rmc/100)  (90) 
so  combining  equations  89  and  90  yields: 

MCOI  =  IC(do)(di)(d2)  =  (IC)(Rmc)/100.  (91) 


Table  11.— Normalization  constants  for  probability  of  ignition  as  a  function  of  moisture  of  extinction  and 
heat  of  Ignition 


Moisture  of 
extinction 
(mx,  %) 

Heat  of 
ignition 
(Qjge,  Btu) 

Constants 

k3 

k4 

Pi 

P2 

15 

268 

0.082164 

2.58 

0.01553 

0.98446 

20 

306 

.008379 

3.15 

.00554 

.99446 

25 

344 

.000932 

3.66 

.00232 

1.99767 

30 

381 

.000120 

4.11 

.00050 

.99950 

35 

425 

.0000124 

4.58 

.00072 

.99928 

40 

470 

.00000135 

5.02 

.00052 

.99948 

iUsed  in  the  1978  NFDRS  IC. 
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Solving  for  Rmc  yields: 

Rmc  =  100(do)(di)(d2).  (92) 

Coefficients  d,  do,  di,  and  d2  are  area  and  unit  (Region, 
Forest,  or  District)  dependent.  A  Regional  scaling  factor, 
d,  equals  the  historical  number  of  man-caused  fires  per 
million  acres,  per  day,  per  unit  of  ignition  component  for 
an  entire  region  over  a  5-year  base  period.  Expressed 
numerically,  it  equals  the  slope  of  the  regression  (forced 
through  zero)  of  the  number  of  fires  per  million  acre- 
days  on  the  daily  IC  for  the  parent  (regional)  unit.  An 
acre-day  is  the  acreage  of  the  regional  protection  unit 
multipHed  by  the  number  of  fire-weather  days  in  the 
5-year  analysis.  The  parent  protection  unit  should  be  at 
least  10  million  acres  in  size,  with  300  or  more  man- 
caused  fires  per  year.  Recalling  that: 

do  =  lOd 

do  will  generally  range  between  0.05  and  1.00. 

di,  a  man-caused  risk  scaling  factor,  accounts  for  dif- 
ferences in  risk  between  the  parent  unit  (Region)  and  the 
protection  unit  (Forest  or  District).  It  serves  as  an  input 
to  both  the  AFFIRMS  and  FIRDAT  processors. 

di  is  the  ratio  of  the  total  number  of  fires  on  the  pro- 
tection unit  divided  by  the  product  of  the  average  IC 
value,  area  of  the  protection  unit  (10^  acres),  and  the 
number  of  fire-weather  days,  all  divided  by  the  total 
number  of  fires  on  the  parent  unit  divided  by  the  prod- 
uct of  the  average  IC  value  and  the  total  acres  in  the 
parent  unit  and  the  total  number  of  fire-weather  days  in 
the  parent  unit. 


R 


SRij 


di  = 


^  (fireSu)/ICu  n^) 
E  (fireSr)/(IC,  A,  n^) 


(93) 


where 

subscripts  u  and  r  denote  the  protection  unit  and 
parent  region,  respectively, 

(fires)  is  the  total  number  of  man-caused  fires  for  the 
computational  period, 

IC  is  the  average  daily  IC  (man-caused)  for  the  period, 

A  is  the  area  in  millions  of  acres,  and 

n  is  the  total  days  used  in  computing  the  IC  values. 

A  correction  factor,  d2,  flows  from  a  subjective  evalua- 
tion of  the  activity  level  of  the  principal  risk  sources  of 
man-caused  fires  on  a  protection  unit.  d2  is  partitioned 
by  day  of  the  week  and  the  eight  statistical  fire  causes 
available  from  standard  fire  report  forms.  d2  adjusts  for 
short-term  fire  problems  (such  as  arson)  or  predictable 
changes  in  man-caused  risk  (weekend,  holidays,  hunting 
season,  etc.): 


d2  =  dj/25 

where  dj  is  the  unnormalized  R^c  defined  by: 

i=l 


(94) 


(95) 


where  Rij  is  the  partial  man-caused  risk  of  source  i  (of 
the  eight  standard  statistical  fire  causes)  on  day  j  of  the 
week,  and 


(96) 


8 

=  7  E 

i=l 


E  firesjj/firesu 

j  =  i 


(97) 


and  Gij  is  the  risk  (number  of  firebrands)  associated  with 
the  daily  activity  level  assigned  by  the  fire  manager  to 
risk  source  i  on  day  j. 

Values  of  G  for  daily  activity  levels  are  presented  in 
the  following  tabulation: 


Daily  activity 
level 

Extreme 

High 

Normal 

Low 

None 


100 
50 
25 
12 
0 


where  Rsrij  is  the  risk  source  ratio  for  day  j  of  the  week. 


Deeming  and  others  (1977)  and  Burgan  and  others 
(1977)  detail  computation  of  d,  do,  d^,  and  d2,  burying 
most  of  the  arithmetic  in  nomograms. 

Man-caused  risk  (Rmc)  determined  from  nomograms  is 
used  to  calculate  the  daily  MCOI  from  equation  90: 

MCOI  =  (IC)Rmc/100-  (98) 

The  average  number  of  fires  for  the  protection  unit 
then  becomes 

E(number  of  fires)  =  (MCOI)A(miUions  of  acres)/10. 

(99) 

APPLYING  MAN-CAUSED  RISK 

The  severity  of  man-caused  fire  problems  on  a  protec- 
tion unit  dictates  the  level  of  man-caused  risk  assess- 
ment. For  example,  a  complete  analysis,  with  eight  risk 
sources  and  a  separate  set  of  daily  risk  source  ratios, 
Rjj,  for  each  month  may  be  needed  for  high  man-caused 
fire  areas  such  as  southern  California  or  southern 
Georgia.  At  the  other  extreme,  managers  of  areas  such 
as  White  Mountain  National  Forest  in  New  Hampshire, 
which  averages  fewer  than  10  fires  per  year,  gain  little 
by  partitioning  risk  source  ratios.  In  such  situations  aU 
risk  sources  may  be  combined  and  it  will  be  an  un- 
necessary complication  to  stratify  by  day  of  week, 
month,  or  even  season.  Between  the  two  extremes  lies 
an  approach  that  can  tailor  man-caused  risk  assessment 
to  match  the  magnitude  of  the  problem. 

Lightning  Risk  Model  and  Lightning 
Occurrence  Index 

The  hghtning  occurrence  index  (LOI)  is  an  average 
estimate  of  the  number  of  reportable  lightning-caused 
fires  per  milUon  acres  protected: 

LOI  =  ICwRl  (100) 
where  ICw  is  the  day's  weighted  lightning-caused  igni- 
tion component,  and  Rl  is  an  estimate  of  fire-starting 
lightning  strikes  in  a  rating  area. 

The  1978  NFDRS  lightning-caused  fire  occurrence  in- 
dex model  is  based  on  an  approach  described  by  Fuquay 
and  others  (1979).  The  development  of  the  LOI  departed 
from  Fuquay's  model  at  the  point  where  the  probability 
of  ignition  of  a  Ughtning-caused  fire  is  computed.  For 
the  sake  of  simpHcity,  the  1978  NFDRS  IC  was  used 
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rather  than  Fuquay's  P(I).  The  two  methods  were  com- 
pared, and  the  resulting  difference  was  judged  to  be 
within  the  range  of  uncertainty  of  the  Fuquay  model. 

Lightning-caused  risk  (Rl)  is  a  function  of  several 
meteorological  parameters  grounded  in  the  lightning  ac- 
tivity level  (LAL).  Rl  is  calibrated  to  a  locality  by  a 
lightning  risk  scaling  factor  (K),  which  compensates  for 
not  using  the  Fuquay  lightning  ignition  model.  The 
lightning  risk  scaling  factor  is  empirically  derived  from 
archived  hghtning-caused  fire  occurrence  and  lightning 
activity  level  data.  The  lightning  risk  scaling  factor  is 
required  because  of  the  differences  between  local  fuel 
types  in  their  susceptibility  to  ignition  by  lightning. 


Thunderstorms,  amounts  of  rain,  and  lightning  discharge 
characteristics  also  differ  among  regions.  The  Fuquay 
lightning  ignition  model  was  derived  using  thunderstorm 
data  from  northwestern  Montana  and  northern  Idaho. 

The  LAL  is  an  index  (forecasted,  observed,  and 
verified)  of  thunderstorm  activity  on  a  scale  from  1  to  6 
for  a  2,500-mi2  rating  zone.  Table  12  tabulates  typical 
thunderstorm  attributes  for  the  six  LAL's.  Four  typical 
storm  characteristics  have  been  designated  for  each  LAL 
from  2  through  5.  Lightning  level  1  signifies  no 
thunderstorms,  and  LAL  6  indicates  a  special  condition 
involving  high-level  dry  thunderstorms. 


Table  12.— Thunderstorm  attributes  for  lightning  activity  levels  (adapted  from  Fuquay  and  others  1979) 


Fraction  of  area  covered  by  radar 
Rel.  freq.       echoes  of  indicated  strength 


Lightning 

activity  on  T/S  Very 

level         Clouds  and  storm  development      days  (%)     light       Light    Moderate   Heavy  0-T 


Percent  of  area  receiving  less 
than  the  amount  of  rain  indicated 


<0.1"     <0.3"  <0.9' 


Typical  Cloud  and  Precipitation  Conditions  (2,500  mi2  or  6  500  km2  area) 
1      No  thunderstorms'"  No  radar  echoes  No  precipitation 


Few  building  cumulus  only  occa- 
sionally reaching  cumulus  con- 
gestus  stage;  single  cumulonimbus 
in  forecast  area.  Visual  tops: 
<30,000  ft  (9  100  m)  m.s.l. 


10 


0.1 


<0.1 


90 


91 


100 


Scattered  cumulus  to  cumulus  con- 
gestus;  widely  scattered 
cumulonimbus  clouds;  cloud-to- 
ground  lightning  averaging  1-2  per 
min  max. 

Growing  cumulus  and  cumulus  con- 
gestus  stage  over  0.1-0.3  of  the 
area;  scattered  cloud-to-ground 
lightning  in  area  averaging  2-3  per 
min  max. 

Cumulus  congestus  common  over 
area,  occasionally  obscuring  the 
sky;  moderate  to  heavy  rain 
associated  with  cumulonimbus 
clouds  light  to  moderate  rain 
preceding  and  following  lightning 
activity.  Lightning  flashes  occurring 
steadily  at  some  place  in  or  during 
storm  period;  maximum  cloud-to- 
ground  flash  rate  greater  than  3  per 
min. 


35 


35 


.2  .2  0.05 


.2  .1  .05 


70  90  98  100 


65  80  95  100 


.3  .1         .05       0.02         50  75  85  100 


Scattered  towering  cumulus  with  a 
few  at  thunderstorm  stage;  very 
limited  horizontal  extent;  high 
bases  (15,000  to  17,000  ft  m.s.l.). 
Virga  in  most  prominent 
hydrometeor  form.  Lightning  flash 
rate  is  low,  averaging  less  than  1-3 
per  5-min  period  each  storm.2. 


<2 


Lighting 
activity 
level 


Maximum  radar  echo  height,  m.s.l. 


Feet 
<  28,000 
26,000-32,000 
30,000-36,000 
>  36,000 


Meters 
<  8,500 
7,900-9,700 
9,100-11,000 
>  1 1 ,000 


Lightning  -  Amount  and  Rate 

Cloud-to-ground  (CG) 
lightning  per  2,500  mi2 
(6  500  km2 


Occurance  rates,  maximum 


CG/5  min 


CG/15  min 


Ave.  rate/min 


20 
40 
80 
160 


0-10 
4-19 
9-32 


0-17 

6-32 
19-77 


1-  2 

2-  3 
3 


i|n  most  general  terms,  2  days  out  of  3  will  not  be  thunderstorm  days  during  a  typical  fire  season  in  the  mountainous  areas  of  the  western 
continental  United  States 
2Used  with  red-flag  warnings  of  extreme  fire  activity. 


29 


LAL  6  conditions,  although  rare  (fewer  than  2  percent 
of  thunderstorm  days),  often  present  extremely  severe 
fire  problems.  They  occur  when  sufficient  moisture  and 
instability  for  thunderstorm  development  are  found  only 
upwards  from  approximately  15,000  feet  (mean  sea 
level).  Virtually  no  wetting  precipitation  reaches  earth  (it 
evaporates— often  producing  strong  and  erratic 
downdrafts),  but  lightning  is  still  present.  When  an  LAL 
of  6  is  forecast,  a  red  flag  or  equivalent  alert  is  usually 
issued.  For  NFDR  purposes  Rl  and  LOI  are  both  set  at 
100. 

The  four  storm  characteristics  grounded  in  LAL's 
2  through  5  are: 

1.  Number  of  cloud-to-ground  lightning  discharges  per 
storm  (Scg), 

2.  Ground  area  covered  by  radar  echoes— rain  area  (Sg), 

3.  Area  intensity  of  rainfall,  and 

4.  Total  storm  size  (s^i). 

Forming  an  idealized  storm  (fig.  15),  these  and  other 
parameters  allow  calculation  of  total  cloud-to-ground 


lightning  strikes  over  a  rating  area.  The  racetrack  with 
an  infield  concept  pictures  a  large  Ughtning  area  encom- 
passing a  smaller  rain  and  lightning  area.  This  storm 
moves  according  to  upper  level  winds  (u^v)  and  rain  falls 
over  a  fraction  of  the  total  storm  area.  The  probability 
of  ignition,  P(I),  is  reduced  in  the  rainfall  area  due  to 
wetting  of  fine  fuels.  Table  13  depicts  typical  storm 
characteristics  by  Ughtning  activity  level.  The  adapta- 
tion of  the  model  to  the  1978  NFDRS  assumes  a  con- 
stant storm  speed  of  30  miles  per  hour  (Uw  =  30)  in  the 
computation  of  other  storm  attributes. 

Precipitation  duration  (p^)  is  computed  from  the  length 
of  the  rainfall  band  (Sg,  mi)  and  the  rate  of  storm  move- 
ment (Uw): 

Pd  =  Sg/Uw  (101) 

Storm  duration  (sd)  is  computed  from  an  empirically 
derived  function  dependent  on  cloud-to-ground  strikes 
(Scg): 

Sd  =  -86.83  +  153.41(Scg)0i437_  (^q2) 


LIGHTNING  AREA 
BOUNDARY 


Figure  15.  — Idealized  shape  of  lightning  and 
rain  areas  from  a  thunderstorm  (from  Fuquay 
and  others  1979). 


Table  13.— Storm  characteristics  by  lightning  activity  level 


Lightning  activity  level 


Attribute 

1 

2 

3 

4 

5 

6 

Cloud-to-ground  discharge 
lightning  strikes 
(strikes/2,500  mi2) 

0 

12.5 

25 

50 

100 

N/A 

Rain  corridor  length  (mi) 

0 

3 

4 

5 

7 

N/A 

Lightning  corridor  length  (mi) 

0 

7 

8 

9 

11 

N/A 

30 


Simple  geometry  defines  the  fractional  areas  with  (f^) 
and  without  (fg)  rainfall: 

f^  =  (4sdUwSs  +  7rs2)/(4sdUwSdi  +  7rs2di)  (103) 

fo  =  1  -  fw  (104) 

Fine  fuel  moisture  (1-hour  timelag)  within  the  rain  area 

is  adjusted  using  a  modified  1-hour  timelag  fuel  moisture 

diffusion  equation: 

mcir  =  mcl  +  (def  -  mci)(l  -  f  exp  (-T/t)).  (105) 
For  1-hour  fuel,  1  =  f  =  t;  thus 

mcir  =  mci  +  (d  -  mci)(l  -  exp(-T))  (106) 
where 

mcij.  is  the  fine  fuel  moisture  in  the  rain  area, 
mci  is  the  computed  1-hour  fuel  moisture, 
T  is  the  storm  precipitation  duration  (T  =  p^),  and 
d  =  (76.0  +  2.7pd)— see  10-hour  fuel  section. 

Heat  of  preignition  (eq.  76)  is  computed  for  the  wetted 
area  using  ambient  fuel  temperature  and  mcir  as  in- 
dependent variables. 

Probability  of  ignition,  P(I),  and  ignition  component 
for  the  rain  area  (ICi)  are  then  computed  from  equations 
80  through  84  as  in  the  MCOI. 

A  weighted  hghtning  ignition  component  (ICw)  is  then 
computed  using  the  rain  fractional  area  and  its 
associated  ignition  component  (ICj),  and  the  nonrain 
fractional  area  and  the  IC  computed  for  man-caused  oc- 
currence index  (IC): 

ICw  =  [fw(ICi)  +  fo(IC)]/100  (107) 

where  both  ICj  and  IC  are  upwardly  limited  at  100. 

If  a  general  rain  is  noted  on  the  fire- weather  observa- 
tion (or  forecast),  ICw  is  set  to  zero. 

LIGHTNING  STRIKES 

A  cloud-to-ground  (CG)  lightning  strike  consists  of 
several  sequential  events.  In  about  20  percent  of  these 
CG  flashes  an  extended  return  stroke,  or  a  "long  contin- 
uing current"  (LCC)  is  present.  Fuquay  and  others 
(1967,  1972)  have  shown  that  this  LCC  stroke  is  highly 
likely  to  ignite  most  lightning-caused  fires.  Therefore,  Rl 
is  modified  to  include  only  LCC  lightning  events: 


Rl  =  0.20(Seg)(RL)(Ki) 
where 


(108) 


Scg  is  the  total  number  of  lightning  strokes,  Rl  is  the 
local  scaling  factor,  and 

Ki  is  a  factor  used  to  scale  the  predictions  to  strikes 
per  million  acres  (model  was  developed  on  Nezperce  and 
Clearwater  National  Forests  in  the  Northern  Region  and 
Ki  =  1.5). 

A  second  factor  (K2)  is  used  to  scale  Rl  such  that 
when: 

LAL  =  5 
Rl  =  1 
Ki  =  1.5 
ICw  =  100 


Rl  will  equal  100.  Recalling  that  when  LAL  =  5, 
Scg  =  100,  and  including  K2,  equation  108  becomes 

100  =  K2(0.20)(100)(l)(Ki).  (109) 

Solving  for  K2  yields  K2  =  5/Kl  =  3.33.  Inserting  these 
back  into  the  Rl  equation: 

Rl  =  K2(0.20)(Scg)(K)  (liO) 

yields  the  number  of  LCC  discharges,  per  milHon  acres, 
scaled  to  a  local  protection  unit. 
The  lightning  occurrence  index  is  then  cEdculated  from: 

LOI  =  [IORl(ICw)]  +  0.25(LOIy)  (111) 
where 

the  multiplier  10  scales  LOI  to  100  when  10  fires  per 

million  acres  are  predicted, 
Rl  is  the  number  of  LCC  discharges  (locally  scaled), 
ICw  is  the  weighted  ignition  component,  and 
0.25(LOIy)  is  a  persistence  function  to  help  account  for 

carryover  from  the  previous  day's  LOI  (LOIy). 
The  lightning  risk  scaling  factor  is  empirically  derived 

as  follows: 

u=N 

10     J2   firesu  (112) 

u  =  l 

u=N 
1.6  LOIu 
u=l 

over  a  recent  3-  or  4-year  period.  The  K  is  then  adjusted 
every  5  years: 

u=N 

10  Koid     E   firesu  (113) 

■^new 


K  = 


u  =  N 

A    i;  LOIu 

u  =  l 

where  A  is  the  rating  area  in  milhons  of  acres,  and  fires 
and  LOI  are  as  defined  before  except  they  are  for  later 
periods  (see  Deeming  and  others  1977  for  more  complete 
operational  examples). 


THE  FIRE  LOAD  INDEX 

The  fire  load  index  measures  the  total  potential  con- 
tainment effort  that  may  be  needed  on  a  given  day.  It 
combines  the  containment  effort  for  a  single  fire  (burn- 
ing index),  with  the  average  number  of  expected  fires: 


FLI  =  V(BI2  +  (MCOI  +  LOI)2)/1.41  (114) 

where  the  BI,  and  the  sum  of  MCOI  and  LOI,  are 
limited  to  a  maximum  value  of  100.  The  normalizing  fac- 
tor of  1.41  limits  the  FLI  to  a  maximum  value  of  100. 
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VARIABLES  AND  SYMBOLS 


Variable  Description 

A  Area  (ft2,  m2) 

Ai  Age  correction  function  for  aging  fuel  sticks 

Aj  Mean  total  surface  area  of  ith  category  per 

unit  fuel  cell 

Ajj  Mean  total  surface  area  of  jth  class  and  ith 

category  per  unit  cell 

At  Mean  total  surface  area  per  unit  fuel  cell 

a  Age  of  fuel  moisture  sticks  (days) 

ah  Y-intercept  on  live  fuel  regression 

a^  Y-intercept  of  woody  fuel  moisture  regression 

B  Secondary  fuel  stick  aging  function 

Bj  Burning  index  (dimensionless) 

bh  Slope  on  live  fuel  regression 

bw  Slope  of  woody  fuel  moisture  regression 

C  Climate  class  dependent  fuel  stick  aging 

function 

Cc  1978  NFDRS  climate  class 

D  Daily  equilibrium  moisture  content  boundary 

conditions  for  100-hour  fuel  moisture  com- 
putations, % 

Dn  Period  boundary  moisture  conditions  for 

forecast  10-hour  fuel  moisture  or  1,000-hour 
computations,  % 

d  Regional  scaling  factor  (dimensionless) 

do  lOd 

di  Man-caused  risk  scaling  factor  (dimensionless) 

d2  Correction  for  risk  sources  (dimensionless) 

dj  Unnormalized  man-caused  risk  for  the  jth  day 

of  the  week  (number  of  person-caused 

firebrands/million  acres) 

EMC  Equilibrium  moisture  content,  % 

ERC  Energy  release  component  (dimensionless) 

fi  Surface  area  weighting  parameter  I  (dimen- 

sionless) ratio  of  surface  area  of  ith  category 
to  total  surface  area,  per  unit  fuel  cell 
Fl  Flame  length  (ft) 

FLI  Fire  load  index  (10  X  number  of  fires/miUion 

acres) 

Fq  Fourier  number  (dimensionless) 

fo  Fractional  area  of  thunderstorm  without 

rainfall 

f  Fraction  of  herbaceous  fuel  to  be  transferred 

to  1-hour  class 

fw  Fractional  area  of  thunderstorm  with  rainfall 

fij  Surface  area  weighting  parameter  II  (dimen- 

sionless) ratio  of  surface  area  of  jth  size  class 
to  total  surface  area  of  ith  category,  per  unit 
fuel  cell 

Gjj  Man-caused  risk  (number  of  firebrands) 

associated  with  the  daily  activity  level 
gu  Elapsed  fraction  of  green-up  period 

H  Heat  content  (Btu/lb) 

IC  Ignition  component  (dimensionless) 


Variable  Description 

ICw  Rainfall  area  weighted  ignition  component 

(dimensionless) 

Ir  Surface  area  weighted  reaction  intensity 

(Btu/ft2-min) 

I  Re  Load  weighted  reaction  intensity  (Btu/ft^-min) 

J  Julian  date 

K  Lightning  risk  scaling  factor  (dimensionless) 

Ki  Lightning  area  scaling  factor  I 

(dimensionless) 
K2  Lightning  area  scaling  factor  II 

(dimensionless) 

kj  Drying  or  wetting  factor  in  live  moisture 

model 

k2  Temperature  factor  in  live  moisture  model 

(dimensionless) 

ks  Ceilibration  constant  for  computing  P(I) 

(dimensionless) 
k4  Calibration  constant  for  computing  P(I) 

(dimensionless) 
L  Length  (ft,  m) 

LAL  Lightning  activity  level  (dimensionless) 

LOI  Lightning-caused  occurrence  index 

(number/million  acres) 

LOIy  Yesterday's  lightning-caused  occurrence  index 

(number/million  acres) 

MCOI  Man-caused  occurrence  index  (man-caused 

fires/ 106  acres) 

mf  Fuel  moisture  content  (fraction) 

mfw  Weighted  "fine"  dead  fuel  moisture  (fraction) 

mx  Moisture  of  extinction,  % 

mxd  Dead  fuel  moisture  of  extinction,  % 

mxi  Live  fuel  moisture  of  extinction,  % 

mc  Ratio  of  characteristic  fuel  bed  moisture  con- 

tent to  dead  or  live  fuel  particle  moisture  of 
extinction 

mch  Live  herbaceous  fuel  moisture,  % 

mCw  Woody  fuel  moisture  content,  % 

mci  1-hour  timelag  fuel  moisture  content,  % 

mcio  10-hour  timelag  fuel  moisture  content,  % 

mciok  Age-corrected  fuel  stick  moisture  content,  % 

mcioo  100-hour  timelag  fuel  moisture  content,  % 

mciooo  1,000-hour  timelag  fuel  moisture  content,  % 

P(I)  Probability  of  ignition 

Pi  Probability  coefficient  for  P(I)  calculations 

P2  Second  probability  coefficient  for  P(I) 

calculations 

Pa  Daily  precipitation  amount  (inches) 

Pd  Precipitation  duration  (h) 

Pdi  Precipitation  duration  for  first  forecast  period 

(h) 

Pd2  Precipitation  duration  for  second  forecast 

period  (h) 

Pr  Precipitation  rate  (inches/h) 
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VARIABLES  AND  SYMBOLS  (con.) 

Variable  Description 


Qig 

Qige 

Qigmx 

R 
Ri 

Rl 

Rmc 

RsRij 
r 

Scg 

Se 

St 

sc 

SCmax 

SC 

Sd 
Sdl 

Ss 
Si 

S2 

T 

Tf 
T 


N 


ig 

Vb 
Vf 

Vv 

W 

Wf 

Wh 


w 


hg 


W- 


Heat  of  preignition  (Btu/lb  fuel) 

Effective  heat  of  preignition  (Btu/lb  fuel) 

Minimum  heat  for  ignition  (Btu/lb  fuel) 

Rate  of  fire  spread  (ft/min,  cm/s) 

Number  of  fire-starting  lightning  strikes  in  a 

rating  class 

Lightning-caused  risk  (number  fire-starting 
strikes/million  acres) 

Man-caused  risk  (man-caused 
firebrands/million  acres) 

Risk  source  ratio  for  man-caused  risk  sources 
for  the  ith  source  on  the  jth  day 

Radius  (ft) 

Number  of  cloud-to-ground  lightning  strikes 
per  2,500  mi2 

Effective  mineral  content  (lb  silica-free 
minerals/lb  ovendry  fuel) 

Total  mineral  content  (lb  minerals/lb  ovendry 
fuel) 

Spread  component  (dimensionless) 

Maximum  probable  spread  component 
(dimensionless) 

Spread  component,  normalized  by  maximum 
probable  value  (dimensionless) 
Thunderstorm  duration  (h) 

Total  storm  length  (with  and  without  rain, 
mi) 

Rainfall  path  length  (mi) 

Lightning  area  scaling  factor  I 
(dimensionless) 

Lightning  area  scaling  factor  II 
(dimensionless) 

Moisture  stress  simulation  time  step  (h) 
Fuel  temperature  (°C) 
Ignition  temperature  (°C) 
Thunderstorm  movement  speed  (mi/h) 
Fuel  bed  total  volume  (ft^,  cm3) 
Fuel  bed  fuel  volume  (ft^,  cm3) 

Fuel  bed  void  volume  (ft^,  cm^) 
Dead-to-live  loading  ratio 

Weight  of  fuel  moisture  sticks  (g) 
Total  herbaceous  fuel  load  (Ib/ft^;  tons/acre) 
Load  of  herbaceous  fuel  that  is  still  green 
(lb/ft2;  tons/acre) 

Net  fuel  load  (for  live  fuels,  i=2  and  j  in- 
dicates 1-hour,  10-hour,  or  100-hour  timelag 
class;  for  dead  fuels,  i=l  and  j  indicates 
1-hour,  10-hour,  or  100-hour  timelag  class) 


Variable 

Wn 

Wo 

Wih 
Wihd 

^1000 

XylOOO 

overbar 

Greek 
Symbol 

q; 
«d 

/3 
X 

Xs 

8 

5d 

Am 

6^ 

Amciooo 


V 

o 

Qh 

ep 

o 

T 


Description 

Net  fuel  mass  (total  fuel  mass  less  inorganic 

mass,  lb/ft2;  tons/acre) 

Total  fuel  mass  (lb/ft2;  tons/acre) 

Load  of  1-h  prior  to  inclusion  of  herbaceous 
fuels  (lb) 

Total  1-hour  load  including  dead  herbaceous 
(lb/ft2;  tons/acre) 

Live  fuel  moisture  recovery  value,  % 

Yesterday's  live  fuel  moisture  recovery  value, 

% 

Mean  value  for  any  indicated  variable 


Description 

Fraction  of  living  fuel  to  dead  in  a  fuel  bed 

Angle  of  daylight  (degrees,  converted  to 
radians) 

Packing  ratio  (dimensionless) 
Probability  of  ignition  subfunction 
Fuel  particle  shape  factor 
Fuel  bed  depth  (ft,  m) 

Angle  of  solar  declination  (degrees,  converted 
to  radians) 

Potential  moisture  content  change  during 
stress  period,  % 

Actual  moisture  content  change  during  stress 
period,  % 

24-hour  change  in  the  1,000-hour  moisture 
content,  % 

Moisture-damping  coefficient  for  spread 
component 

Moisture-damping  coefficient  for  energy 
release  component 

Relative  moisture  content,  % 

Diffusivity  (cm^/s) 

Fuel  bed  porosity  (dimensionless) 

Fire-weather  station  latitude  (degrees) 

Slope  factor 

Wind  factor 

Fuel  bed  density  (Ib/ft^,  g/cm^) 

Fuel  particle  density  (Ib/ft^,  g/cm^) 

Surface  area-to- volume  ratio  (1/ft,  1/cm) 

Fuel  particle  timelag  (h) 

Flaming  residence  time  (min) 

Timelag  similarity  coefficient  (dimensionless) 


33 


REFERENCES 

Albini,  Frank  A.  Computer-based  models  of  wildland  fire 
behavior:  a  user's  manual.  Ogden,  UT:  U.S. 
Department  of  Agriculture,  Forest  Service, 
Intermountain  Forest  and  Range  Experiment  Station; 
1976a.  68  p. 

Albini,  Frank  A.  Estimating  wildfire  behavior  and  ef- 
fects. Gen.  Tech.  Rep.  INT-30.  Ogden,  UT:  U.S. 
Department  of  Agriculture,  Forest  Service, 
Intermountain  Forest  and  Range  Experiment  Station; 
1976b.  92  p. 

Albini,  F.  A.;  Baughman,  R.  G.  Estimating  windspeeds 
for  predicting  wildland  fire  behavior.  Res.  Pap. 
INT-221.  Ogden,  UT:  U.S.  Department  of  Agriculture, 
Forest  Service,  Intermountain  Forest  and  Range 
Experiment  Station;  1979.  12  p. 

Albini,  F.  A.;  Latham,  Don  J.;  Baughman,  R.  G. 
Estimating  upslope  convective  windspeeds  for  predic- 
ting wildland  fire  behavior.  Res.  Pap.  INT-257.  Ogden, 
UT:  U.S.  Department  of  Agriculture,  Forest  Service, 
Intermountain  Forest  and  Range  Experiment  Station; 
1982.  19  p. 

Anderson,  Hal  E.  Heat  transfer  and  fire  spread.  Res. 
Pap.  INT-69.  Ogden,  UT:  U.S.  Department  of 
Agriculture,  Forest  Service,  Intermountain  Forest  and 
Range  Experiment  Station;  1969.  20  p. 

Anderson,  Hal  E.;  Schuette,  Robert  D.;  Mutch,  Robert  W. 
Timelag  and  equilibrium  moisture  content  of 
ponderosa  pine  needles.  Res.  Pap.  INT-202.  Ogden, 
UT:  U.S.  Department  of  Agriculture,  Forest  Service, 
Intermountain  Forest  and  Range  Experiment  Station; 
1978.  28  p. 

Andrews,  Patricia  L.  Testing  the  fire  behavior  model.  In: 
Proceedings,  sixth  conference  fire  and  forest 
meteorology;  1980  April  22-24;  Seattle,  WA. 
Washington,  DC:  Society  of  American  Foresters;  1980: 
70-77. 

Bevins,  Collin.  An  evaluation  of  the  slash  fuel  model  of 
the  1972  National  Fire-Danger  Rating  System.  Seattle: 
University  of  Washington;  1976.  104  p.  Master's 
thesis. 

Blackmarr,  W.  H.  Moisture  content  influences  ignitabil- 
ity  of  slash  pine  litter.  Res.  Note  SE-173.  Asheville, 
NC:  U.S.  Department  of  Agriculture,  Forest  Service, 
Southeastern  Forest  Experiment  Station;  1972.  7  p. 

Blackmarr,  W.  H.;  Planner,  W.  B.  Seasonal  and  diurnal 
variation  in  moisture  content  of  six  species  of  pocosin 
shrubs.  Res.  Pap.  SE-33.  Asheville,  NC:  U.S.  Depart- 
ment of  Agriculture,  Forest  Service,  Southeastern 
Forest  Experiment  Station;  1968.  lip. 

Brackebusch,  Arthur  P.  Gain  and  loss  of  moisture  in 
large  forest  fuels.  Res.  Pap.  INT-173.  Ogden,  UT:  U.S. 
Department  of  Agriculture,  Forest  Service,  Intermoun- 
tain Forest  and  Range  Experiment  Station;  1975.  50  p. 

Brown,  James  K.  Physical  fuel  properties  of  ponderosa 
forest  floors  and  cheatgrass.  Res.  Pap.  INT-74.  Ogden, 
UT:  U.S.  Department  of  Agriculture,  Forest  Service, 
Intermountain  Forest  and  Range  Experiment  Station; 
1970a.  16  p. 

Brown,  James  K.  Ratios  of  surface  area  to  volume  for 
common  fine  fuels.  For.  Sci.  16:  101-105;  1970b. 


Brown,  James  K.  Vertical  distributions  of  fuel  in  spruce- 
fir  logging  slash.  Res.  Pap.  INT-81.  Ogden,  UT:  U.S. 
Department  of  Agriculture,  Forest  Service,  Intermoun- 
tain Forest  and  Range  Experiment  Station;  1970c.  9  p. 

Brown,  James  K.  Field  test  of  rate-of-fire  spread  model 
in  slash  fuels.  Res.  Pap.  INT- 116.  Ogden,  UT:  U.S. 
Department  of  Agriculture,  Forest  Service,  Intermoun- 
tain Forest  and  Range  Experiment  Station;  1972.  24  p. 

Burgan,  Robert  E.  The  effect  of  latitude  and  season  on 
index  values  in  the  1977  NFDR  system.  In:  Pro- 
ceedings, 4th  national  conference  on  fire  and  forest 
meteorology;  1976  November  16-18;  St.  Louis,  MO. 
Gen.  Tech.  Rep.  RM-32.  Fort  Collins,  CO:  U.S.  Depart- 
ment of  Agriculture,  Forest  Service,  Rocky  Mountain 
Forest  and  Range  Experiment  Station;  1976:  70-74. 

Burgan,  Robert  E.  Estimating  live  moisture  for  the  1978 
National  Fire-Danger  Rating  System.  Res.  Pap. 
INT-226.  Ogden,  UT:  U.S.  Department  of  Agriculture, 
Forest  Service,  Intermountain  Forest  and  Range 
Experiment  Station;  1979.  16  p. 

Burgan,  Robert  E.;  Cohen,  Jack  D.;  Deeming,  John  E. 
Manually  calculating  fire-danger  ratings— 1978  Na- 
tional Fire-Danger  Rating  System.  Gen.  Tech.  Rep. 
INT-40.  Ogden,  UT:  U.S.  Department  of  Agriculture, 
Forest  Service,  Intermountain  Forest  and  Range  Ex- 
periment Station;  1977.  49  p. 

Byram,  George  M.  Combustion  of  forest  fuels.  In:  Forest 
fire  control  and  use.  2d  ed.  New  York:  McGraw-Hill 
Book  Co.;  1959.  175  p. 

Byram,  George  M.  An  analysis  of  the  drying  process  in 
forest  fuel  material.  Unpublished  paper  presented  at 
the  1963  international  symposium  on  humidity  and 
moisture;  1963  May  20-23;  Washington,  DC.  38  p. 

Cohen,  Jack  D.;  Deeming,  John  E.  The  1978  National 
Fire-Danger  Rating  System— basic  equations.  Gen. 
Tech.  Rep.  PSW-.  Berkeley,  C A:  U.S.  Department  of 
Agriculture,  Forest  Service,  Pacific  Southwest  Forest 
and  Range  Experiment  Station.  In  preparation. 

Countryman,  Clive  M.;  Philpot,  C.  W.  Physical 
characteristics  of  chamise  as  a  wildland  fuel.  Res.  Pap. 
PSW-66.  Berkeley,  CA:  U.S.  Department  of 
Agriculture,  Forest  Service,  Pacific  Southwest  Forest 
and  Range  Experiment  Station;  1970.  16  p. 

Cramer,  Owen  P.  Predicting  moisture  contents  of  fuel 
moisture  indicator  sticks.  Res.  Pap.  PNW-41.  Portland, 
OR:  U.S.  Department  of  Agriculture,  Forest  Service, 
Pacific  Northwest  Forest  and  Range  Experiment 
Station;  1961.  17  p. 

Deeming,  J.  E.;  Lancaster,  J.  W.;  Fosberg,  Michael  A.; 
Furman,  William  R.;  Schroeder,  Mark  J.  The  National 
Fire-Danger  Rating  System.  Res.  Pap.  RM-84.  Fort 
Collins,  CO:  U.S.  Department  of  Agriculture,  Forest 
Service,  Rocky  Mountaiin  Forest  and  Range  Experi- 
ment Station;  1972,  revised  1974.  165  p. 

Deeming,  John  E.;  Burgan,  Robert  E.;  Cohen,  Jack  D. 
The  National  Fire-Danger  Rating  System— 1978.  Gen. 
Tech.  Rep.  INT-39.  Ogden,  UT:  U.S.  Department  of 
Agriculture,  Forest  Service,  Intermountain  Forest  and 
Range  Experiment  Station;  1977.  63  p. 

Dell,  John  D.;  Philpot,  Cheirles  W.  Variations  in  the 
moisture  content  of  several  fuel  size  components  of 
live  and  dead  chamise.  Res.  Note  PSW-83,  Berkeley, 


34 


CA:  U.S.  Department  of  Agriculture,  Forest  Service, 
Pacific  Southwest  Forest  and  Range  Experiment 
Station;  1965.  7  p. 
Fischer,  William  C;  Hardy,  Charles  E.  Fire- weather 
observers'  handbook.  Agric.  Handb.  494.  Washington, 
DC:  U.S.  Department  of  Agriculture,  Forest  Service; 
1976.  152  p. 

Fosberg,  Michael  A.  Drying  rates  of  heartwood  below 
fiber  saturation.  For.  Sci.  16:  57-63;  1970. 

Fosberg,  Michael  A.  Climatological  influences  on 
moisture  characteristics  of  dead  fuel:  theoretical 
analysis.  For.  Sci.  17:  64-72;  1971a. 

Fosberg,  Michael  A.  Moisture  calculations  for  the 
100-hour  timelag  fuels  for  fire-danger  rating.  Res.  Note 
RM-199.  Fort  Collins,  CO:  U.S.  Department  of 
Agriculture,  Forest  Service,  Rocky  Mountain  Forest 
and  Range  Experiment  Station;  1971b.  7  p. 

Fosberg,  Michael  A.  Theory  of  precipitation  effects  on 
dead  cylindrical  fuels.  For.  Sci.  18:  98-108;  1972. 

Fosberg,  Michael  A.;  Deeming,  John  E.  Derivation  of 
the  1-  and  10-hour  timelag  fuel  moisture  calculations  of 
fire-danger.  Res.  Note  RM-207.  Fort  Collins,  CO:  U.S. 
Department  of  Agriculture,  Forest  Service,  Rocky 
Mountain  Forest  and  Range  Experiment  Station;  1971. 
8  p. 

Fosberg,  Michael  A.;  Rothermel,  Richard  C;  Andrews, 
Patricia  L.  Moisture  content  calculations  for  1000-hour 
timelag  fuels.  For.  Sci.  27:  19-26;  1981. 

Fosberg,  Michael  A.;  Schroeder,  Mark  J.  Fine  her- 
baceous fuels  in  fire-danger  rating.  Res.  Note  RM-185. 
Fort  Collins,  CO:  U.S.  Department  of  Agriculture, 
Forest  Service,  Rocky  Mountain  Forest  and  Range 
Experiment  Station;  1971.  7  p. 

Frandsen,  William  H.  Fire  spread  through  porous  fuel 
through  the  conservation  of  energy.  Combustion  and 
Flame.  16:  9-16;  1971. 

Frandsen,  William  H.  Using  effective  heating  number  as 
a  weighting  factor  in  Rothermel' s  fire  spread  model. 
Gen.  Tech.  Rep.  INT-20.  Ogden,  UT:  U.S.  Department 
of  Agriculture,  Forest  Service,  Intermountain  Forest 
and  Range  Experiment  Station;  1973.  7  p. 

Fujioka,  Francis,  M.  Personal  communication.  Riverside, 
CA:  U.S.  Department  of  Agriculture,  Forest  Service, 
Pacific  Southwest  Forest  and  Range  Experiment 
Station;  n.d. 

Fuquay,  D.  M.;  Baughman,  R.  G.;  Latham,  D.  J.  A 
model  for  predicting  Hghtning  fire  ignition  in  wildland 
fuels.  Res.  Pap.  INT-217.  Ogden,  UT:  U.S. 
Department  of  Agriculture,  Forest  Service, 
Intermountain  Forest  and  Range  Experiment  Station; 
1979.  21  p. 

Fuquay,  D.  M.;  Taylor,  A.  R.;  Howe,  R.  G.;  Schmid,  C. 

W.,  Jr.  Lightning  discharges  that  caused  forest  fires. 

J.  Geophys.  Res.  77(12):  2156-2158;  1972. 
Fuquay,  Donald  M.;  Baughman,  R.  G.;  Taylor,  A.  R.; 

Howe,  R.  G.  Characteristics  of  seven  lightning 

discharges  that  caused  forest  fires.  J.  Geophys.  Res. 

72(24):  6371-6373;  1967. 
Furman,  R.  William.  Personal  communication. Fort 

Collins,  CO:  U.S.  Department  of  Agriculture,  Forest 

Service,  Rocky  Mountain  Forest  and  Range 

Experiment  Station;  n.d. 


Furman,  R.  William,  An  aid  to  streamlining  fire-weather 
station  networks.  Gen.  Tech.  Rep.  RM-17. 
Fort  Collins,  CO:  U.S.  Department  of  Agriculture, 
Forest  Service,  Rocky  Mountain  Forest  and  Range 
Experiment  Station;  1975.  4  p. 

Furman,  R.  William;  Brink,  Glen  E.  The  National  Fire- 
Weather  Library:  what  it  is  and  how  to  use  it.  Gen. 
Tech.  Rep.  RM-19.  Fort  Collins,  CO:  U.S.  Department 
of  Agriculture,  Forest  Service,  Rocky  Mountain  Forest 
and  Range  Experiment  Station;  1975.  8  p. 

Furman,  R.  William;  Helfman,  Robert  S.  A  computer 
program  for  processing  historic  fire-weather  data  for 
the  National  Fire-Danger  Rating  System.  Res.  Note 
RM-234.  Fort  Collins,  CO:  U.S.  Department  of 
Agriculture,  Forest  Service,  Rocky  Mountain  Forest 
and  Range  Experiment  Station;  1973.  12  p. 

Gisborne,  H.  T.  Measuring  forest  fire-danger  in  northern 
Idaho.  Misc.  Publ.  29.  Washington,  DC:  U.S.  Depart- 
ment of  Agriculture;  1928.  63  p. 

Haines,  Donald  A.;  Frost,  John  S.  Weather  effects  on 
fuel  moisture  sticks:  corrections  and  recommendations. 
Res.  Pap.  NC-154.  St.  Paul,  MN:  U.S.  Department  of 
Agriculture,  Forest  Service,  North  Central  Forest 
Experiment  Station;  1978.  7  p. 

Haines,  Donald  A.;  Main,  William  A.;  Frost,  John  S.; 
Simard,  Albert  J.  Fire-danger  rating  and  wildfire  oc- 
currence in  the  northeastern  United  States.  In 
preparation. 

Haltiner,  Jeff.  Environmental  influences  on  fuel  moisture 
sticks.  Fort  Collins,  CO:  Colorado  State  University; 
1975.  104  p.  M.S.  thesis. 

Harrington,  Michael  G.  Estimating  ponderosa  pine  fuel 
moisture  using  National  Fire-Danger  Rating  moisture 
values.  Res.  Pap.  RM-233.  Fort  Collins,  CO:  U.S. 
Department  of  Agriculture,  Forest  Service,  Rocky 
Mountain  Forest  and  Range  Experiment  Station;  1982. 
7  p. 

Helfman,  Robert  S.;  Deeming,  John  E.;  Straub,  Robert  J.; 
Furman,  R.  Wilham.  User's  guide  to  AFFIRMS: 
time-share  computerized  processing  for  fire-rating. 
Gen.  Tech.  Rep.  RM-15.  Fort  Collins,  CO:  U.S. 
Department  of  Agriculture,  Forest  Service,  Rocky 
Mountain  Forest  and  Range  Experiment  Station;  1975. 
107  p. 

Helfman,  Robert  S.;  Straub,  Robert  J.;  Deeming,  John  E. 
User's  guide  to  AFFIRMS:  time-share  computerized 
processing  for  fire-danger  rating.  Gen.  Tech.  Rep. 
INT-82.  Ogden,  UT:  U.S.  Department  of  Agriculture, 
Forest  Service,  Intermountain  Forest  and  Range 
Experiment  Station;  1975,  revised  1980.  150  p. 

Hough,  W.  A.;  Albini,  F.  A.  Predicting  fire  behavior  in 
palmetto-gallberry  fuel  complexes.  Res.  Pap.  SE-174. 
Asheville,  NC:  U.S.  Department  of  Agriculture,  Forest 
Service,  Southeastern  Forest  Experiment  Station; 
1978.  44  p. 

Jackson,  Craig.  Predicting  man-caused  fire  incidence: 
man-caused  risk  and  the  NFDRS.  UnpubHshed  paper 
presented  at  the  Smoky  Bear/Fire  Prevention 
Workshop;  1977  January  17-21;  Shreveport,  LA.  St. 
Paul,  MN:  U.S.  Department  of  Agriculture,  Forest 
Service,  North  Central  Forest  Experiment  Station. 
14  p. 


35 


Keetch,  John  J.  Circular  Letter  4400.  Washington,  DC: 
U.S.  Department  of  Agriculture,  Forest  Service;  1960 
April  19. 

Keetch,  John  J.;  Byram,  G.  M.  A  drought  index  for 
forest  fire  control.  Res.  Pap.  SE-38.  Asheville,  NC: 
U.S.  Department  of  Agriculture,  Forest  Service, 
Southeastern  Forest  Experiment  Station;  1968.  32  p. 

Kelsey,  Rick  G.;  Shafizadeh,  Fred;  Lowery,  David  P. 
Heat  content  of  bark,  twigs,  and  foliage  of  nine  species 
of  western  conifers.  Res.  Note  INT-261.  Ogden,  UT: 
U.S.  Department  of  Agriculture,  Forest  Service, 
Intermountain  Forest  and  Range  Experiment  Station; 
1979.  7  p. 

King,  Rudy  M.;  Furman,  R.  WiUiam.  Fire-danger  rating 
network  density.  Res.  Pap.  RM-177.  Fort  Collins,  CO: 
U.S.  Department  of  Agriculture,  Forest  Service,  Rocky 
Mountain  Forest  and  Range  Experiment  Station;  1976. 
4  p. 

Lancaster,James  W.  Timelag  useful  in  fire  danger  rating. 

Fire  Control  Notes.  31(3):  6-8,  10;  1970. 
Lettau,  H.  H.;  Davidson,  B.  Exploring  the  atmosphere's 

first  mile.  V.2.  New  York:  Pergamon  Press;  1957: 

377-578. 

Loomis,  Robert  M.;  Main,  William  A.  Comparing  jack 
pine  slash  and  forest  floor  moisture  contents  and 
National  Fire-Danger  Rating  System  predictions.  Res. 
Pap.  NC-189.  St.  Paul,  MN:  U.S.  Department  of 
Agriculture,  Forest  Service,  North  Central  Forest 
Experiment  Station;  1980.  10  p. 

Loomis,  Robert  M.;  Blank,  Richard  W.  Summer  moisture 
content  of  some  northern  lower  Michigan  understory 
plants.  Res.  Note  NC-263.  St.  Paul,  MN:  U.S. 
Department  of  Agriculture,  Forest  Service,  North 
Central  Forest  Experiment  Station;  1981.  4  p. 

Maiin,  William  A.;  Straub,  Robert  J.;  Paananen,  Donna  M. 
FIREFAMILY:  fire  planning  with  historic  weather 
data.  Gen.  Tech.  Rep.  NC-73.  St.  Paul,  MN:  U.S. 
Department  of  Agriculture,  Forest  Service,  North 
Central  Forest  and  Range  Experiment  Station;  1982. 
31  p. 

OttmEir,  Roger  D.  Evaluation  and  adjustment  of  large 
fuel  moisture  models  from  the  National  Fire-Danger 
Rating  System  for  prescribed  fire  planning.  Seattle: 
University  of  Washington;  1980.  126  p.  Master's 
thesis. 

Palmer,  William  C.  Meteorological  drought.  Res.  Pap.  45. 

Washington,  DC:  U.S.  Department  of  Commerce, 

Weather  Bureau;  1965.  58  p. 
Philpot,  C.  W.  Mineral  content  and  pyrolysis  of  selected 

plant  materials.  Res.  Note  INT-84.  Ogden,  UT:  U.S. 

Department  of  Agriculture,  Forest  Service, 

Intermountain  Forest  and  Range  Experiment  Station; 

1968.  4  p. 

Philpot,  C.  W.  Influence  of  mineral  content  on  the 
pyrolysis  of  plant  materials.  For.  Sci.  16:  461-471; 
1970. 

Rothermel,  Richard  C.  A  mathematical  model  for  predic- 
ting fire  spread  in  wildland  fuels.  Res.  Pap.  INT- 115. 
Ogden,  UT:  U.S.  Department  of  Agriculture,  Forest 
Service,  Intermount2dn  Forest  and  Range  Experiment 
Station;  1972.  40  p. 


Rothermel,  Richard  C.  How  to  predict  the  spread  and  in- 
tensity of  forest  and  range  fires.  Gen.  Tech.  Rep. 
INT- 143.  Ogden,  UT:  U.S.  Department  of  Agriculture, 
Forest  Service,  Intermountain  Forest  and  Range 
Experiment  Station;  1983.  161  p. 

Rothermel,  R.  C.  Proposed  live  fuel  moisture  model  for 
the  1978  National  Fire-Danger  Rating  System. 
Missoula,  MT:  U.S.  Department  of  Agriculture,  Forest 
Service,  Intermountain  Forest  and  Range  Experiment 
Station,  Northern  Forest  Fire  Laboratory.  Unpubhshed 
report. 

Rothermel,  R.  C.  Unpublished  data  on  file.  Missoula, 
MT:  U.S.  Department  of  Agriculture,  Forest  Service, 
Intermountain  Forest  and  Range  Experiment, 
Northern  Forest  Fire  Laboratory. 

Rothermel,  Richard  C;  Philpot,  Cheirles  W.  Predicting 
changes  in  chaparral  flammability.  J.  For.  71(10): 
640-643;  1973. 

Schroeder,  Mark  J.  Ignition  probability.  Fort  Collins, 
CO:  U.S.  Department  of  Agriculture,  Forest  Service, 
Rocky  Mountain  Forest  and  Range  Experiment 
Station;  1969.  Unpublished  report. 

Schroeder,  Mark  J.;  Fosberg,  Michael  A.;  Deeming,  John 
E.;  Lancaster,  James  W.;  Furman,  R.  William. 
Technical  development  of  the  National  Fire-Danger 
Rating  System.  Washington,  DC:  U.S.  Department  of 
Agriculture,  Forest  Service;  1973.  151  p.  Unpublished 
report. 

Simard,  Albert  J.  The  moisture  content  of  forest  fuels  - 
I.  Inf.  Rep.  FF-X-14.  Ottawa,  Ontario:  Forest  Fire 
Research  Institute;  1968.  47  p. 

Simard,  Albert  J.;  Eenignburg,  James;  Blank,  Richard. 
Predicting  fuel  moisture  in  jack  pine  slash.  St.  Paul, 
MN:  U.S.  Department  of  Agriculture,  Forest  Service, 
North  Central  Forest  Experiment  Station.  In 
preparation. 

Simard,  Albert  J.;  Main,  William  A.  Comparing  methods 
of  predicting  jack  pine  slash  moisture.  Can.  J.  For. 
Res.  12(4):  793-802;  1982. 

Sneeuwjagt,  Richard  J.  Evaluation  of  the  grass  fuel 
model  of  the  National  Fire-Danger  Rating  System. 
Seattle:  University  of  Washington;  1974.  104  p. 
Master's  thesis. 

Sneeuwjagt,  R.  J.;  Frandsen,  W.  H.  Behavior  of  ex- 
perimental grass  fires  vs.  predictions  based  on 
Rothermel's  fire  model.  Can.  J.  For.  Res.  7:  357-367; 
1977. 

Stamm,  A.  J.  Wood  and  cellulose  science.  New  York: 
Ronald  Press;  1964:  283-284. 

Thornthwaite,  C.  W.  The  climates  of  North  America  ac- 
cording to  a  new  classification.  Geogr.  Rev.  4: 
633-655;  1931. 

U.S.  Department  of  Agriculture,  Forest  Service,  Forest 
Products  Laboratory.  Wood  handbook:  wood  as  an 
engineering  material.  Agric.  Handb.  72.  Washington, 
DC:  U.S.  Department  of  Agriculture,  Forest  Service; 
1955,  revised  1974.  Paged  by  chapter. 

U.S.  Department  of  Agriculture,  Forest  Service.  The 
wildland  fire-danger  rating  system.  Berkeley,  CA: 
Pacific  Southwest  Forest  and  Range  Experiment 
Station;  1958,  revised  1968.  108  p. 


36 


U.S.  Department  of  Agriculture,  Forest  Service.  Hand- 
book on  National  Fire-Danger  Rating  System.  Hand- 
book 5109.11.  Washington,  DC:  U.S.  Department  of 
Agriculture;  1964. 

Van  Gelder,  Randall  J.  Unpublished  report.  Berkeley, 
CA:  U.S.  Department  of  Agriculture,  Forest  Service, 
Pacific  Southwest  Forest  and  Range  Experiment  Sta- 
tion; 1976. 

Van  Wagner,  C.  E.  Drying  rates  of  some  fine  forest 

fuels.  Fire  Control  Notes.  30(4):  5;  1969. 

Notes.  30(4):  5;  1969. 
Warren,  John  R.;  Vance,  Dale  L.  Remote  automatic 

weather  stations  for  resource  and  fire  management 

agencies.  G«n.  Tech.  Rep.  INT-116.  Ogden,  UT:  U.S. 

Department  of  Agriculture,  Forest  Service, 

Intermountain  Forest  and  Range  Experiment  Station; 

1981.  11  p. 

Yancik,  Richard  F.;  Roussopoulos,  Peter  J.  User's  guide 
to  the  National  Fire  Occurrence  Library.  Fort  Collins, 
CO:  U.S.  Department  of  Agriculture,  Forest  Service, 
Rocky  Mountain  Forest  and  Range  Experiment 
Station;  1982.  25  p. 


APPENDIX  A:  ORGANIZATION  AND 
LOCATIONS  OF  NFDRS  PROJECTS 

The  National  Fire-Danger  Rating  System  research 
work  unit  was  organized  in  1968  by  the  U.S. 
Department  of  Agriculture,  Forest  Service,  Rocky 
Mountain  Forest  and  Range  Experiment  Station,  Fort 
Collins,  Colo.  It  consisted  of  the  following  scientists: 
Mark  J.  Schroeder,  project  leader;  Michael  A.  Fosberg, 
research  meteorologist;  and  James  W.  Lancaster, 
forester. 

In  the  spring  of  1970,  Schroeder  returned  to  the  fire 
meteorology  project  at  the  Riverside  Fire  Laboratory, 
Riverside,  Cahf.  James  W.  Lancaster  became  project 
leader  and  John  E.  Deeming  joined  the  project, 
relocating  from  the  Southern  Forest  Fire  Laboratory  at 
Macon,  Ga. 

Other  persons  eventually  involved  included  R.  William 
Furman,  meteorologist;  David  Rainey,  forest  resource 
technician;  and  Bruce  McCammon,  forest  resource 
technician. 

In  1973,  Lancaster  and  Deeming  relocated  to  the  Boise 
Interagency  Fire  Center  (BIFC),  and  in  1974  were  joined 
by  Robert  J.  Straub,  computer  specialist. 

In  1975,  the  fire-danger  research  work  unit  was 
relocated  under  the  Intermountain  Forest  and  Range 
Experiment  Station  at  the  Northern  Forest  Fire 
Laboratory  in  Missoula,  Mont.  The  research  work  unit 
was  led  by  John  Deeming,  with  assistance  from  research 
foresters  Robert  E.  Burgan  and  Jack  D.  Cohen. 

The  1971  NFDRS  was  initially  tested  on  the  following 
National  Forest  System  units: 

Arizona:  Tonto,  Gila,  and  Coronado  National  Forests. 

California:  Klamath  National  Forest. 

Idaho:  Boise  National  Forest. 

Montana:  Lolo  National  Forest. 

New  Mexico:  Carson  National  Forest. 

Wyoming:  Medicine  Bow  National  Forest. 

Additionally,  it  was  tested  on  the  Chiricahua  and 
Saguaro  National  Monuments  in  Arizona;  the  Shoshone, 
Idaho,  District  of  the  Bureau  of  Land  Management;  and 
finally,  by  the  Georgia  State  Forestry  Commission. 
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APPENDIX  B:  FUEL  MODELS 

Sections  titled  "Forest  Fuels"  and  "Fuel  Models"  in 
the  main  text  defined  eight  fuel  parameters  required  to 
completely  describe  a  fuel  complex  for  fire  modeling  pur- 
poses. These  parameters  are: 

1.  H  -  heat  content,  Btu/lb, 

2.  TTp-  fuel  particle  density,  Ib/ft^, 

3.  St  -  total  mineral  content,  %, 

4.  Se  -  effective  mineral  content,  %, 

5.  Wo  -  total  fuel  load,  lb/ft2, 

6.  6  -  fuel  bed  depth,  ft, 

7.  o  -  surface  area  to  volume  ratio,  1/ft,  and 

8.  mx  -  moisture  of  extinction,  %. 

The  section  titled  "Spread  Component"  described  the 
use  of  wind  reduction  factors  used  to  reduce  the  20-foot 
windspeed  to  the  midflame  windspeed  required  for  solu- 


tion of  the  fire  spread  model,  and  the  section  titled 
"Ignition  Component"  discussed  the  maximum  probable 
spread  component  used  in  computing  the  ignition  compo- 
nent. For  calculating  SC^ax'  each  fuel  model  was  run 
through  the  fire  model  with  these  parameter  conditions: 


1-hour  fuel  moisture  4% 

10-hour  fuel  moisture  6% 

100-hour  fuel  moisture  8% 

1,000-hour  fuel  moisture  11% 

Herbaceous  fuel  moisture  65% 

Woody  fuel  moisture  75% 

Windspeed  (20-ft)  20  mi/h 

1978  NFDRS  slope  class  1  (22.5% 

slope) 

1978  NFDRS  climate  class  3 


Table  14  delineates  the  physical  attributes  of  the  fuel 
models  on  the  1978  NFDRS. 


Table  14.— Physical  attributes  of  each  of  the  fuel  models  on  the  1978  National  Fire-Danger  Rating  Systenn 


 Fuel  Model  

Attribute  ABCDEFGH  IJKLNOPQRSTU 


1-hour  dead 

0.2 

3.5 

0.4 

2.0 

1.5 

2.5 

2.5 

1.5 

12.0 

7.0 

2.5 

0.25 

1.5 

2.0 

1.0 

2.0 

0.5 

0.5 

10-hour  dead 

4.0 

1.0 

.5 

2.0 

2.0 

2.0 

1.0 

12.0 

7.0 

2.5 

1.5 

3.0 

1.0 

2.5 

.5 

.5 

100-hour  dead 

.5 

.25 

1.5 

5.0 

2.0 

10.0 

6.0 

2.0 

3.0 

.5 

2.0 

.5 

.5 

1,000-hour  dead 

12.0 

2.0 

12.0 

5.5 

2.5 

2.0 

1.0 

.5 

Woody 

11.5 

.5 

3.0 

.5 

9.0 

.5 

.5 

2.0 

7.0 

.5 

4.0 

.5 

.5 

Herbaceous 

.3 

.8 

.75 

.5 

.5 

.5 

.5 

.5 

.5 

.5 

.5 

Load  (tons/acre) 

1.0  1.5 

.5  1.5 

-  1.0 

2.5  .5 
.5  .5 

Surface-area-to-volume  ratio  (1/ft) 

1-hour  dead  3,000  700  2,000  1,250  2,000  700  2,000  2,000  1,500  1,500  1,500  2,000  1,600  1,500  1,750  1,500  1,500  2,500  2,500  1,750 
10-hour  dead  —     109    109    109    109    109    109    109    109    109    109      —     109    109    109    109    109    109    109  109 

100-hour  dead  —  30  —  —  30  30  30  30  30  30  30  —  —  30  30  30  30  30  —  30 
1,000-hour  dead  —  8  —  —  —  —  8  8  8  8  8  —  —  8  —  8  —  8  —  — 
Woody  —  1,250  1,500  1,500  1,500  1,250  1,500  1,500      —      _      _      _  1,500  1,500  1,500  1,200  1,500  1,200  1,500  1,500 

Herbaceous       3,000      —  2,500  1,500  2,000      —  2,000  2,000      _      _      _  2,000      —      —  2,000  1,500  2,000  1,500  2,000  2,000 

Heat  content  (all  fuels) 

(Btu/lb)  8,000  9,500  8,000  9,000  8,000  9,500  8,000  8,000  8,000  8,000  8,000  8,000  8,700  9,000  8,000  8,000  8,000  8,000  8,000  8,000 
Moisture  of  extinction  (%) 

Dead                     15      15      20      30      25      15      25      20      25      25      25      15      25      30      30      25  25      25  15  20 

Fuel  bed  depth  (ft)      .8     4.5     .75     2.0       .4     4.5     1.0       .3     2.0     1.3       .6     1.0     3.0     4.0       .4     3.0  .25       .4  1.25  .5 

SC^ax                     301      58      32      68      25      24      30        8      65      44      23     178     167      99      14      59  6      17  96  16 


Constant  fuel  particle  values  for  all  fuels: 

Fuel  particle  density  (Tr^):  32  Ib/ft^ 

Total  mineral  content  (S,):  0.0555 

Effective  mineral  content  (Se):  0.01 
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APPENDIX  C:  EQUILIBRIUM 
MOISTURE  CONTENT 

Equilibrium  moisture  content  (EMC)  is  an  important 
aspect  of  the  NFDRS  dead  fuel  moisture  models.  The 
EMC  is  the  moisture  content  (%)  of  a  fuel  particle  allow- 
ed sufficient  time  to  reach  equilibrium  with  its  environ- 
ment (no  net  moisture  exchange).  EMC  calculations  are 
functions  of  temperature  and  relative  humidity  and  are 
computed  from  regression  equations  developed  by 
Simard  (1968).  Simard  used  tables  from  the  Wood  Hand- 
book published  by  the  U.S.  Department  of  Agriculture, 
Forest  Service,  in  1955,  revised  1974,  as  the  basis  for  his 
equations.  The  equations  for  equilibrium  moisture  con- 
tent are: 

'  0.03299  +  0.281073h  -  0.000578hT 

h  <  11  (C-1) 

EMC  =/  2.22749  +  0.160107h-0.01478T 

10  <  h  <   51  (C-2) 


23.5  is  the  maximum  solar  decHnation  (at  the  solstices). 

The  distance  from  the  polar  axis  to  the  circle  of  il- 
lumunation  is  calculated  from 

Z  =  d(tan  6)  (D-3) 

where  d  is  the  distance  from  the  equator  to  the  latitude 
circle  defined  by  the  station  latitude  (0),  and  r  is  the 
radius  of  the  earth: 


d  =  r  sin  0 
S  =  r  cos  0. 

Thus 

Z  =  r  sin  0  tan  h 

and 


a.  =  COS 


-1 


(Z/S). 


Substituung  for  S  and  Z 
a  =  cos-1  r  sin  0  tan  a 

which  simplifies  to 

a.  =  cos~i  (tan  0  tan  5). 


(D-4) 
(D-5) 


(D-6) 


(D-7) 


(D-8) 


(D-9) 


21.06060  +  0.005565h2 

-0.00035hT  -  0.483199h        h  >  50  (C-3) 
where 

h  is  the  fuel-atmosphere  interface  relative  humidity 

(%), 

T  is  the  fuel-atmosphere  interface  temperature  ( °F), 
and 

EMC  is  the  fuel  particle  equilibrium  moisture  content 
(%). 

APPENDIX  D:  WEIGHTING  24-HOUR 
BOUNDARY  CONDITIONS  BY  DAY 
LENGTH 

Latitude,  date,  and  the  earth's  angle  of  decHnation  in 
its  orbit  about  the  sun  control  day  length  cycles.  In 
calculating  the  100-hour  and  1,000-hour  fuel  moisture, 
day  length  is  used  to  weight  the  24-hour  maximum  and 
minimum  EMC's  to  compute  the  boundary  EMC. 

The  day's  maximum  EMC  is  estimated  from  the 
minimum  temperature  and  maximum  relative  humidity 
(assumed  to  occur  simultaneously  at  night).  Minimum 
EMC  is  computed  from  the  day's  minimum  relative 
humidity  and  maximum  temperature  (assumed  to  occur 
simultaneously  during  the  day). 

These  values  are  then  weighted  by  the  hours  of  either 
nighttime  or  daytime  and  averaged  for  the  24-hour 
period: 

EMC24  =  (hdEMCinax  +  hiEMCn^n)/24  (D-1) 

where  h^  and  h^  are  hours  of  dark  and  hours  of  daylight, 
respectively.  Hours  of  daylight  are  computed  from  sta- 
tion latitude  and  date  of  observation. 

For  the  following  discussion,  refer  to  figure  16.  The 
day's  angle  of  solar  declination  (8)  is  computed  from: 

6  =  (23.5)sin(k(J-82))  (D-2) 

where 

J  =  Julian  date 

k  =  degree/day  constant  (0.9863  =  1.00),  and 


-EQUATOR 


Figure  16.— Geometry  for  computing  a  day's 
angle  of  solar  declination.  The  earth's  axis  is 
inclined  66.5°  from  the  plane  of  orbit. 
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a  is  then  proportional  to  day  length.  That  is 
hi  =  24h  -  2a24h 


(D-10) 


where  the  angle  measure  is  in  radians. 
Substituting  for  a,  day  length  is 

hi  =  24h  (l-(cos-i  (tan  0  tan  5))/x)  in  hours  (D-11) 

Using  the  data  from  stations  in  Lytle  Creek,  Calif. 
(34°  N),  Libby,  Mont.  (48°  N),  and  Fairbanks,  Alaska 
(65°  N),  Burgan  (1976)  tested  the  effects  of  daylight) 
weighted  boundary  conditions  on  the  MCOI,  BI,  and 
ERC  indexes  of  the  1978  NFDRS. 

Day  length  has  little  effect  on  MCOI  and  BI  because 
they  are  most  sensitive  to  the  fine  fuel  moisture  content 
as  expressed  through  the  spread  component  (used  in 
computation  of  both  the  MCOI  and  BI).  Here  daily 
values  affect  the  fine  fuel  moistures  more  than  seasonal 
trends. 

The  ERC,  on  the  other  hand,  is  strongly  influenced  by 
the  moisture  contents  of  the  heavier  100-hour  and 
1,000-hour  fuels,  which  are  affected  by  seasonal  drying 
trends.  Figures  17,  18,  and  19  illustrate  these 
differences. 
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Figure  18.— Energy  Release  Component  com- 
puted for  Libby,  Mont,  (latitude  48°  N,  1973 
data),  with  a  constant  day  length  value  (dot- 
ted lines)  and  with  variable  day  length  (solid 
line)  (from  Burgan  1976). 


Figure  17.  — Energy  Release  Component  com-  Figure  19.  — Energy  Release  Component  com 

puted  for  Lytle  Creek,  Calif,  (latitude  34°  N,  puted  for  Fairbanks,  Alaska  (latitude  65°  N, 

1974  data),  with  a  constant  day  length  value  1975  data)  with  a  constant  day  length  value 

(dotted  line)  and  with  variable  day  lengths  (dotted  lines)  and  with  variable  day  length 

(solid  lines)  (from  Burgan  1976).  (solid  lines)  (from  Burgan  1976). 
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APPENDIX  E:  REFLECTIONS  ON  THE 
DEVELOPMENT,  APPLICATION,  AND 
FUTURE  OF  THE  NATIONAL  FIRE- 
DANGER  RATING  SYSTEM 

John  E.  Deeming 

Pacific  Northwest  Forest  and  Range  Experiment  Station 
USDA  Forest  Service  Portland,  Oreg. 

Introduction 

Five  years  have  passed  since  the  NFDRS  was  updated 
(Deeming  and  others  1977).  I  remain  interested  in  and 
am  frequently  consulted  about  the  general  subject  of 
fire-danger  rating  and  the  NFDRS  specifically.  What 
follows  are  my  views  of  the  NFDRS— its  technical  ade- 
quacy, its  use,  and  the  obstacles  to  overcome  to  improve 
utilization  of  the  1978  or  any  future  version  of  the 
System. 

The  NFDRS  has  some  technical  flaws.  The  develop- 
ment team  was  working  against  a  deadline,  and 
therefore  utilized  existing  technology  or  the  technology 
that  could  be  developed  in  the  time  available.  Our  "best 
guess"  prevailed  when  all  else  failed.  I  am  convinced, 
however,  that  the  usefulness  of  the  System  is  not  com- 
promised by  technological  "splints",  but  by  inap- 
propriate application. 

System  Administration 

NFDRS  APPLICATIONS 

Fire  managers  expected,  and  in  most  instances  still  ex- 
pect, the  NFDRS  to  provide  the  high  spatial  and  tem- 
poral resolution  of  a  fire  behavior  prediction  system. 
Fire-danger  rating  areas  are  typically  greater  than 
100,000  acres  and  the  weather  is  observed  and  predicted 
for  one  specific  time  during  the  day  at  one  specific  loca- 
tion. Using  this  "worst  case"  approach  reduces  the 
uncertainty  concerning  the  meaning  of  a  rating  derived 
from  such  a  modest  sampling  effort. 

Many  decisions  made  at  the  different  administrative 
levels  of  fire  management  have  not  been  carefully 
analyzed  and  the  requirements  of  supporting  information 
have  not  been  developed.  Managers  must  recognize  that 
the  scale  of  the  information  used  in  making  a  decision 
must  match  the  scale  of  that  decision.  The  NFDRS  was 
designed  for  low  resolution,  medium-to-large-scale  ap- 
plications; the  Fire  Behavior  System  (FBS)  (Rothermel 
1983)  was  designed  for  high  resolution,  small-scale 
application. 

We  must  all  recognize  that  managing  wildfires  is  much 
more  complex  than  when  the  NFDRS  development  pro- 
gram began.  Managed  fires,  natural  prescribed  fires, 
delayed  initial  attack,  economic  efficiency,  and  cost-plus- 
net-value  change  were  not  everyday  considerations  in 
fire  control  offices  in  1968.  The  scope  of  the  modern  fire 
management  job  requires  new  approaches  for  using 
weather  and  fire-danger  ratings.  Techniques  for  blending 
the  NFDRS  and  FBS  have  not  been  developed.  What 
has  happened?  The  NFDRS  has  been  frequently  misap- 
plied and  has  been  used  where  it  has  failed  to  provide 


the  high  resolution  information  needed  by  management. 
Such  failures  have  stimulated  unjustified  criticism  and 
have  undermined  the  credibihty  of  the  System. 

PREPAREDNESS  CLASSES 

The  bottom  line  of  fire-danger  rating  in  the  day-to-day 
operation  of  a  fire  program  is  the  manning  class.  This  is 
sometimes  called  the  manning  and  specific  action  class, 
preparedness  class,  adjustive  class,  or  precautions  class. 
The  idea  is  to  divide  the  continuum  of  fire  danger  into 
discrete  classes  to  which  preplanned  management  ac- 
tions are  keyed.  The  designations  for  the  classes  are 
commonly  numerical,  I  (one)  through  V  (five);  or  adjec- 
tive. Low  through  Moderate,  High,  and  Very  High  to 
Extreme. 

The  NFDRS  does  not  include  the  procedures  for 
delineating  preparedness  classes;  the  NFDRS  produces 
numerical  indexes  of  the  fire  danger.  Manning  class 
definition  is  the  fire  manager's  job  and  one  that  has  not 
been  done  satisfactorily  in  most  instances. 

Manning  classes  should  be  based  on  incremental  needs 
for  specific  levels  of  suppression  action.  For  instance, 
using  wildfire  reports  (Yancik  and  Roussopoulos  1982) 
and  historical  fire  danger  (Furman  and  Brink  1975;  Main 
and  others  1982),  it  could  be  determined  at  what  level  of 
fire  danger  the  probability  of  failure  of  initial  attack 
with  hand  crews  is  unacceptably  high.  That  level  of  fire 
danger  might  then  be  designated  Class  III,  and  a  retar- 
dant  aircraft  preplanned  for  dispatch  on  Class  III  days. 
There  is  little  hope  the  NFDRS  or  any  updated  version 
of  the  NFDRS  will  be  satisfactory  in  the  eyes  of  fire 
managers  until  the  manning  class  problem  is  resolved. 

TRAINING 

The  NFDRS  is  a  flexible,  adaptable  system.  If  the 
user  understands  the  NFDRS  technology,  the  unique  re- 
quirements of  the  myriad  of  agencies  can  be  served  and 
the  evolving  range  of  fire  management  tasks  supported. 
The  situation  facing  management  is  that  a  continuing  ef- 
fort must  be  made  to  train  new  personnel  as  the  ranks 
of  experienced  NFDRS  users  thin  over  time. 

Training  is  a  continuing  requirement  for  the  proper  ap- 
plication of  NFDRS.  The  NFDRS  is  not  like  a  bell  buoy 
anchored  to  a  reef.  The  requirement  for  the  bell  buoy  to 
provide  visual  and  audio  warnings  of  a  navigation 
hazard  does  not  change;  the  requirements  for  the  infor- 
mation the  NFDRS  is  capable  of  producing  are  constant- 
ly evolving.  Hence,  up-to-date,  knowledgeable  users  are 
essential 

Technical 

INPUT  DATA 
Weather 

Without  a  program  that  assures  uninterrupted,  quality 
weather  observations,  the  System  will  not  work.  Fischer 
and  Hardy  (1976)  published  an  excellent  guide  for  set- 
ting up  and  maintaining  weather  stations  and  making 
observations.  If  this  guide  were  followed,  and  provisions 
made  for  taking  data  every  day  starting  30  days  before 
the  ratings  are  needed,  many  of  the  difficulties  with  the 
NFDRS  would  disappear. 
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A  network  of  weather  stations  designed  around  the 
NFDRS  would  decidedly  improve  predictions.  Until 
1980,  however,  such  a  system  was  not  feasible  because 
instruments  had  to  be  located  where  people  were 
available  to  read  and  record  data.  With  the  development 
of  solar-powered,  automatic  instrumentation,  and 
satellite  communications  (Warren  and  Vance  1981),  the 
"people"  limitation  no  longer  exists.  The  issue  of 
weather  station  network  design  is  now  timely  and  is 
being  pursued  at  both  Rocky  Mountain  and  Pacific 
Southwest  Forest  and  Range  Experiment  Stations, 
Forest  Sevice,  U.S.  Department  of  Agriculture  (Furman 
1975;  King  and  Furman  1976). 

The  higher  the  density  of  weather  stations,  the  more 
accurate  the  fire-danger  information.  Because  weather 
stations  are  expensive  to  set  up  and  operate,  the  design 
must  be  cost  effective.  What  information  is  vital  to  the 
decision  making  process?  How  accurate  must  the  fire- 
danger  information  be?  Network  design  guidelines  would 
enable  users  to  get  the  most  for  their  money.  This 
returns  us  to  a  subject  already  addressed,  the  require- 
ment for  fire  management  to  systematize  decisionmaking. 

Fuel  Models 

Fuel  model  selection  must  be  based  on  the  capabilities  of 
the  model  to  emulate  seasonal  trends  of  fire  potential  in 
a  fire-danger  rating  area.  The  suitability  of  a  fuel  model 
is  a  function  of  how  well  the  fuel  classes  and  the  relative 
proportions  of  those  classes  in  the  fuel  model  match 
what  is  on  the  ground.  If  it  can  be  shown  that  none  of 
the  20  available  fuel  models  are  suitable,  efforts  to  im- 
prove or  develop  new  models  should  be  considered. 

Before  building  new  fuel  models,  however,  one  should 
be  certain  that  is  where  the  fault  lies.  A  new  fuel  model 
will  not  eliminate  poor  performance  in  situations  where 
the  NFDRS  is  being  applied  incorrectly  or  there  exists 
an  inadequate  weather  observation  program.  A  great 
deal  of  art  is  involved  in  the  development  and  evaluation 
of  fuel  models.  Guidelines  and  a  training  program  are 
needed  for  fuel  model  developers. 

DERIVED  INPUTS 

Dead  Fuel  Moisture  Models 

The  unrealistic  indication  of  the  recovery  of  fire 
danger  after  a  precipitation  event  is  a  serious  short- 
coming of  the  NFDRS  in  areas  where  the  primary  fuel  is 
litter  and  duff.  Contrary  to  the  assumptions  made  in  the 
System,  the  atmosphere  is  not  the  only  moisture 
source/sink  (as  determined  by  measurement  of  the 
ambient  relative  humidity,  temperature,  and  precipita- 
tion). The  soil  and  duff,  after  a  precipitation  event,  are  a 
significant  source  of  moisture  for  dead  plant  material  in 
the  forest  floor. 

Another  factor  contributing  to  the  overrating  of  fire 
potential  in  forested  areas  is  the  way  the  atmospheric 
data  that  drive  the  dead  fuel  moisture  models  are  col- 
lected and  interpreted.  In  the  United  States,  fire  danger 
has  always  been  evaluated  at  open  sites:  the  practice  is 
consistent  with  the  "worst  case"  approach  to  rating  fire 
danger.  This  may  be  the  time  to  introduce  some  flexi- 
bility into  our  fire-danger  rating  policy  and  use  in-stand 
conditions  where  open  areas  are  atypical.  The  Canadians 
have  done  this  for  years  for  forested  areas. 


Evidence  that  the  current  dead  fuel  moisture  predic- 
tive models  can  be  improved  comes  from  Forest  Service 
research  units  at  East  Lansing,  Mich.  (Loomis  and  Main 
1980;  Simard  and  Main  1982;  Simard  and  others,  in 
preparation),  Seattle,  Wash.  (Ottmar  1980),  and  Tempe, 
Ariz.  (Harrington  1982).  Better  NFDRS  fuel  moisture 
models  would  greatly  benefit  prescribed  burning  and  the 
fire  behaviour  prediction  and  fire  management  planning 
systems.  The  NFDRS  fuel  moisture  predictions  are  used 
in  those  applications  because  the  NFDRS  is  the  only 
weather  "bookkeeping  system"  available.  Improved  dead 
fuel  moisture  models  would  be  easily  incorporated  into 
the  NFDRS,  but  one  must  be  careful  to  assess  the  "bot- 
tom line"  effects  on  the  ratings  (see  section  titled 
"System  Tuning"). 

Live  Fuel  Moisture  Models 

The  critical  step  in  applying  the  NFDRS  where  the 
condition  of  the  live  vegetation  dominates  the  fire 
damger  is  the  proper  keying  of  the  System's  live  fuel 
moisture  models  when  new  growth  commences  (Burgan 
1979).  It  is  a  stated  requirement  that  the  NFDRS  com- 
putations start  a  month  before  green-up,  but  that  it 
seldom  done  in  areas  with  a  mid-  to  late-summer  fire 
season. 

The  importance  of  responsive,  accurate,  live  fuel 
moisture  increases  as  the  ratio  of  live-to-dead  fuels  in 
the  fuel  complex  increases.  In  open  hardwood  and  con- 
ifer forests,  brush  and  chaparral,  and  range  vegetation 
types,  the  live-to-deal  fuel  ratio  is  high;  hence,  live  fuel 
moisture  predictions  must  at  least  parallel  actual  condi- 
tions. The  current  models  do  a  reasonable  job  if  green-up 
is  triggered  on  time  and  if  the  temperature  adheres  to  a 
typical  seasonal  pattern. 

What  happens  when  conditions,  principally  the 
temperature,  do  not  follow  the  typical  seasonal  patterns? 
The  models  do  not  work.  How  well  or  how  poorly  the 
current  NFDRS  live  fuel  models  work  has  not  been 
documented  except  for  a  study  conducted  in  the  North- 
east by  the  Forest  Service  research  group  at  East 
Lansing,  Mich.  (Loomis  and  Blank  1981).  It  is  my  view 
that  new  models  incorporating  temperature  are  required. 

The  NFDRS  provides  the  option  for  directly  entering 
live  fuel  moisture  values,  bypassing  the  models.  The 
flammability  of  chaparral  is  very  much  a  function  of  the 
moisture  content  of  the  foliage.  In  California,  therefore, 
for  years  fuel  moisture  has  been  directly  sampled. 
Because  a  general  model  for  predicting  the  moisture  con- 
tent of  live  plants  is  unlikely,  a  foUage-sampling  pro- 
gram may  be  warranted  where  the  condition  of  live 
plants  is  very  important  to  fire-danger  rating 
Drought 

The  effects  of  the  intermediate-term  (up  to  6  weeks) 
meteorological  drought  is  reflected  by  the  NFDRS  in 
both  the  live  and  dead  fuel  moisture  predictions.  Intro- 
duction of  the  live  fuel  model  to  account  for  curing  of 
live  herbaceous  vegetation,  and  the  1,000-hour  timelag 
class  for  large  dead  fuels  in  1978  ameliorated  a  short- 
coming the  1972  NFDRS,  but  it  did  not  eliminate  the 
drought  problem  in  those  areas  where  organic  soils  or 
very  deep  duff  and  litter  are  found.  For  those  regions 
provisions  must  be  made  to  incorporate  a  measure  of 
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long-term  meteorological  drought  such  as  that  developed 
by  Palmer  (1965)  or  Keetch  and  Byram  (1968). 
Midflame  Windspeed 

In  the  NFDRS,  the  factors  used  to  reduce  the  20-foot 
standard  windspeed  to  the  midflame  height  depend  on 
the  fuel  model  and  range  from  0.6  for  the  "open"  fuels 
to  0.4  for  the  "sheltered"  fuels.  Albini  and  Baughman 
(1979)  have  provided  a  better  set  of  reduction  factors 
that  range  from  0.1  to  0.6. 

A  great  deal  of  "noise"  in  the  fire-danger  ratings  is 
caused  by  the  literal  use  of  windspeeds  in  the  lower  (non- 
significant) speed  ranges.  Wind  varies  tremendously 
from  point  to  point  and  from  time  to  time.  Windspeed 
measured  at  one  point  at  one  time  is  not  necessarily  a 
reliable  indicator  of  the  wind  over  a  large  fire-danger 
rating  area. 

I  recommend  that  all  winds  less  than  10  mi/h  be 
treated  as  a  constant  6  mi/h.  Why  10  mi/h  for  the 
cutoff?  The  reasoning  is  that  winds  less  than  about  10 
mi/h  are  dominated  by  local  effects  such  as  differential 
heating  (Albini  and  others  1982).  Above  that  windspeed, 
the  odds  are  good  that  the  wind  field  is  being  dominated 
by  a  meso-  or  synoptic-scale  weather  process  that  can  be 
described  and  predicted. 

THE  FIRE-DANGER  RATING  PROCESSOR 
Moisture  of  Extinction 

The  moisture  of  extinction  is  the  fuel  moisture  level 
above  which  a  fire  will  not  spread.  In  the  NFDRS  the 
moisture  of  extinction  varies  by  fuel  model,  ranging 
from  a  low  of  15  percent  in  annual  grasses  to  40  percent 
in  southern  pine  litter.  Those  values  are  supported,  in 
some  cases  by  studies  (Brown  1972;  Sneeuwjagt  1974; 
Bevins  1976;  and  Sneeuwjagt  and  Frandsen  1977),  but 
they  have  been  established  subjectively  for  most  fuel 
models.  The  fact  is  that  the  moisture  of  extinction  not 
only  varies  between  fuel  types  but  also  within  fuel  types 
as  the  windspeed  and  slope  change.  Also,  the  moisture 
of  extinction  for  the  initial  igition  is  lower  than  the 
moisture  of  extinction  for  a  going  fire  because  of  the 
limited  energy  contained  in  the  typical  firebrand. 

Weighted  (Characteristic)  Fuel  Moisture 

In  the  Rothermel  fire  spread  model  (Rothermel  1972), 
the  heterogeneous  fuel  (more  than  one  fuel  class  in  a  fuel 
complex)  situation  is  addressed  by  using  a  weighting 
process  to  calculate,  for  instance,  a  characteristic 
surface-area-to-volume  ratio  for  the  fuel  complex.  In  the 
NFDRS,  the  weighting  for  the  Spread  Component  is 
done  exactly  as  Rothermel  (1972)  does,  using  as  the 
basis  of  weighting  the  proportions  of  the  total  surface 
area  of  all  fuels  contributed  by  the  individual  fuel 
classes.  For  the  Energy  Release  Component,  the 
weighting  is  based  on  the  contribution  of  the  individual 
class  load  to  the  total  load  of  the  fuel  complex. 
Weighting  by  surface  area  causes  fire  behavior  to  be 
underrated  when  the  larger  fuels  are  dry  enough  to  burn, 
whereas  weighting  by  loading  causes  the  fire  behavior  to 
be  underrated  when  the  larger  fuels  are  too  wet  to  burn. 
Weighting  by  loading  was  introduced  in  the  1972 
NFDRS  and  retained  in  the  1978  NFDRS  to  increase  the 
influence  of  the  condition  of  the  larger  fuels  on  the 


rating  in  the  upper  range  of  fire  danger. 

Using  the  weighted  fuel  moisture  causes  unrealistic 
ratings  when  the  100-hour  and  1,000-hour  timelag  fuels 
have  moisture  contents  well  above  the  fixed  dead  fuel 
moisture  of  extinction.  Put  another  way,  even  though 
the  fine  and  intermediate  fuels  are  dry  enough  to  burn 
and  carry  fire,  the  NFDRS  says  "zero"  because  the 
weighted  average  fuel  moisture  exceeds  the  moisture  of 
extinction.  Experienced  fire  managers  know  that  a  fire 
will  burn  through  one  fuel  stratum  without  involving  the 
other  fuel  strata.  The  NFDRS  should  reflect  that 
behavior. 

SYSTEM  TUNING 

If  changes  are  made  to  any  of  the  items  I  have  men- 
tioned, the  System  will  have  to  be  recahbrated.  The  1972 
and  1978  versions  of  the  NFDRS  were  subjectively 
tuned  to  compensate  for  limitations  remaining  in  the 
component  models.  How  is  that  done?  The  final  tuning 
is  done  by  manipulating  the  fuel  model  parameters  until 
the  ratings  properly  reflect  conditions  for  selected  cases 
(a  severe  fire  period  and  a  benign  fire  period,  for  in- 
stance) Recall  my  comment  that  fabricating  fuel  models 
involves  a  great  deal  of  art.  What  constitutes  a  proper 
rating  of  any  situation  is  very,  very  subjective.  We 
spent  hundreds  of  hours  with  users  adjusting  fuel  model 
parameters  to  cause  the  NFDRS  ratings  to  match  sub- 
jective, nonquantitative  appraisals  of  case  situations. 

SYSTEM  VALIDATION 

Currently,  there  is  no  common  measurement  of  any 
fire  phenomenon  to  correlate  with  NFDRS  ratings. 

Don  Haines  and  his  coworkers  at  East  Lansing  will 
soon  pubHsh  the  results  of  an  evaluation  using  (1)  prob- 
ability of  a  fire  day,  (2)  probability  of  a  large  fire  day  (a 
day  with  one  or  more  fires  10  acres  and  larger),  (3) 
number  of  fires  per  day  and  (4)  number  of  fires  per  fire 
day  (Haines  and  others,  in  preparation).  The  results  are 
encouraging  and  show  that  the  Ignition  Component  and 
the  Spread  Component  are  highly  correlated  with  those 
four  measures  of  fire  acitivity. 

Most  of  the  pubhshed  work  has  been  directed  at 
validating  components  of  the  System  such  as  the  fuel 
moisture  models  (Simard  and  Main  1982;  Simard  and 
others,  in  preparation;  Forsberg  and  other  1981; 
Harrington  1982;  Loomis  and  Main  1980;  Loomis  and 
Blank  1981;  Ottmar  1980).  Much  work  has  been  done  to 
evaluate  the  Rothermel  spread  model  which  is  the  basis 
of  the  NFDRS  Spread  Component  (Andrews  1980; 
Bevins  1976;  Brown  1972;  Sneeuwjagt  1974;  Sneeuwjagt 
and  Frandsen  1977). 

An  effort  is  under  way  at  the  Northern  Forest  Fire 
Laboratory  to  develop  software  to  access  the  Fire 
Report  Library  (Yancik  and  Roussopoulos  1981)  and  Fire 
Weather  Library  (Furman  and  Brink  1975)  for  com- 
parisons similar  to  those  done  at  East  Lansing.  The 
problem  remains,  however,  that  without  an  agreed-upon 
set  of  measures  of  fire  phenomena  and  standard  methods 
of  analysis,  it  is  impossible  to  determine  how  good  the 
System  is  or  how  much  proposed  changes  to  the  NFDRS 
would  improve  its  performance. 
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SUMMARY 

A  technical  revision  of  the  NFDRS  alone  will  not  cure 
all  the  problems  with  the  fire-danger  rating  programs. 
What  is  needed  is  a  national  program  that  will  em- 
phasize research  in  NFDRS  appHcation,  management, 
and  validation,  and  that  will  revise  the  System  as 
required. 

Piecemeal  or  regional  revisions  would  be  unwise 
because  every  modification  requires  an  extensive 
checkout  of  the  impact  on  System  performance.  It  would 
be  much  more  efficient  if  all  the  contemplated  changes 
were  made  at  once  by  one  group.  A  research,  develop- 
ment, and  applications  program  would  be  appropriate  for 
this  task. 

Training  at  the  national,  regional,  and  local  levels  is  a 
continuing  need  and  must  be  provided  for.  Fire  manage- 
ment is  constantly  evolving.  It  takes  knowlegeable,  ex- 
perienced people  to  redefine  the  role  of  NFDRS  and  to 
properly  apply  the  System. 

Fire  management  is  a  demanding  and  increasingly 
complex  job.  The  need  to  especially  good  decisionmaking 
is  increased  by  the  high  cost  of  wildfire  suppression  and 
the  forces  that  must  be  available  to  do  the  job.  If  pro- 
perly supported  and  implemented,  the  NFDRS  can  con- 
tribute to  efficient  fire  management. 
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The  Intermountain  Station,  headquartered  in  Ogden,  Utah,  is  one 
of  eight  regional  experiment  sfations  charged  with  providing  scien- 
tific knowledge  to  help  resource  managers  meet  hwman  needs  and 
protect  forest  and  range  ecosystems.  X'  > 

The  Intermquntajn-x^tation  included  the;  States  of  Montana, 
Idaho,'  Utah,^  Nev^ada,  arij  western  Wydr^ing."  About  231  million 
acres,  or  85  f^rc^t^  of  Jhe^nd  are^  'n  jt>e  Station  territory  are 
classified  as  fo*^#aQd -rat^eland:  These-,  lands  include  grass- 
lands, deserts,  shRiSiagds,  alpine  areas,  and  well-stocked  forests. 
They  supply  fiber  for  forest  industries;  minerals  for  energy  and  in- 
dustrial development;  and  water  for  domestic  and  industrial  con- 
sumption. They  also  provide  recreation  opportunities  for  millions 
of  visitors  each  year. 

Field  programs  and  research  work  units  of  the  Station  are  main- 
tained in: 

Boise,  Idaho 

Bozeman,  Montana  (in  cooperation  with  Montana  State 
University) 

Logan,  Utah  (in  cooperation  with  Utah  State  University) 

Missoula,  Montana  (in  cooperation  with  the  University 
of  Montana) 

Moscow,  Idaho  (in  cooperation  with  the  University  of 
Idaho) 

Provo,  Utah  (in  cooperation  with  Brigham  Young  Univer- 
sity) 

Reno,  Nevada  (in  cooperation  with  the  University  of 
Nevada) 


